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ABSTRACT 

Quantum algebras are a mathematical tool which provides us with a class of sym- 
metries wider than that of Lie algebras, which are contained in the former as a 
special case. After a self-contained introduction to the necessary mathematical 
I tools (g-numbers, g-analysis, g-oscillators, g-algebras), the su g (2) rotator model 

Q\ ■ and its extensions, the construction of deformed exactly soluble models (u(3)Dso(3) 

model, Interacting Boson Model, Moszkowski model), the 3-dimensional (/-deformed 
O { harmonic oscillator and its relation to the nuclear shell model, the use of de- 

QJ ' formed bosons in the description of pairing correlations, and the symmetries of 

the anisotropic quantum harmonic oscillator with rational ratios of frequencies, 
which underly the structure of superdcformed and hyperdeformed nuclei, are dis- 
cussed in some detail. A brief description of similar applications to the structure 
of molecules and of atomic clusters, as well as an outlook are also given. 
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1. Introduction 

Quantum algebras Bfflii (also called quantum groups) are deformed versions of the usual Lie al- 
gebras, to which they reduce when the deformation parameter q is set equal to unity. From the 
mathematical point of view they are Hopf algebras Q. Details about their mathematical properties can 
be found in Bat Their use in physics became popular with the introduction of the g-deformed 

harmonic oscillator (sec. 10) as a tool for providing a boson realization of the quantum algebra su 9 (2) 
(sec. 14), although similar mathematical structures had already been known (sec. 11). Initially used for 
solving the quantum Yang-Baxter equation i, quantum algebras have subsequently found anplications 
in several branches of physics, as, for example, in the description of spin chains0-i, anyons0O@, quan- 
tum optics^@§3lll§l, rotational and vibrational nuclear and molecular spectra, and in conformal field 
theories. By now several kinds of generalized deformed oscillators (sec. 12) and generalized deformed 
algebras (sec. 17) have been introduced. 

It is clear that quantum algebras provide us with a class of symmetries which is richer than the class 
of Lie symmetries, which are contained in the former as a special case. It is therefore conceivable that 
quantum algebras can turn out to be appropriate for describing symmetries of physical systems which 
are outside the realm of Lie algebras. 

Here we shall confine ourselves to applications of quantum algebras in nuclear structure physics. 
The structure of this review is as follows: In order to make this review self-contained, we are going 
first to give a brief account of the necessary tools: g-numbers and g-analysis (sees 2-9), g-deformed 
oscillators (sees 10-13), g-deformed algebras (sees 14-18). The remainder will be devoted to specific 
applications in nuclear structure problems, starting with phenomenology and advancing towards more 



microscopic subjects. The su 9 (2) rotator model (sees 19-23) and its extensions (sees 24-26), as well as 
the formulation of deformed exactly soluble models (u(3)Dso(3) model (sees 28, 29), Interacting Boson 
Model (sees 27, 30, 31), Moszkowski model (sec. 32)) will be covered in some detail. In addition the 
3-dimensional g-deformed harmonic oscillator and its relation to the nuclear shell model (sec. 28) will 
be extensively discussed. Subsequently, the use of quantum algebraic techniques for the description 
of pairing correlations in nuclei (sees 33-35), as well as the symmetries of the anisotropic quantum 
harmonic oscillators with rational ratios of frequencies (sec. 36) will also be considered in some detail. 
The latter are of current interest in connection with the symmetries underlying superdeformed and 
hyperdeformed nuclear bands (sec. 36). Finally, a brief account of applications of the same techniques 
to the structure of molecules (sec. 37) and of atomic clusters (sec^ 38), as well as an outlook (sec. 39) 
will be given. An earlier review of the subject has been given in cJ. 



2. q- numbers 

The q-number corresponding to the ordinary number x is defined as 

q x — q~ x 



x 



q-q 



-i 



(2.1) 



where q is a parameter. The same definition holds if x is an operator. We remark that g-numbers 
remain invariant under the substitution q — > g -1 . 

If q is real, g-numbers can easily be put in the form 



sinh(rx) 
sinh(r) 



(2.2) 



where q = e T and r is real. 

If g is a phase factor, g-numbers can be written as 



x 



sin tx 



suit 



(2.3) 



where g = e %T and r is real. 

In both cases it is clear that in the limit g — > 1 (or, equivalently, r — > 0) g-numbers (or operators) 
tend to the ordinary numbers (or operators): 



limlx] = x. 

A few examples of g-numbers are given here: 

[0] = 0, [1] = 1, [2]=g + g-\ [3] = g 2 + 1 + g~ 2 . 



(2.4) 



(2.5) 



Identities between g-numbers exist. They are, however, different from the familiar identities between 
usual numbers. As an exercise one can show (using the definition of g-numbers) that 



[a][b + l]-[b][a + l] = [a-b]. 
The q-factorial of an integer n is defined as 

[n]! = [n][n-l]...[2][l]. 
The q-binomial coefficients are defined as 

[ml! 



m 
n 



[m — n] ! [n] ! 



while the q-binomial expansion is given by 



fe=0 



m 
k 



a m - k (±by 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



In the limit q — > 1 we obviously have 



[nl! -> n! and m - ( m ), (2.10) 



where n! and (™) are the standard factorial and binomial coefficients respectively. 

It should be noticed that two-parameter deformed numbers have also been introduced 

M M = 5L -^r- (2.H) 

g — p 1 

In the special case p = q they reduce to the usual g-numbers. 

3. g-deformed elementary functions 

In addition to g-deformed numbers and operators, g-deformed elementary functions can be intro- 
duced. The q-exponential function is defined as 

00 a n 



e q {ax) ^ w . 



EOT^ (3.1) 



while the q-trigonometric functions are defined as 



x 2n+l oo x 2n 



rtBf (, ) = E(-ir^ rn5f , <«,(*) = E(-i)" M - 



It should also be noticed that g-deformed polynomials, such as g-deformed Hermite polynomials E3c3EZl 
and g-deformed Laguerre polynomials also exist (see also subsec. 36.3). 

The definitions of the g-exponential function and the g-trigonometric functions given above are not 
unique; for a different set of definitions, based on the Tsallis statistics, see 0. 

4. g-derivatives 

Proceeding along this path one can build a new differential calculus, based on g-deformed quantities 
(see for concise expositions). For this purpose the q-derivative is defined as 

D , (x) = /M-ZOrM (41) 

(g — g 1 )x 

The similarity between the present definition and the one of g-numbers (eq.(2.1)) is clear. 
Using the definition of the g-derivative one can easily see that 

D q x (ax n ) = a[n}x n -\ (4.2) 

D q x c q [ax) = ae q (ax). (4.3) 
^From the definition of the g-derivative one can further derive the sum rule 

D x (f(x) + g(x)) = D q J(x) + D*g(x), (4.4) 

as well as the rule 

D q X2 [a Xl ± bx 2 } m = ±[m]6[aa:i ± bx 2 ] mr ~ 1 , (4.5) 

where a and b are costants and [ax% ± bx2\ m is given by the g-binomial expansion (eq. (2.9)). One can 
also prove the q-integration by parts formula 

DlU{x)g{x)) = /M^W-Z^M^g-M ^ (46) 



From this, the following two forms of the Leibnitz rule can be derived 

Dl(f(x)g(x)) = (Dlf(x))g(q^x) + f(qx)(Dlg(x)), (4.7) 

Dl(f(x)g(x)) = (Dlgix^fiq-'x) + g(qx)(Dlf(x)). (4.8) 
In addition one can show the property 

D«J(qx) = qDlf(x)\ x=qx , (4.9) 

and the chain rules 

DlJ(x) = -Dlf(x) , (4.10) 
a 

Dlf(x n ) = [n]*"- 1 !* (4.11) 
where a is a constant. Another useful result is 

a — i 



-r n— 1 

Dff(x) = — J2 Dlf(q 2k -^x). (4.12) 



k=0 



n=0 



5. g-integration 

The ^-integration (see E1B for concise expositions) in the interval [0, a] is defined by 

/ f(x)d q x = a(q~ 1 -q)J2q 2n+1 f(q 2n+1 a), (5.1) 

„ n 

while for the interval [0, oo) one has 

/ f(x)d q x=(q- 1 -q) £ q 2n+1 f(q 2n+1 ). (5.2) 

n=-oo 

The indefinite q-integral is defined as 

„ oo 

/ f(x)d q x = (q- 1 -q)Y, q 2n+1 xf{q 2n+1 x) + constant, (5.3) 

J n=0 

where < q < 1. For entire functions /(#) one can easily see that this q- integral approaches the 
Riemann integral as q — > 1, and also that the operators of ^-differentiation and g- integration are inverse 
to each other 

Dl J f(x)d q x = f(x) = J Dlf{x)d q x. (5.4) 

One can also easily see that 

/ ax n ~ l d q x = -7-^ax n + constant, (5.5) 
J [n] 

J e q {ax)d q x = —e q {ax) + constant. (5.6) 
From (4.6) one can also prove the following formulae of integration by parts 

f(qx)(Dlg(x))d q x = f{x)g{x)\lzl ~ f {Dlf{x))g{q~ l x)d q x, (5.7) 

Jo 

* f(q- 1 x)(D^g(x))d q x = f(x)g(x)\ x x Z a ~ j* {Dlf{x))g{qx)d q x. (5.8) 
The following formulae can also be proved 

f{x)d aq x = a f{x)d q x, (5.9) 



f(x n )d gnX n = [n] I x n - L f(x n )d q x, 

i n— 1 



f(x)d q x 



J_ V- „2fc-(n-l) 



fc=0 



2fc-(n-l) 



The q-analogue for Euler's formula for the function T(x) is 

K 



e q (-x)x n d q x = [n\[n - l\[n - 2] . . . [1] = [n]!. 



(5.10) 
(5.11) 

(5.12) 



A proof of this formula can be found in 



6. Q-numbers 



The definition of g-numbers given in sec. 2 is not the only possible one. We have already seen the 
two-parameter deformed numbers, defined in eq. (2.11). Furthermore, a different definition of quantum 
numbers (the Q-numbers) has beenused in mathematics since the early ninetenth century, with rich 
literature existing on this subject eHe3. Q-numbers are defined as 



Q 3 



Q-l ' 



(6.1) 



where x can be a number or an operator and Q is a deformation parameter. Q is a real number 
(Q 7^ 0, 1). The notation Q = e T , where T a real number, will be often used. The subscript Q will be 
used in this review in order to distinguish deformed numbers defined as in eq. (6.1) from these defined 
by eq. (2.1). It is clear that in the limit Q — > 1 (or, equivalently, T — > 0) Q-numbers become ordinary 
numbers, i.e. [x]q — > x. 

A few examples of Q-numbers are given here: 



[0]c 



0. 



[1] 



Q 



[2] Q = Q + 1, [3} Q = Q 2 + Q + l. 



(6.2) 



Q-numbers clearly do not remain invariant under the substitution Q — > Q . One can easily prove 



that 



Q 



x-l 



Q-numbers are connected to g-numbers through the relation 

with Q = q 2 . 



(6.3) 



(6.4) 



(21 



The definitions of Q-factorials and Q-binomial coefficients still look like the ones given in eqs. (2.7)- 



\m\n\ 



[1]< 



[m-n] Q \[n] Q \ 

As it can be easily seen from eq. (6.1) under the substitution Q 

[n\ Q \ = Q^-W 

and 

Qk(n—k) 



Q one obtains 



n 



n 
k 



Q-factorials are connected to g-factorials by 



n 
k 



i/Q 



n 



g -n(n-l)/2 [n]g!) 



with 



Q = q 2 



(6.5) 
(6.6) 

(6.7) 
(6.8) 

(6.9) 



7. Q-deformed elementary functions 



The definitions of Q-deformed elementary functions look similar to these given in sec. 3. The 
Q-deformed exponential function is defined as 



00 a n 



n=o i n \Q- 

and satisfies the property 

eQ(x)ey Q (-x) = 1. (7.2) 

(Notice that eQ{x)eq{— x) 7^ 1.) 

The Q-deformed trigonometric functions are defined as 

00 T 2n+1 

sin Q (*) = £(-ir— - - 1 , (7.3) 

n=0 l Zn + iQ 

00 2ti 

cos Q (x) = £(-!)" T^n- (7-4) 

n=0 i 2n \Q l 



One can easily show that 



2 

Instead of the familiar identity sin 2 (x) + cos 2 (a;) = 1 one has 



sin Q (x) = — (e Q (ix) - e Q (-ix)) , (7.5) 
cosq(x) = i (eq(ix) + e Q (-ix)) . (7.6) 



sing(x) sin 1/Q (x) + cos Q (a;) cos 1/Q (x) = 1. (7.7) 

The above defined Q-deformed functions are examples of Q-deformed hypergeometric functions0l!>@. 

lynomif 
lials fa 



In addition Q-deformed polynomials, which are counterparts of the ordinary non deformed polynonijals, 
can be defined, such as Q-deformed Hermite polynomials and Q-deformed Laguerre polynomials 



8. Q-derivative 

Given the function f(x) one defines its Q-derivative Ell by the expression 

The similarity between this definition and the one of Q- numbers (eq. (6.1)) is clear. 
One can easily prove that 

s 1RT = [ " hx ^' (8 ' 2) 

which looks exactly like eq. (4.2). In addition one has 

D®e Q (ax) = ae Q (ax), (8.3) 

D®e 1/Q (ax) = ae 1/Q (aQx), (8.4) 

D^sing(ax) = acosg(ax), (8.5) 

D® cosq(o,x) = —asiiiQ(ax). (8.6) 

One can also easily see that sing (ax) and cosg(ax) are the linearly independent solutions of the Q- 
differential equation 

(D®) 2 u(x) + a 2 u{x) = 0, (8.7) 
while the functions sini/Q (aa;) and cosi/Q(ax) satisfy the equation 

{D^) 2 u(x) + a 2 u(Q 2 x) = 0. (8.8) 



The following Leibnitz rules can also be shown: 

D Q X (h(x)f 2 (x)) = ip^hix)) f 2 (Qx) + h{x) (D$f 2 (x] 



D Q X (h(x)f 2 (x)) = (D^\(x)) f 2 (x) + f x {Qx) (D^f 2 (x)) . 



One can further obtain 



dS 



h{x) _ (D2fi(xj)f 2 (x)-f 1 (x) (lg/ 2 ( 



h{x) f2(Qx)f 2 (x) 
For the second derivative of f(x) one has 

(D$ff(x) = (Q- iy 2 Q- l x- 2 {f(Q 2 x) - (Q + l)f(Qx) + Qf(x)} 

and by mathematical induction we obtain the general formula 



(D®) n f(x) = {Q- iyn Q -n{n-X)/2 x -n J- 



k=0 



n 
k 



{-ifQ^- 1 ^ 2 f{Q n ~ k x). 



(8.9) 
(8.10) 

(8.11) 
(8.12) 
(8.13) 



Q 



9. Q-integration 

In a wav analogous to that of sec. 5 the definite Q-integral of the function f(x) in the interval [0, 1] 
is defined B as follows 

(9.1) 



/ f(x)d Q x = (l-Q)Y / f(Q s )Q 



assuming that Q is real and \Q\ < 1, while for the definite integral of f(x) in the interval [0, oo], we 
have 

POO °° 

/ f(x)d Q x = (1 - Q) Y, Q s f(Q s )- (9-2) 
For the indefinite Q- integral of f(x) one has 

/oo 
f(x)d Q x = (1 - Q)x Q S f(Q Sx ) + constant. (9.3) 
s=-oo 

One can easily check that (^-differentiation and Q-integration are operations inverse to each other 

D<*J f(x)d Q x = f(x). (9.4) 
The formula for Q-integration by parts reads 

J (D^Mx)) f 2 (x)d Q x = h(x)f 2 (x) - J f x {Qx) (D^f 2 (x)) d Q x. (9.5) 



10. The g-deformed harmonic oscillator 

The interest for possible applications of quantum algebras in physics has been triggered in 1989 
by the introduction of the g-deformed harmonic oscillator 1100, of which earlier equivalent versions 
existed 0H. 

The g-deformed harmonic oscillator sME3EI'EI'e3 is defined in terms of the creation and annihilation 
operators a) and a and the number operator N, which satisfy the commutation relations 

[N,at] = cf, [N,a] = -a, (10.1) 

aa ] - g T Va = q ±N . (10.2) 

In addition the following conditions of hermitian conjugation hold (in the case of q being a real number 
or q being a root of unity) 

(a f ) f = a, iV f = N. (10.3) 

Eq. (10.1) is the same as in ordinary quantum mechanics, while eq. (10.2) is modified by the presence of 
the deformation parameter q. For g — ^ 1 it is clear that eq. (10.2) goes to the usual boson commutation 
relation [a,a^] = 1. An immediate consequence of (10.2) is that 

a ] a = [N], aa ] = [N + 1]. (10.4) 

Thus the number operator N is not equal to a^a, as in the ordinary case. The operators a) and a are 
referred to as q-deformed boson creation and annihilation operators respectively. 

The basis of the Fock space is defined by repeated action of the creation operator a) on the vacuum 
state, which is annihilated by a: 

(a r ) n 

a|0) = 0, |n>=^=|0>. (10.5) 



The action of the operators on the basis is given by 

N\n>=n\n>, (10.6) 



a)\n >= ^J[n + l]|n+ 1 >, (10.7) 
a\n >= J\n]\n - 1 > . (10.8) 



We remark that these equations look very similar to the ones of the ordinary case, the only difference 
being that g-numbers appear under the square roots instead of usual numbers. 
The Hamiltonian of the g-deformed harmonic oscillator is 

H = ^(aa f + a f a), (10.9) 

and its eigenvalues in the basis given above are 

E(n) = T ^-([n] + [n + l}). (10.10) 

One can easily see that for q real the energy eigenvalues increase more rapidly than the ordinary case, 
in which the spectrum is equidistant, i.e. the spectrum gets "expanded". In contrast, for q being a 
phase factor (g = e lT with r real) the eigenvalues of the energy increase less rapidly than the ordinary 
(equidistant) case, i.e. the spectrum is "compressed" . In particular, for q real (g = e r ) the eigenvalues 
can be written as 

fiuj sinh (t (n + |) ) 

E(n) = — M — , (10.11) 

v ; 2 sinhf ' v ; 



while for q being a phase factor (q = e lT ) one has 



. . ^sin (r (n + \ 
E(n) = v > - 2/7 . (10.12) 



2 sin | 



In both cases in the limit q — > 1 (r — > 0) the ordinary expression 

£7(n) = no; (n + (10.13) 

is recovered. 

In addition, the following commutation relation holds 

[a,a f ] = [N + 1] - [N]. (10.14) 
For q being a phase factor, this commutation relation takes the form 

(2jV+l)r 

[a,a^] = — 1 . (10.15) 

cos I 

It is useful to notice that the g-deformed boson operators a) and a can be expressed in terms of 
usual boson operators a+ and a (satisfying [a, a+] = 1 and N = a 1 a) through the relations 



The square root factors in the last equation have been called q-deforming functionals. 

For q being a primitive root of unity, i.e. q = e 2m l k (k = 2, 3, . . . ), it is clear the the representation 
of eqs (10.5)-(10.8) becomes finite-dimensional and has dimension k, since only the vectors |0 >, |1 >, 
. . . , \k — 1 > can be present. This case has been related to the system of two anyons @. In what follows 
we are going to assume that q is not a primitive root of unity. 

A discussion on the position and momentum operators of the g-deformed oscillator can be found in 

11. The Q-deformed harmonic oscillator 

A different version of the deformed harmonic oscillator can be obtained by defining the 
operators b, b + through the equations 

a = q^bq-^ 2 , a) = q^q-Wtf. (11.1) 

Eqs. (10.1) and (10.2) then give 

[N, 6+1 = 6+, [N,b] = -b, (11.2) 

bb ] - q 2 b ] b = 1. (11.3) 

This oscillator has been first introduced by Arik and Coon @ and later considered also by Kuryshkin 
One then easily finds that 

tfb=[N] Q , 66+ = [iV+l] Q , (11.4) 
where Q = q 2 and Q-numbers are defined in (6.1). The basis is defined by 



6|0) = 0, |n >= -^===|0 >, (11.5) 



n 



while the action of the operators on the basis is given by 

N\n>=n\n>, (11.6) 



&t|n>= y/[n + l] Q \n + l>, (11.7) 



b\n>=^[n] Q \n-l> . (11.8) 
The Hamiltonian of the corresponding deformed harmonic oscillator has the form 

H = — (66 f + 6+6), (11.9) 

the eigenvalues of which are 

E(n) = ^([n] Q + [n + l] Q ). (11.10) 

One can easily see that for Q = e T , where T > and real, the spectrum increases more rapidly than 
the ordinary (equidistant) spectrum, while for Q = e T , with T < and real, the spectrum is increasing 
less rapidly than the ordinary (equidistant) case. 

From the above relations, it is clear that the following commutation relation holds 

[b,tf) = Q N . (ii.il) 



12. The generalized deformed oscillator 

In addition to the oscillators described in the last two sections, many kinds of deformed oscillators 
have been introduced in the literature (see e3 for a list). All of them can be accommodated within the 
common mathematical framework of the generalized deformed oscillator e£E2I, which is defined as the 
algebra generated by the operators {1, a, a\ N} and the structure function $(x), satisfying the relations 

[a,N]=a, [a\N] = -a\ (12.1) 

a) a = $(iV) = [N], aa ] = <$>(N + 1) = [N + 1], (12.2) 

where <&(x) is a positive analytic function with $(0) = and iV is the number operator. ^From eq. 
(12.2) we conclude that 

AT = $- 1 (a t a), (12.3) 
and that the following commutation and anticommutation relations are obviously satisfied: 

[a,<f] = [N+l]-[N], {a,a f } = [N + 1} + [N]. (12.4) 

The structure function $(x) is characteristic to the deformation scheme. In Table 1 the structure 
functions corresponding to different deformed oscillators are given. They will be further discussed at 
the end of this section. 

It can be proved that the generalized deformed algebras possess a Fock space of eigenvectors 
|0 >, |1 >,..., \n >, ... of the number operator iV 

N\n >= n\n >, < n\m >= d~ nm , (12.5) 

if the vacuum state |0 > satisfies the following relation: 

a|0>=0. (12.6) 

These eigenvectors are generated by the formula: 

|n >= -^=fa t ) n |0 >, (12.7) 



where 



k=l k=l 



The generators a) and a are the creation and annihilation operators of this deformed oscillator algebra: 

,tl 



a\n >= \J[n\a\n — 1 >, a'\n >= J [n + l]a\n + 1 > . (12.9) 
These eigenvectors are also eigenvectors of the energy operator 



H = ^(aa f + a) a), (12.10) 



corresponding to the eigenvalues 



E{n) = (n) + $(n + 1)) = ^([n] + [n + 1]). (12.11) 

For 

$(n) = n (12.12) 
one obtains the results for the ordinary harmonic oscillator. For 

$( n ) = = [ n ] (12.13) 

one has the results for the g-deformed harmonic oscillator, while the choice 

$(n) = = [n] g (12.14) 

leads to the results of the Q-deformed harmonic oscillator. Many more cases are shown in Table 1, on 
which the following comments apply: 

i) Two-parameter deformed oscillators have been introduced , in analogy to the one-parameter 
deformed oscillators. 

ii) Parafermionic oscillators " of order p represent particles of which the maximum number which 
can occupy the same state is p. Parabosonic oscillators ^ can also be introduced. The deformed 
oscillator realization of the parafermionic and parabosonic oscillator has been studied in e3, while in 
an equivalent realization for the parabosonic oscillator has been introduced. 

iii) g-deformed versions of the parafermionic and parabosonic oscillators have also been introduced 



iv) g-deformed versions of the fermionic algebra @ have also been introduced as well 

as g-deformed versions of generalized g-deformed fermionic algebras H It has been proved, however, 
that g-deformed fermions are fully equivalent to the ordinary fermions §111. 

13. The physical content of deformed harmonic oscillators 

In order to get a feeling about the physical content of the various deformed harmonic oscillators it 
is instructive to construct potentials giving spectra similar to these of the oscillators. 

13.1. Classical potentials equivalent to the q-oscillator 

Let us consider the g-deformed harmonic oscillator first. For small values of r one can take Taylor 
expansions of the functions appearing there and thus find an expansion of the q- number [n] of eq. (2.1) 
in terms of powers of r 2 . The final result is 

r 2 t 4 t 6 

\n]=n± —in - n 3 ) + (7n - 10n 3 + 3n 5 ) ± (31n - 49n 3 + 21n 5 - 3n 7 ) + . . . , (13.1) 

[J 6 v ) t 36Q v 15120 v ' ' v ' 

where the upper (lower) sign corresponds to q being a phase (real). Using this expansion the energy of 
the g-deformed harmonic oscillator of eq. (10.10) can be rewritten as 

E(n)/(huj) = (n + l)(l±^) T ^ („ + i) 3 + . . . (13.2) 



Table 1: Structure functions of special deformation schemes 







Reference 


i 


X 


harmonic oscillator, bosonic 
algebra 


ii 


q x —q~ x 
q-q- 1 


g-deformed harmonic oscil- 
lator SB 


iii 


f-1 
q-1 


Arik— Coon, Kuryshkin or 






Q-deformed oscillator @H 


iv 


q-p- 1 


2-parameter deformed oscil- 
lator iii 


V 


x(p + 1 — x) 


parafermionic oscillator S 


vi 


sinh(rx) sinh(r(p+l— x)) 
sinh 2 (r) 


g-deformed parafermionic 
oscillator HH 


vii 


xcos {nx/z) + [x +p — 1) sin (irx/z) 


parabosomc oscillator c-a 


viii 


smh(-nr) cosh(r(a;+2JVo — 1)) qq^2/—^, /o\ i 
sinh(r) cosh(r) V / / 
sinh(T(!r+2Aro-l)) cosh(ra:) ^2/^ /Q\ 
sinh(T) cosh(r) V / / 


g-deformed parabosonic os- 
cillator EBEd 


ix 


sin 2 nx/2 


fermionic algebra S 


X 


y Sill 71 J, / Z 


g-ueioimeu ieimiomc aige- 
bra SfflfflS 


xi 


!-(-<?)* 
1+9 


generalized g-deformed 
fermionic algebra 


xii 


X n 




xiii 







On the other hand, one can consider the potential 

V(x) = V + kx 2 + Xx 4 + fix 6 + £x 8 + . . . . (13.3) 

If A, fi, £ are much smaller than k, one can consider this potential as a harmonic oscillator potential 
plus some perturbations and calculate the corresponding spectrum through the use of perturbation 
theory S (see also subsec. 34.2). In order to keep the subsequent formulae simple, we measure x in 
units of {h/ (2muj)) l l 2 . Using standard first order perturbation theory one finds that the corresponding 
spectrum up to the order considered is 

E(n) = E + {2k + 25/i)(n + \) + (6A + 245£)(n + \f + 20/x(n + i) 3 + 70£(n + -) 4 (13.4) 

By equating coefficients of the same powers of (n + \) in eqs. (13.2) and (13.4), we can determine the 
coefficients appearing in the expansion of the potential given in eq. (13.3). The final result for the 
potential, up to the order considered, is 

^) = (^±y)* 2 T^ 6 . (13.5) 

We see therefore that to lowest order one can think of the g-oscillator with a small value of the parameter 
r as a harmonic oscillator perturbed by a x 6 term. It is clear that if we go to higher order, the next 
term to appear will be proportional to r 4 x 10 . 

The results of this subsection are corroborated by an independent study of the relation between the 
g-deformed harmonic oscillator and the ordinary anharmonic oscillator with x 6 anharmonicities 



A different interpretation of the g-deformed oscillator as the nonlinear oscillator with the frequency 
dependent on the amplitude has been given in 0. 



13.2. Classical potentials equivalent to the Q-deformed oscillator 

Similar considerations can be made for the Q-oscillator. Defining Q = e T it is instructive to construct 
the expansion of the Q- number of eq. (6.1) in powers of T. Assuming that T is small and taking Taylor 
expansions in eq. (6.1) one finally has 



Nq = n + ~(n 2 -n) + y^(2™ 3 - 3n 2 + 1) + — (n 4 - 2n 3 + n 2 ) + 



;i3.6) 



Then the corresponding expansion of the energy levels of the oscillator of eq. (11.10) is 
E{n)/{huj) = E' + (1 - 



T T 2 T 3 w 1 T T 

1 h . . .)(n + -) + ( 

2 8 16 A 2 ; v 2 4 



. . .)(n+ -Y 

48 A 2 ; 



+ ( 



6 



h . . .)(n + -) 3 + ( — 

12 A 2 J V 24 



■){n+-f. 



(13.7) 



Comparing this expansion to eq. (13.4) and equating equal powers of (n + 1) we arrive at the following 
expression for the potential 



T 2 

V(x) = — 
v ' 12 



4 T T 2 T 3 2 T T 2 T 3 4 T 2 
v 2 ~ I ~ 24 + 48 ^ + ^12 ~ 24 ~ 144 ^ + ^120 



T 3 T 3 

)X + : 3T 



13.8) 



240' 1680 

Thus to lowest order one can think of the Q-oscillator with small values of the parameter T as a 
harmonic oscillator perturbed by a x 4 term. A similar expression is found, for example, by Taylor 
expanding the modified Poschl-Teller potential Ell, which, among other applications, has been recently 
used in the description of hypernuclei The modified Poschl-Teller potential (see also subsec. 

34.2) has the form 

V(x)pt = 7^-^. (13.9) 



Its Taylor expansion is 



V(x) PT = A(-l + aV 



cosh(ax) 2 



-a'.r ' + —a"/' - 
3 45 



;i3.io) 



We remark that this expansion contains the same powers of x as the expansion (13.8). Furthermore, 
the signs of the coefficients of the same powers of x in the two expansions are the same for T < 0. 



13.3. WKB-EPs for the q-deformed oscillator 

The potentials obtained above are only rough lowest order estimates. More accurate methods exist 
for constructing WKB equivalent potentials (WKB-EPs) giving (within the limits of WKB approxima- 
tion) the same spectrum as the above mentioned oscillators. A method bywhich this can be achieved 
has been given by Wheeler @ and is described by Chadan and Chabatier 0. Applyingthis method to 
the g-deformed oscillator (sec. 10) with q being a phase factor one finds the potential E§E1 



V{x) 
where 



> 2sin(r/2) < 



-x 



1 - 



/ X 



15 \2K 



4448 



1575 \2R 



x 



R, = - — 



h 2 V /2 /2sin(r/2) > 



t \2m J \ huo 
while for the g-deformed oscillator with q real one has 



V{x) 



T 



2sinh(r/2) 



-x 



1 + TT 



/ X 



15 \2R h 



4 4448 
+ 



1575 \2R 



x 



345344 / x 



675675 V2iL 



12 



+ 



1/2 



345344 / x \ 12 
+ 675675 \2R~J + 



13.11) 



'13.12) 



13.13) 



where 



The results of this subsection are corroborated by an independent study of the relation between the 
g-deformed harmonic oscillator and the ordinary anharmonic oscillator with x 6 anharmonicities!!. 

13. 4- WKB-EPs for the Q-deformed oscillator 

Using the same technique one finds that the WKB equivalent potential for the Q-deformed harmonic 
oscillator (sec. 11) takes the form 



1 2 / x \ 2 23 / x \ 4 134 / x \ 6 . r )297 / ./■ \ s 



3 \R'J 45 \R'J 315 \R'J 14172 \R' 



'13.15) 



where 



and 



^«- 2V g (V g + l)I 

\ cum(Q + l) )^/2 y ' 

We remark that this WKB-EP contains all even powers of x, in contrast to the WKB-EPs for the 
g-oscillator (eqs (13.11), (13.13)), which contain only the powers x 2 , x 6 , x 10 , .... This is in agreement 
with the lowest order results obtained in subsecs 13.1 and 13.2. 

In this section approximate expressions for some potentials having similar energy spectrum with some 
deformed oscillators have been given. Another point of view is the search for phase space coordinates 
having exact symmetries, like the operators involved in g-deformed algebras in the case of real q 0. In 
the case of q being a root of unity the corresponding phase space is discrete 

14. The quantum algebra su g (2) 

Quantum algebras are generalized versions of the usual Lie algebras, to which they reduce when the 
deformation parameter q is set equal to unity. A simple example of a quantum algebra is provided by 
su„ (2) US, which is generated by the operators J + , Jq, J_, satisfying the commutation relations 

[Jo, J±] = ± J±, (14.1) 

[J + ,J_] = [2J ], (14.2) 

with Jq = J , (J+Y = J-. We remark that eq. (14.1) is the same as in the case of the ordinary su(2) 
algebra, while eq. (14.2) is different, since in the usual su(2) case it reads 

[J+, J-) = 2J . (14.3) 

In the rhs of eq. (14.3) one has the first power of the Jo operator, while in the rhs of eq. (14.2) one has 
the g-operator [2 Jo], defined in sec. 2. Because of eqs (2.2) and (2.3) it is clear that if one writes the 
rhs of eq. (14.2) in expanded form all odd powers of Jq will appear: 



'2rJ (2rJ ) 3 (2rJ, 



^• J -^^H7 ) [^r + ^ + ^r L + ^) fM q ^ eT ' {UA) 

„ ,. 1 (2tJ (2rJ ) 3 (2tJ„) 5 \ , T 



Thus su g (2) can be loosely described as a nonlinear generalization of the usual su(2): While in usual 
Lie algebras the commutator of two generators is always producing a linear combination of generators, 
in the case of quantum algebras the commutator of two generators can contain higher powers of the 
generators as well. 

The irreducible representations (irreps) D J of su g (2) (which have dimensionality 2 J + 1) are deter- 
mined by highest weight states with J = 0, |, 1, . . .. The basic states \J,M> (with — J < M < J) are 
connected with highest weight states \J, J > as follows 



\J,M>-- 



[J + M}\ 
\ [2J)\[J-M)\ 



{J- 



, J-M\ 



J,J>, 



(14.6) 



with J + \ J, J >= and < J, J| J, J >= 1. The action of the generators of the algebra on these basic 
vectors is given by 

J \J,M >= M\J,M >, (14.7) 



J±\J,M >= J[Jt M}[J±M + 1]|J,M±1 > . 



:i4.8) 



These expressions look similar to the ones of the usual su(2) algebra, the only difference being that 
g-numbers appear under the square root instead of ordinary numbers. 

The second order Casimir operator of su g (2) is determined from the condition that it should commute 
with all of the generators of the algebra. The resulting operator is 



a = + [j ] [j + 1] = j+j- + [j ] [j - 1]. 

Its eigenvalues in the above mentioned basis are given by 

C%\J,M>= [J][J + 1]\J,M>, 



;i4.9) 



(14.10) 



while for the usual su(2) the eigenvalues of the Casimir operator are J( J+ 1). One can easily check that 
for real q (q = e T with r real) the eigenvalues [J][J+1] produce a spectrum increasing more rapidly than 
J(J+ 1) (an expanded spectrum), while for q being a phase factor (q = e lT with r real) the eigenvalues 
[J] [J + 1] correspond to a spectrum increasing less rapidly than J(J + 1) (a compressed spectrum). 

It should be noticed that the generators J+, Jo, J_ of su g (2) are connected to the generators j + , jo, 
2- of the usual su(2), which satisfy the commutation relations 



[jo , j±] = ±j± , []+ ,J-] = 2 jo , 



through the relations iH 



Jo = Jo, 



J4 



bo + j][jo 



J_ 



bo + j] [m - 1 - j] 
\ bo + j]bo-i-j] 



\ Uo • ./)i./o- 1-./V 
where j is determined by the relation for the second order Casimir operator of su(2) 

C = 3-3+ + Jo(jo + 1) = j+j- + J0U0 - 1) = j{3 + !)• 



;i4.ii) 



14.12) 



(14.13) 



15. Realization of su q (2) in terms of g-deformed bosons 

Realizations of Lie algebras in terms of (ordinary) bosons are useful not only as a convenient math- 
ematical tool, but also because of their applications in physics Ea. In the case of quantum algebras 
it turns out that boson realizations are possible in terms of the g-deformed boson operators already 
introduced in sec. 10. 



In the case of su g (2) the generators can be mapped onto g-deformed bosons as follows 
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J + = a\a 2 , J_ = ajai, J = ^(JVi - N 2 ), (15.1) 



where a\, ai and iVj are g-deformed boson creation, annihilation and number operators as these intro- 
duced in sec. 10. One can easily prove that the boson images satisfy the commutation relations (14.1) 
and (14.2). For example, one has 



[J+, J_] = a\a 2 a\ ai - a\a x a\a 2 = [iVx][iV 2 + 1] - [JV X + 1][N 2 ] = [N t - N 2 ] = [2J ] 

where use of the identity (2.6) has been made. 

In the g-boson picture the normalized highest weight vector is 



(15.2) 



\JJ >-- 



[2J]! 



0>, 



;i5.3) 



while the general vector \JM > is given by 

\JM >= 



[a\) J+M 



10 > 



[J + M]\ J[J-M] 



(15.4) 



It should be noticed that it was the search for a boson realization of the su g (2) algebra that led to 
the introduction of the g-deformed harmonic oscillator in 1989 EH. 

Starting from su g (2) one can formulate a g-deformed version of angular momentum theory. Some 
references are listed here: 

i) Clebsch-Gordan coefficients for su g (2) can be found in@ffl@S§§il@. 

ii) 6-j symbols for su g (2) can be found in [H^'EH. 
hi) 9-j symbols for su g (2) can be found in E. 
iv) 3n-j symbols (with n=l, 2, 3, 4, 5) for su g (2) can be found in 



IOC 



101 



102 



L03 



v) The g-deformed version of the Wigner-Eckart theorem can be found in 

vi) Irreducible tensor operators for su g (2) can be found in E2I. 

In addition, it should be noticed that a two-parameter deformation of su(2), labelled as su Pi? (2) 

has beeri intrnrhi ceH E3Elll£Hl2l'E2l'E£3'E2l . Clebsch-Gordan coefficients for su Pi g(2) have been discussed 
j n |i 09U1 1 o| Ji i ilJ 1 1 2I. Fit^ 
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The way in which the algebra su q (2) can be realized in terms of the Q-deformed bosons of sec. 11 
is given in subsec. 35.1 (see eqs (35.13)-(35.16)). 



16. The quantum algebra su g (l,l) 

In this section we shall give a brief account of the algebra su q (l,l) ElBllll. 
In the classical case the so(2,l) generators satisfy the commutation relations 



[K u K 2 ] = -iK 3 , [K 2 , K 3 ] = %K U [K 3 , K x 



iK. 



2- 



(16.1) 



which differ from the classical so(3) commutation relations in the sign of the r.h.s. of the first commu- 
tator. Defining 

K + = K x + iK 2 , K„ = K x -iK 2 , K 3 = K z , (16.2) 
one obtains the su(l,l) commutation relations 



[K z ,K d 



±K 4 



[K + ,K„ 



-IK, 



(16.3) 



which differ from the familiar su(2) commutation relations in the sign of the r.h.s. of the last commu- 
tator. The generators of su(l,l) accept the following boson representation til 



a\a 2 , 



a±a 2 , 



(a^ai + a 2 a 2 + 1), 



(16.4) 



where a\, a\, a 2 , a 2 satisfy usual boson commutation relations. 
The second order Casimir operator of so(2,l) is llll 



C 2 [so(2,l)] = -(K 2 1 +K 2 2 -K 2 3 ), 



(16.5) 



while for su(l,l) one has 

C 2 [su(l, 1)] = [K ][K - 1] - K + K_ = [K Q ][K + 1] - K„K + . 
In the quantum case, the generators of su q (l,l) satisfy the commutation relations K u 1 ! 1 1 ' 

[Kq, K ± ] = ±K ± , [K+, K-] = -[2K \. 
The generators of su 9 (l,l) accept the following boson representation 



K 4 



K_= ai a 2 , K = -{Ni + N 2 + l) 



(16.6) 



;i6.7) 



16.8) 



where the bosons a\, di (i = 1, 2) satisfy the usual g-boson commutation relations. 
The second order Casimir operator of su g (l,l) is 



Its eigenvalues are 
where 



C 2 [su,(l, 1)] = [K ][K - 1] - K + K_ = [K }[K + 1] - K_K, 
nil 

C 2 [sug(l, >= [k][k - 1]\k/i >, 
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K 



1 + |m — n 2 | 



1 + rii + n 2 



since the basis has the form >= \rii > \n 2 >, with 



7li > 



«l)" l |o > . 



(16.9) 
16.10) 
16.11) 

16.12) 



71; 



In this basis the possible values of /i are given by 

li = k, k + 1, k + 2, . . . , 



16.13) 



up to infinity, while k may be any positive real number. The action of the generators on this basis is 
given by 

K \nfi) = (16.14) 



K ± \kh) = ^±K][/iTK±l]|¥±l>' (16.15) 

Clebsch-Gordan and Racah coefficients for su g (l,l) can be found in EiBBBlOJIIi Fiir lhernnje, a 
two-parameter deformed version of su g (l,l), labelled as su Pi9 (l,l), has been introduced & 



106 



12C ,12 



rmre i 
12 jp.22 



17. Generalized deformed su(2) and su(l,l) algebras 

In the same way that in addition to the g-deformed oscillators one can have generalized deformed 

oficial case 



oscillators, it turns out that generalized deformed su(2) algebras, containing su g (2) as 



and having representation theory similar to that of the usual su(2), can be constructed l 123 H 124 l. It has 
been proved that it is possible to construct an algebra 



±J4 



[j+, J-] = $(J (Jo + 1)) - HMJo - 1)), 



(17.1) 



where Jq, J + , J_ are the generators of the algebra and $(x) is any increasing entire function defined 
for x > —1/4. Since this algebra is characterized by the function $, we use for it the symbol su$(2). 
The appropriate basis \l,m > has the properties 



J \L,M>= M\L,M>, 



(17.2) 



J+\L, M >= yJ$(L(L + 1)) - $(M(M + 1))|L, M + 1 >, 



(17.3) 



where 



and 



J-\L,M >= J$(L(L + 1))-$(M(M-1))\L,M-1 >, 



1 3 5 
2' 2' 2' ^' 



M 



-L, -L + 1. -L + 2, . . . , L - 2, L - 1. L. 



The Casimir operator is 



C = J_ J + + $(J (J + 1)) = J+J- + $(Jo(Jo - 1)), 
its eigenvalues indicated by 

C\L,M >= + 1))\L,M > . 

The usual su(2) algebra is recovered for 



$(x(x + 1)) = x(x + 1), 



while the quantum algebra su 9 (2) 



[Jo, J± 



±J4 



[J + ,J_] = [2J ( 



Ok, 



(17.4) 
(17.5) 
(17.6) 

(17.7) 
(17.8) 

(17.9) 
(17.10) 



occurs for 

&{x(x + l)) = [x][x + l]. (17.11) 

The su$(2) algebra occurs in several cases, in which the rhs of the last equation in (17.1) is an odd 
function of Jo It can be seen that other algebraic structures, like the quadratic Hahn algebra QH(3) 
^3 and the finite W algebra Wo tH3 can be brought into the su$(2) form, the advantage being that the 
representation theory of su$(2) is already known. It can also be proved that several physical systems, 
like the isotropic oscillator in a 2-dim curved sp ace wi th constant curvature the Kepler system 

in a 2-dim curved space with constant curvature t23'E^I, and the system of two identical particles in two 
dimensions t23 can also be put into an su$(2) form. More details can be found in H 25 
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A unified description of several deformed su(2), su(l,l) and osp(l,2) algebras has been given in 
using Verma modules techniques. Nonlinear deformed su(2) algebras involving two deforming functions 
have been constructed in and have been endowed with a double Hopf algebraic structure, referred 
to as a two-colour quasitriangular Hopf algebra, in 



Nonlinear si (2) algebras subtending generalized angular momentum theories have been introduced 
in and their representation theory has been developed. In particular the quadratic sl(2) and cubic 
sl(2) algebras have been studied in detail Bf. 

A general method for developing nonlinear deformations of the su(2) and su(l,l) algebras has been 
developed in E1M'EHIE13. Using this method nonlinear deformations of the su(l,l) algebra related to the 
Poschl-Teller potential have been obtained in Infill. 

18. Generalized deformed parafermionic oscillators 

It turns out that the generalized deformed su$(2) algebras mentioned in the last section are related 
to generalized deformed parafermionic oscillators, which we will therefore describe here. 

It has been proved E1PJ that any generalized deformed parafermionic algebra of order p can be written 
as a generalized oscillator (sec. 12) with structure function 



F(x) = x(p + 1 — x) (A + fj,x + vx 2 + px 3 + ax 4 + 
where A, /x, u, p, a, . . . are real constants satisfying the conditions 

A + px + vx 2 + px 3 + ax 4 + . . . > 0, x G {1, 2, . 



,P}- 



:i8.D 



[18.2) 



Considering an su$(2) algebra with structure function 



$(J (Jo + 1)) = AJ (J Q + 1) + B(J (J + l)) 2 + C(J (Jo + l)) 3 , 



(18.3) 



and making the correspondence 



J + -> A\ J_ -> A, J -> iV, (18.4) 

one finds that the su$(2) algebra is equivalent to a generalized deformed parafermionic oscillator of the 
form 

F(N) = N(p + 1 - N) 

h(p 2 (p + 1)C + pB) + (p 3 C + (p - 1)B)N + ((p 2 - p + 1)C + B)N 2 + (p - 2)CN 3 + CW 4 ], (18.5) 
if the condition 

A + p(p + l)5 + p 2 (p + 1) 2 C = (18.6) 

holds. The condition of eq. (18.2) is always satisfied for B > and C > 0. 

In the special case of C = one finds that the su$(2) algebra with structure function 

$( J ( J + 1)) = AJ ( Jo + 1) + fl( Jq( Jo + l)) 2 (18.7) 

is equivalent to a generalized deformed parafermionic oscillator characterized by 

F(N) = BN(p + 1 - N)(-p + (p - 1)JV + N 2 ), (18.8) 

if the condition 

A + p(p+ 1)5 = (18.9) 

is satisfied. The condition of eq. (18.2) is satisfied for B > 0. 

Including higher powers of Jo(Jo + 1) in eq. (18.3) results in higher powers of N in eq. (18.5) and 
higher powers of p(p+l) in eq. (18.6). If, however, one sets B = in eq. (18.7), then eq. (18.8) vanishes, 
indicating that no parafermionic oscillator equivalent to the usual su(2) rotator can be constructed. 



It turns out that several other mathematical structures, like the finite W algebras Wq £23 and W3 



(see subsec. 36.5) can be put into the generalized deformed parafermionic oscillator form. The same 
is true for several physical systems, su ch as t he isotropic oscillator and the Kepler problem in a 2-dim 
curved space with constant curvature l 128 H 129 l, and the Fokas-Lagerstrom tli, Smorodinsky-Winternitz 
Ell, and Holt Ell potentials. Further details can be found in llM 

A detailed discussion of the representation theory of several deformed oscillator algebras can be 
found in ES. 



19. The su g (2) rotator model 

It has been suggested by Raychev, Roussev and Smirnov 111! and independently by Iwao 0@ that 
rotational spectra of deformed nuclei can be described by the g-deformed rotator, which corresponds 
to the 2nd order Casimir operator of the quantum algebra su 9 (2), already studied in sec. 14. We shall 
show here that this assumption works and discuss the reasons behind this success, as well as the relation 



ila between the su g (2) model and the Variable Moment of Inertia (VMI) model (see sec. 20). 
The g-deformed rotator corresponds to the Hamiltonian 

H = l } C 2 {m q (2)) + E , (19.1) 

where I is the moment of inertia and Eq is the bandhead energy (for ground state bands Eq = 0). For 
q real, i.e. with q = e T with r real, the energy levels of the g-rotator are 

E(J) - + 1] + E = + 1)} + (W.2) 

For q being a phase, i.e. q = e lT with r real, one obtains 

EW = kj]U + 1] + E = l 3i ^ J > *W + £„. (19.3) 
21 21 sin it) 



Raychev et al. £13 have found that good fits of rotational spectra of even-even rare earths and 
actinides are obtained with eq. (19.3). It is easy to check that eq. (19.2) fails in describing such 
spectra. In order to understand this difference, it is useful to make Taylor expansions of the quantities 
in the numerator of eq. (19.2) (eq. (19.3)) and collect together the terms containing the same powers 
of J(J + 1) (all other terms cancel out), finally summing up the coefficients of each power. In the first 



case the final result is 



E{J) = E + - ( ^-j 1/a ( T) J(j + 1) + r^- h/2{ r)(J(J + 1)) 



2r 2 nr 

+ 



3 V 2r 



h/2(r)(J(J + l)) 3 + y yj/ 7/2 (r)( J( J + l)) 4 + • • • (19.4) 
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where J^pl n+ i(r) are the modified spherical Bessel functions of the first kind 
In the second case (eq. (19.3)) following the same procedure one obtains 

E{J) = E + ^—L^(j ( T )J(J + 1) - rj x {T){J{J + l)) 2 
+It*j 2 (t)(J(J + l)) 3 - \t*j,{t){J{J + l)) 4 + ^t%(t)(j(j + l)) 5 -...), (19.5) 
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where j n ij) are the spherical Bessel functions of the first kind 
Both results are of the form 

E(J) = AJ(J + 1) + B(J{J + l)) 2 + C(J(J + l)) 3 + D(J{J + l)) 4 + . . . . (19.6) 

Empirically it is known that nuclear rotational spectra do show such a behaviour, the coefficients A, 
B, C, D, ... having alternating signs (starting with A positive) and magnitudes dropping by about 3 
orders of magnitude each time one moves to the next higher power of J(J + 1) El^lill. 

It is interesting to check if the empirical characteristics of the coefficients A, B, C, D are present 
in the case of the expansions of eqs. (19.2), (19.3), especially for small values of r. (Since we deal 
with rotational spectra, which are in first order approximation described by the usual algebra su(2), we 
expect r to be relatively small, i.e. the deviation of su g (2) from su(2) to be small. This is in agreement 
to the findings of ll 4 !, where r is found to be around 0.03.) 

One can easily see that in eq. (19.2) it is impossible to get alternating signs, while in eq. (19.3) the 
condition of alternating signs is readily fulfilled. This fact has as a result that the energy levels given by 
eq. (19.2) increase more rapidly than the levels given by the J(J + 1) rule, while the levels given by eq. 
(19.3) increase less rapidly than J(J +1). In order to check the order of magnitude of the coefficients 
for small values of r, it is useful to expand the spherical Bessel functions appearing in eq. (19.3) and 
keep only the lowest order term in each expansion. The result is 

E(J) = E + 1(J(J + 1) - j(J(J + l)) 2 + ^(J(J + l)) 3 

-—(J(J + l)) 4 + -^^(J(J + I)) 5 -•••)• (19.7) 
315 V V 14175 V V ;; ; V ; 

We remark that each term contains a factor r 2 more than the previous one. For r in the area of 0.03, 
t 2 is of the order of 10~ 3 , as it should. We conclude therefore that eq. (19.6) is suitable for fitting 
rotational spectra, since its coefficients have the same characteristics as the empirical coefficients of eq. 
(19.6). 

Extended comparisons of the su g (2) pre diction s to e xperimental data for ground state bands of rare 
earth and actinide nuclei can be found in I 146 h 149 h 153 h 154 I . More recently, the su g (2) formalism has been 



used for the description of (3- and 7-bands of deformed rare earths and actinides, with satisfactory 



results 



It is necessary for E(J) to be an increasing function of J. In order to guarantee this in eq. (19.3) 
one must have 

r(J+l)<|. (19.8) 



In the case of r = 0.036 (as in 232 U in Eia), one finds J < 42, this limiting value being larger than the 

Similarly, for r = 0.046 (as in 178 Hf 



highest observed J in ground state bands in the actinide region 
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one finds J < 32, this limiting value being again higher than the highest observed J in ground 



state bands in the rare earth region 
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It should be mentioned that a method of fixing the deformation parameter by demanding that the 



relevant Hilbert space be large enough to contain all the irreps of the system has been suggested 
the framework of Barnett-Pegg theory, with good results. 
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20. Comparison of the su q (2) model to other models 



20.1. The Variable Moment of Inertia (VMI) model 

In lowest order approximation rotational nuclear spectra can be described by the formula 

EW - (20.1) 

where B is the moment of inertia of the nucleus, which is assumed to be constant. However, in order to 
get closer agreement to experimental data, one finds that he has to include higher order terms in this 
formula, as shown in eq. (19.6). 

Another way to improve agreement with experiment is to let the moment of inertia to vary as a 
function of the angular momentum J. One thus obtains the Variable Moment of Inertia (VMI) model 
In this model the levels of the ground state band are given by 
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E(J) = + \C{Q{J) - e ) 2 , (20.2) 

where Q(J) is the moment of inertia of the nucleus at the state with angular momentum J, while C 
and 0o are the two free parameters of the model, fitted to the data. The parameter Oo corresponds to 
the ground state moment of inertia, while instead of the parameter C it has been found meaningful to 
use the parameter combination 

ff = (20 ' 3) 

which is related to the softness of the nucleus. The moment of inertia at given J is determined through 
the variational condition 

dE(J) 

de(jy 

which is equivalent to the cubic equation 



0, (20.4) 



6(J) 3 - 6(J) 2 0o - ^±12 = o. (20.5) 



This equation has only one real root, which can be written as 

0(J) 



j(j + i) eg / (j(j + i)) 2 ep{J + K 

\ AC 27 V 16C 2 54C 



+ "i^C— + 27~i 16C 2 + 54C + T- (m6) 
Expanding the roots in this expression one obtains 

0(J) = e (i + a J(J + 1) - 2a 2 (J(J + l)) 2 + 7a 3 ( J( J + l)) 3 - 30a 4 ( J( J + l)) 4 + ...). (20.7) 



Using eq. (20.7) in eq. (20.2) one obtains the following expansion for the energy 



E(J) = ±-{J{J + 1) - \a{J{J + l)) 2 + a 2 ( J( J + l)) 3 - 3a 3 ( J( J + l)) 4 + ...). (20.* 



Empirically it is known that for rotational nuclei the softness parameter o is of the order of 10 -3 . 
Therefore the expansion of eq. (20.8) has the same characteristics as the expansion of eq. (19.6) 
(alternating signs, successive coefficients falling by about 3 orders of magnitude). 

20.2. Comparison of the su q (2) model to the VMI and related models 

We now turn to the comparison of the expansion of eq. (19.5) to the Variable Moment of Inertia 
(VMI) model, discussed in the previous subsection. Comparing eqs (19.5) and (20.8) we see that both 
expansions have the same form. The moment of inertia parameter I of (19.5) corresponds to the ground 
state moment of inertia 0o of (20.8). The small parameter of the expansion is r 2 in the first case, while 
it is the softness parameter 1/(2COq) in the second. However, the numerical coefficients in front of 
each power of J( J + 1) are not the same. 

In Eil a comparison is made between the parameters obtained by fitting the same spectra by the 
su g (2) and VMI formulae. The agreement between 1/(2 7) a nd l/(20 o ) is very good, as it is the 
agreement between r 2 and a as well. Therefore the known li^l smooth variation of Go and a with the 
ratio i? 4 = E(4)/E(2) is expected to hold for the parameters I and r 2 as well. This is indeed seen in 



ii<: 



The difference between the expansions of eqs (19.5) (or (19.7)) and (20.8) is also demonstrated by 
forming the dimensionless ratios AC/(AB 2 ) and A 2 D/(2AB 3 ) in eq. (19.6). In the case of eq. (20. 



both quantities are equal to 1, as expected, since it is known that the VMI is equivalent to the 

Harris expansion in which both quantities are known to be equal to 1. In the case of eq. (19.7) 
the corresponding values are AC /(AB 2 ) = 1/10 and A 2 D/(2AB 3 ) = 1/280. According to the empirical 



values of these ratios given in l 152 H 162 l, the ratios given from eq. (19.7) are better than the ratios given 



by eq. (20.8), especially the second one. 
20.3. The hybrid model 



The hybrid model of nuclear collective motion l 163 H 164 H 165 H 166 l has been introduced in order to provide 
a link between the two successful w ays of describing low-lying nuclear excitations: th e exten ded form of 



the Bohr-Mottelson model (BMM) and the Interacting Boson Model (IBM) EafiatZa (see sees 27, 

31 for more details). The hybrid model combines the advantages of both models, i.e. the geometrical 
significance of the collective coordinates inherent in the extended BMM, and the use of group theoretical 
concepts characterizing IBM. In the framework of Hie rotational limit of the hybrid model, associated 
with the chain u(6)Dsu(3), Partensky and Quesne E3EZ3, starting from the fact that in the geometrical 
description the square of the deformation is proportional to the moment of inertia of the ground state 
band, proved that the energy levels of the ground state band are given by 

Eif > = WTl)TB JiJ + 1) ' (20 ' 9) 

with A being a free parameter and B given by 

B = 8N 2 + 22N - 15, (20.10) 

where N is the sum of the number of valence proton pairs (or proton-hole pairs, when more than half 
of the proton valence shell is filled) N n and the number of the valence neutron pairs (or neutron-hole 
pairs, when more than half of the neutron valence shell is filled) N v . We remark that in the framework 
of this model the moment of inertia 

Q(J,N)= J(J \ l } + B (20.11) 
depends on both the angular momentum J and the valence pair number N. 



It is instructive to expand E(J) in powers of J(J+1) 



E{J)=A±{-lf 



fe+i 



k=0 



+ 1) - + I)) 2 + ^(J(J + I)) 3 - + I)) 4 + •••)• (20.12) 



Comparing the present expansion to the one of eq. (19.7) for the su g (2) model one has 




87V 2 + 22/V- 15 : 



(20.13) 



obtaining thus a connection between the r parameter of the su g (2) model and a microscopic quantity, 
the valence nucleon pair number. ^From this equation it is clear that r decreases for increasing N. Thus 
the minimum values of r are expected near the midshell regions, where the best rotators are known to 
be located. This is a reasonable result, since r describes the deviations from the su(2) (rigid rotator) 
limit. It is clear that the minimum deviations should occur in the case of the best rigid rotators. 

From eq. (20.13) one can obtain for each nucleus the value of rfrom the number of valence nucleon 
pairs present in it. These predictions for r have been compared in tla to the values of r found empirically 
by fitting the corresponding spectra, with good results in both the rare earth and the actinide regions. 
In addition eq. (20.13) indicates that in a given shell nuclei characterized by the same valence nucleon 
pair number N will correspond to the same value of r. Such multiplets have been studied in the 
framework of the hybrid model by liZ3. 

It is worth noticing that fits of 71-bands in the rare earth region (Er, Yb isotopes) EH give r 
parameters very similar to the ones coming from fitting the corresponding ground state bands, in 
addition exhibiting the same as the one mentioned above behaviour of r as a function of N. 

Taking further advantage of the above connection between the su g (2) model and the hybrid model, 
the parameter r has been connected to the nuclear deformation parameter /3, as well as to the 
electromagnetic transition probabilities B(E2:2+ — > 0^). Since both /3 and B(E2:2+ — > 0^) are known to 
increase with increasing collectivity, r is expected to decrease with increasing (3 or increasing B(E2:2 1 f 
Of). This expectation is corroborated by the results reported in 
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20.4- Other models 



It should be noticed that an empirical formula very similar to that of eq. (19.3) has been recently 
proposed by Amal'skii EZ3 on completely different physical grounds. The formula reads 



(20.14) 



where A and B are free parameters. This formula should be compared to eq. (19.3), with which it is 
almost identical. 

A different formula, also giving very good fits of rotational spectra, has been introduced by Celeghini, 
Giachetti, Sorace and Tarlini E^l, based on the g-Poincare rotator. 



21. Electromagnetic transitions in the su g (2) model 

We have already seen that the su 9 (2) formalism provides an alternative to the VMI model, the 
deformation parameter q being connected to the softness parameter of the VMI model. 

The stretching effect present in rotational energy levels, which can equally well be described in 
terms of the VMI model and the su g (2) symmetry, should also manifest itself in the B(E2) transition 
probabilities among these levels. If deviations from the su(2) symmetry are observed in the energy levels 
of a band, relevant deviations should also appear in the B(E2) transitions connecting them. In the case 
of the VMI model no way has been found for making predictions for the B(E2) transition probabilities 
connecting the levels of a collective band. The su g (2) symmetry naturally provides such a link. Before 
studying the su g (2) case, though, it is useful to recall the predictions of other models on this matter. 



21.1. The collective model of Bohr and Mottelson 
In rotational bands one has III! 



B(E2 : J + 2 - J) = j^Ql\C^\ 2 , (21.1) 

i.e. the B(E2) transition probability depends on the relevant Clebsch-Gordan coefficient of su(2), while 
Qq is the intrinsic electric quadrupole moment and K is the projection of the angular momentum J on 
the symmetry axis of the nucleus in the body-fixed frame. For K — bands, as the ground state bands, 
one then has 

V ; 167r °2(2J + 3)(2J + 5) V ; 

It is clear that the B(E2) values should saturate with increasing J. 

21.2. The Interacting Boson Model (IBM) 

It is also instructive to mention what happens in the case of the Interacting Boson Model (IBM) 
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the successful algebraic model of nuclear structure with which we are going to deal in more detail 
later (see sections 27, 31). In the case of the su^ limit of the IBM, which is the limit applicable to 
deformed nuclei, the corresponding expression is 
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2 3 (J + l)(J + 2) (2N- J)(2N + J + 3) 



B(E2: J + 2^ J) = — — — —^r- >— rr , (21.3) 

V ; 16vr °2(2J + 3)(2J + 5) (2N + 3/2) 2 ' v ; 

where N is the total number of bosons. Instead of saturation one then gets a decrease of the B(E2) 
values at high J, which finally reach zero at J = 2N. This is a well known disadvantage of the simplest 
version of the model (IBM-1) due to the small number of collective bosons (s (J = 0) and d (J = 2)) 
taken into account. It can be corrected by the inclusion of higher bosons (g (J = 4), i (J = 6), etc), 



which approximately restore saturation (see 1 169 H 180 | for full list of references) 



Another way to avoid the problem of decreasing B(E2)s in the su(3) limit of IBM at high J is the 
recently proposed lllll transition from the compact su(3) algebra to the noncompact sl(3,R) algebra. The 
angular momentum at which this transition takes place is fitted to experiment. In this way an increase 
of the B(E2) values at high J is predicted, which agrees well with the experimental data for 236 U. 

21.3. The su q (2) model 

In order to derive a formula similar to (21.2) in the su g (2) case, one needs to develop an su q (2) 
angular momentum theory. As mentioned in sec. 15, this has already been achieved. It turns out 
that an equation similar to (21.1) holds in the q-deformed case, the only difference being that the 
Clebsch-Gordan coefficient of the su 9 (2) algebra must be used instead. These coefficients have the form 

E 



n J+2,2,J _ 2K 
<V^K,0,-K — 1 



[3] [4] [J - K + 2}[J - K + 1][J + K + 1][J + K + 2] 
\ [2] [2 J + 2] [2 J + 3] [2 J + 4] [2 J + 5] 



(21.4) 



For K = bands one then has 



qK ' ; 16tt °[2][2J + 2][2J + 3][2J + 4][2J + 5] V ; 



For q = e lT this equation takes the form 



B,( £ 2;J + 2^J) = A f^>^> 



167T sin(2r) sin(r) 
(sin(r(J+l))) 2 (sin(r(J + 2))) 2 



sin(r(2 J + 2)) sin(r(2 J + 3)) sin(r(2 J + 4)) sin(r(2 J + 5)) ' 



(21.6) 



It is useful to get an idea of the behaviour of this expression as a function of J, especially for the small 
values of r found appropriate for the description of ground state spectra. Expanding all functions and 
keeping corrections of the leading order in r only, one has 



Bq{E2 , J + ^J^ ^ y^^ l + ±(W + 22, + 12)). (21.7) 

We see that the extra factor, which depends on r 2 , contributes an extra increase with J, while the usual 
su(2) expression reaches saturation at high J and IBM even predicts a decrease. 

21.4- Comparison to experiment 

Is there any experimental evidence for such an increase? In order to answer this question one should 
discover cases in which the data will be consistent with the su g (2) expression but inconsistent with 
the classical su(2) expression. (The opposite cannot happen, since the classical expression is obtained 
from the quantum expression for the special parameter value r = 0.) Since error bars of B(E2) values 
are usually large, in most cases both symmetries are consistent with the data. One should expect the 
differences to show up more clearly in two cases: 

i) In rare earth nuclei not very much deformed (i.e. with an R± = E(A)/ E(2) ratio around 3.0). 
These should be deformed enough so that the su g (2) symmetry will be able to describe them, having, 
however, at the same time values of r not very small. Since in several of these nuclei backbending (or 
upbending) occurs at J = 14 or 16, one can expect only 5 or 6 experimental points to compare the 
theoretical predictions with. 

ii) In the actinide region no backbending occurs up to around J = 30, so that this is a better test 
ground for the two symmetries. However, most nuclei in this region are well deformed, so that small 
values of r should be expected, making the distinction between the two theoretical predictions difficult. 



A few characteristic examples (4 rare earths and an actinide) are given in E£3. In all cases it is clear 
that the su,j(2) curve follows the experimental points, while the su(2) curve has a different shape which 
cannot be forced to go through all the error bars. Several comments are now in place: 

i) For a given nucleus the value of the parameter r obtained from fitting the B(E2) values among 
the levels of the ground state band should be equal to the value obtained from fitting the energy levels 



of the ground state band. In 1153 it is clear that both values are similar, although in most cases the value 
obtained from the B(E2)s is smaller than the value obtained from the spectra. It should be taken into 
account, however, that in most cases the number n' of levels fitted is different (larger) than the number 
n of the B(E2) values fitted. In the single case ( 184 W) in which n = n', the two r values are almost 
identical, as they should. 

ii) One can certainly try different fitting procedures. Using the value of r obtained from the B(E2) 
values for fitting the spectrum one gets a reasonably good description of it, although the squeezing of 
the spectrum is not as much as it should have been (with the exception of 184 W). Using the value of r 
obtained from the spectrum for fitting the B(E2) values one obtains an increase more rapid than the 
one shown by the data (again with the exception of 184 W). One can also try to make an overall fit of 
spectra and B(E2)s using a common value of r. Then both the squeezing of the spectrum and the rise 
of the B(E2)s can be accounted for reasonably well although not exactly. One should notice, however, 
that the experimental uncertainties of the B(E2)s are much higher than the uncertainties of the energy 
levels. 

iii) Concerning energy levels, the rigid rotator model and the su(3) limit of the IBM predict a J( J+l) 
increase, while the su ? (2) model and the VMI model predict squeezing, which is seen experimentally. 

iv) Concerning the B(E2) values, the VMI makes no prediction, the rigid rotator predicts saturation 
at high J, the su(3) limit of the IBM predicts decrease, while the su ? (2) model predicts an increase. The 
evidence presented in 1153 supports the su 5 (2) prediction, but clearly much more work, both experimental 
and analytical is needed before final conclusions can be drawn. The modified su(3) limit of IBM 
described in cM, as well as the su(3) limit of the sdg-IBM 1123 also support the increase of the B(E2) 
values at high J. Increasing BE(2) values are also predicted i n th e framework of the Fermion Dynamical 



Symmetry Model II2 4 ], as well as by a binary cluster model t 8 ^. There is also empirical evidence for 



increasing B(E2) values in the recent systematics by Zamfir and Casten 
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v) It is clear that much further work is needed as far as comparisons of the su g (2) predictions to 
experimental BE(2) values are concerned for safe conclusions to be reached. 

vi) Since the quadrupole operator is not a member of the symmetry algebra su 9 (2) under consider- 
ation, it is clear that the B(E2) values studied here do not contain any dynamical deformation effects, 
but only the kinematical ones (through the use of the g-deformed Clebsch-Gordan coefficients). A more 
complete approach to the problem will be the construction of a larger algebra, of which the quadrupole 
operator will be a member and it will also be an irreducible tensor under su g (2) or so 9 (3). This task 
will be undertaken in sec. 29. 



22. Superdeformed bands 

One of the most impressive experimental discoveries in nuclear physics during the last decade was 
that of superdeformation (see tH^HE^ for relevant reviews). The energy levels of superdeformed 
bands follow the J(J + 1) rule much more closely than the usual rotational bands. Levels with J larger 



than 60 have been observed. A compilation of superdeformed bands has been given in E£il. The best 
examples have been found in the A^i50 mass region, while additional examples have been found in 
the A^80, A^130 and A^190 regions. It is understood that the superdeformed bands in the A^150 
region correspond to elongated ellipsoidal shapes with an axis ratio close to 2:1, while in the A^130 
and Af»190 regions the ratios 3:2 and 1.65:1 respectively appear closer to reality. 

Since the su g (2) model has been found suitable for describing normal deformed bands, it is plausible 
that it will also be successful in describing superdeformed bands as well. A test has been performed in 
The su g (2) model has been found to give good results in all mass regions, the deformation parameter 
r being smaller than in the case of normal deformed bands, thus indicating smaller deviations from the 
su 9 (2) symmetry. In particular, r has been found to obtain values about 0.01 in the A?»130 and A^190 
regions, while it obtains even smaller values, around 0.004, in the A?»150 region, which contains the 
best examples of superdeformed bands observed so far. These results should be compared to the values 
of r around 0.03 obtained in the case of normal deformations. 

Concerning the corresponding B(E2) values, the experimental information is still quite poor for 
allowing a meaningful comparison of the su g (2) predictions to experiment. 

A new effect, called AI = 2 staggering, has been recently observed in a few superdeformed bands 
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I9j,|i9^|i96[|i97| jf /\j = 2 staggering is present, then, for example, the levels with / = 2, 6, 10, 14, 
. . . are displaced relatively to the levels with / = 0, 4, 8, 12, . . . , i.e. the level with angular momentum 
I is displaced relatively to its neighbours with angular momentum / ± 2. In most cases, though, the 
effect is smaller tha n t h e experimental error. AI = 2 staggering has also been seen in some normal 
deformed bands t£^l'£23'E29. Although many theoretical attempts for the explanation of the effect have 
been made (see E2il for a, r el evan t list), it seems plausible that what is seen in most cases is simply a 
series of badcbendings EMEH The same effect has been observed in rotational bands of diatomic 
molecules c23. In this case the effect is clearly larger than the experimental error and it is also clear 
that it is due to interba nd int eractions E2II. Furthermore in diatomic molecules a AI = 1 staggering 
effect has been observed r 05 tr 06 l, i.e. a relative displacement of the levels with even I with respect to the 
levels with odd /. This effect is clearly larger than the experimental error. A description of the AI = 1 
staggering effect in terms of two g-deformed rotators with slightly different parameter values has been 
attempted 
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23. The physical content of the su g (2) model 

From the above it is clear that the su g (2) model offers a way of describing nuclear stretching, i.e. 
the departure of deformed nuclei from the su(2) symmetry of the rigid rotator, similar to the one of 
the VMI model. The parameter r describes this departure quantitatively, vanishing in the rigid rotator 
limit. Therefore the deformation parameter r should not be confused with nuclear deformation; it is in 
fact related to nuclear softness, as already discussed in sec. 20. 

On the other hand, the increase of the moment of inertia with increasing J in the framework of 



the VMI model means that collectivity gets increased EZa. The su g (2) model is an alternative way for 



describing this increase in collectivity. But increased collectivity implies increased B(E2) transitions. 



Therefore it is not surprising that the su g (2) model predicts B(E2) values increasing with J. 

Given the su 5 (2) generators J + , J_, Jo, it is instructive to define as usual the operators J x , J y , J z 

by 

J+ = Jx + i Jy, J- = J x — iJy, Jo = J z - (23.1) 
The su 9 (2) commutation relations can then be rewritten in the form 

i 

\Jxi Jy\ 2 P J ^ ' \Jy> J?\ iJx) \Jzi Jx\ iJy, (23.2) 

which is a generalization of the so (3) commutation relations, obtained in the limit q — > 1. We remark 
that while in the classical so(3) case the three commutation relations have exactly the same form, in the 
quantum case the first commutation relation differs (in the right hand side) from the other two, thus 
indicating that in the framework of the problem under study the z-direction is not any more equivalent 
to the x- and y- directions. This is of course a phenomenological way to describe the softness of deformed 
nuclei by adding the appropriate perturbations to the pure su(2) Hamiltonian and has nothing to do 
with the isotropy of space, as implied in 
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It is worth remarking at this point that the Casimir operator of su 9 (2) is also invariant under 
the usual su(2) lH2a. Therefore the quantum number J characterizing the irreps of su g (2), and as a 
result the nuclear levels, is exactly the same as the quantum number J used in the case of the usual 
su( 2). Therefore there is no reason for a "total reformulation of quantum mechanics", as implied in 
the su ? (2) generators being connected to their su(2) counterparts by the g-deforming functionals 
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of sec. 14. In other words, one continues to believe in usual angular momentum theory and usual 
quantum mechanics. All what is done in the framework of the su g (2) model is to add to the usual su(2) 
Hamiltonian several perturbations which have a special form making them suitable to be summed up, 
including the original su(2) term, into the form of the su 9 (2) Hamiltonian. 

24. The u p(? (2) rotator model 



An extension of the su 9 (2) model is the u Pj(? (2) model E2j'E13'Eij'Bl3, which is based on a two-parameter 
deformed algebra (see sec. 15 for a list of references). For p = q (using the definition of (p, g)-numbers 
of eq. (2.11)) this model reduces to the su g (2) one. This model has been successfully applied to 
superdeformed nuclear bands When Taylor expanded, it becomes clear that the eigenvalues of 

the Hamiltonian of this model (which is the second order Casimir operator of u Pi5 (2)) contain terms of 
the form J(J(J + l)) n , in addition to the (</(</ + l)) n ones. It is therefore closer to the modification of 
su g (2) which will be discussed in subsec. 26.3. 

25. Generalized deformed su(2) models 

Another formula giving very good results for rotational spectra has been introduced by Holmberg 
and Lipas EH, and rediscovered in Ell. In this case the energy levels are given by 



E(J) = a 



l + bJ(J + 1) - 1 



(25.1) 



Taylor expansion of the square root immediately shows that the present formula is a special case of eq. 
(19.6). This formula can be derived from the collective model of Bohr and Mottelson tla. 



It has been argued in 11^3 that the Hamiltonian of eq. (25.1) gives better agreement to rotational 
nuclear spectra than the one coming from the su g (2) symmetry. Using the techniques described in 
detail in sec. 17 one can construct a generalized deformed algebra su$(2), characterized by a function 
$( J( J+ 1), giving the spectrum of eq. (25.1) exactly. In this particular case the algebra is characterized 
by the structure function 

$(J(J+1)) 



1 + bJ(J+ 1) - 1 



(25.2) 



It is of interest to check if this choice of structure function also improves the agreement between theory 
and experiment in the case of the electromagnetic transition probabilities connecting these energy levels. 



In order to study this problem, one has to construct the relevant generalized Clebsch-Gordan coefficients. 
This problem is still open. 

26. Quantum algebraic description of vibrational and transitional nuclear spectra 

We have already seen that the su g (2) model describes successfully deformed and superdeformed 
bands. It is not surprising that the applicability of the su g (2) formalism is limited to the rotational 
region (where the ratio R4 = E(4)/E(2) obtains values between 3.0 and 3.33), since it is based on a 
deformation of the rotation algebra. For describing nuclear spectra in the vibrational (2.0 < R4 < 2.4) 
and transitional (2.4 < R4 < 3.0) regions it is clear that an extension of the model is needed. In order 
to be guided towards such an extension, we briefly review the existing experience of other successful 
models. 

26.1. The Interacting Boson Model 

In the rotational (su(3)) limit of the Interacting Boson Model (IBM) (see sees 27, 31 for more 
details) the spectrum is described by a J(J + 1) expression, while in the vibrational (u(5)) cl3 and 
transitional (0(6)) Bl3 limits expressions of the form J(J + c) with c > 1 appear. In the u(5) limit, in 
particular, the energy levels are given by 

E(N, n d , v, n A , J, Mj) = E + en d + an d (n d + 4) + ft2v(v + 3) + 7 2 J( J + 1), (26.1) 

where N is the total number of bosons, n d is the number of d-bosons, v is the seniority, n& is the 
"missing" quantum number in the reduction from o(5) to o(3), Mj is the third component of the 
angular momentum J, while E , e, a, ft, 7 are free parameters. The ground state band, in particular, 
is characterized by quantum numbers n d = 0, 1, 2, . . . , v = n d , n A = 0, J = 2n d , so that the energy 
expression for it reads 

E(J) = E + € -J+ jJ{J + 8) + I J( J + 6) + 2 7 J( J + 1). (26.2) 

In the 0(6) limit the energy is given by 

E(N, a, t, u A , J, Mj) = E + (32t(t + 3) + 7 2J( J + 1) + r]2a(a + 4), (26.3) 

where a is the quantum number characterizing the irreducible representations (irreps) of 0(6), r is the 
quantum number characterizing the irreps of o(5), v A is the missing quantum number in the reduction 
from o(5) to o(3), while E , ft, 7, r\ are free parameters. The ground state band is characterized by the 
quantum numbers a — N, r = 0, 1, 2, . . . , u A — 0, J — 2r, so that the relevant energy expression takes 
the form 

E{J) =E + ^J(J + 6) + 7 2J( J + 1) + r/2N(N + 4). (26.4) 

The message from eqs (26.2) and (26.4) is that nuclear anharmonicities are described by expressions 
in which J and J 2 appear with different coefficients, and not with the same coefficient as in J(J + 1). 
The earliest introduction of this idea is in fact the Ejiri formula EiZ 

E(J) = k.J(J + 1) + aJ, (26.5) 

which has been subsequently justified microscopically in £12. 

26.2. Generalized VMI 



The two-parameter VMI model is known to continue giving good fits in the transitional and even in 
the vibrational region. In these regions, however, the accuracy of the model is substantially improved 
hy .a d di n g a third parameter, which essentially allows for treating J and J 2 with a different coefficient 



The usual VMI model has been briefly reviewed in subsec. 20.1. One of the (essentially equivalent) 
three-parameter extensions of the model, which give improved fits of vibrational and transitional spectra, 
is the generalized VMI (GVMI) model ELSE!] in which the energy levels are described by 

E(J) = J + ; J ( ( ^~ 2) + - $ ) 2 , (26.6) 

which can be easily rewritten in the form 

E{J) = 4i^7p + ^' ($,(J) ^ $ ° )2 ' (26 ' 7) 

where x' = x^ 1 — 2. It is clear that for x = 1/3 the GVMI reduces to the usual VMI, while for 
transitional and vibrational nuclei x obtains lower values, so that x' becomes greater than 1. The 
variational condition determining the moment of inertia still has the form of eq. (20.4), while the 
expansion of the energy turns out to be 

E{J) = i( J(J + *') - °-{J{J + x')f + (a') 2 (J(J + *')) 3 - 3(a') 3 (J(J + x')f + ...), (26.8) 

where 

°' = wW- (26 ' 9) 

We remark that an expansion in terms of J(J + x') is obtained, as compared to an expansion in terms 
of J(J + 1) in the case of the usual VMI. The physical content of the parameters is clear: the centrifugal 
stretching effect is accounted for by the softness parameter a', as in the case of the usual VMI, while 
anharmonicities, important in the vibrational region, are introduced by i' > 1. Since centrifugal 
stretching and anharmonicities are two effects of different origins, it is reasonable to describe them by 
two different parameters. 



26.3. Modification of the su q (2) model 

The evidence coming from the IBM and the generalized VMI model described above, suggests a 
model in which the spectrum is given by 

E(J) = jj[JW + c] q , (26.10) 

which contains 3 parameters: the moment of inertia /, the deformation parameter q and the new 
parameter c, which is expected to be 1 in the rotational limit and larger than 1 in the vibrational and 
transitional regions. This energy expression can be expanded as 

E(J) = YlJJ^rW^ {T)J{J + C) " T ^ J ( J + c )) 2 

+^ 2 j 2 (r)(J(J + c)) 3 - lr 3 j 3 (r)(J(J + c)) 4 + ^r%(r)(J(J + c)) 5 -...), (26.11) 

which is similar to eq. (19.5) with J(J + 1) replaced by J(J + c). 

It is expected that the deformation parameter r, which plays the role of the small parameter in 
the expansion, as the softness parameter does in the case of the VMI, will describe the centrifugal 
stretching effect, while the parameter c will correspond to the an harm on i city effects. These expectations 
are corroborated from fits of the experimental data reported in The following comments can be 
made: 

i) The anharmonicity parameter c is clearly decreasing with increasing i? 4 , i.e. with increasing 
collectivity. It obtains high values (8-18) in the vibrational region, while in the rotational region it 
stays close to 1. (It should be noted that by fixing c = 1 in the rotational region the fits are only very 
slightly changed, as expected.) In the transitional region its values are close to 3. 



ii) The deformation parameter r, which corresponds to the centrifugal stretching, is known from 
the su g (2) model to obtain values close to 0.3-0.4 in the rotational region, a fact also seen here. The 
same range of values appears in the vibrational region as well, while in the transitional region r reaches 
values as high as 0.6. It is not unreasonable for this parameter, which is connected to the softness of 
the nucleus, to obtain its highest values in the region of 7-soft nuclei. 

iii) It is worth remarking that eq. (26.10) coincides for q = 1 and c= integer =Nwith the eigenvalues 
of the Casimir operator of the algebra o(N+2) in completely symmetric states EH3£H3. In the rotational 
region the fits gave N=l, which corresponds to o(3), as expected, while in the transitional region the 
fits gave approximately N=3, which corresponds to o(5), which is a subalgebra contained in both the 
u(5) and 0(6) limits of the IBM. A unified description of 927 low lying states of 271 nuclei from five 
major shells along these lines has been given in c2a. 

iv) It is also worth remarking that a special case of the expression of eq. (26.10) occurs in the 
q-deformed version of the 0(6) limit of the Interacting Boson Model, which will be reported below (see 
eq. (31.7)). 

v) The su q (2) symmetry is known to make specific predictions for the deviation of the behaviour 
of the B(E2) values from the rigid rotator model (subsec. 21.3). It will be interesting to connect the 
spectrum of eq. (26.10) to some deformed symmetry, at least for special values of c, and examine the 
implications of such a symmetry for the B(E2) values. Such a study in the framework of the q-deformed 
version of the 0(6) limit of IBM, mentioned in iv), is also of interest. 

vi) It is worth noticing that an expansion in terms of J( J+c) can also be obtained from a generalized 
oscillator (sec. 12) with a structure function 

F(J) = [J(J + c)] Q , (26.12) 

where [x]q stands for the Q-numbers introduced in sec. 6 and Q = e T , with T real. This is similar to 
an oscillator successfully used for the description of vibrational spectra of diatomic molecules (see sec. 
37). It can also be considered as a deformation of the oscillator corresponding to the Morse potential 
(see sec. 37). 

We have therefore introduced an extension of the su g (2) model of rotational nuclear spectra, which 
is applicable in the vibrational and transitional regions as well. This extension is in agreement with the 
Interacting Boson Model and the Generalized Variable Moment of Inertia model. In addition to the 
overall scale parameter, the model contains two parameters, one related to the centrifugal stretching 
and another related to nuclear anharmonicities. In the rotational region the model coincides with the 
usual sUq(2) model, while in the transitional region an approximate o(5) symmetry is seen. These results 
give additional motivation in pursueing the construction of a deformed version of the Interacting Boson 
Model. This problem will be discussed in the next section as well as in sec. 31. 



27. A toy Interacting Boson Model with su,j(3) symmetry 

The Interacting Boson Model (IBM) EH is a very popular algebraic model of nuclear structure. In 
the simplest version of IBM low lying collective nuclear spectra are described in terms of s ( J = 0) and 
d (J = 2) bosons, which are supposed to be correlated fermion pairs. The symmetry of the simplest 
version of the model is u(6), which contains u(5) (vibrational), su(3) (rotational) and 0(6) (7-unstable) 
chains of subalgebras (see also sec. 31). A simplified version of the model, having the su(3) symmetry 
with su(2) and so(3) chains of subalgebras also exists 0. It can be considered as a toy model for 
two-dimensional nuclei, but it is very useful in demonstrating the basic techniques used in the full IBM. 

In the present section we will construct the g-deformed version of this toy model. Since this project 
requires the construction of a realization of su g (3) in terms of g-deformed bosons, we will also use this 
opportunity to study su g (3) in some detail. 



27.1. The su q (3) algebra 

In the classical version of the toy IBM IhI one introduces bosons with angular momentum m = 0, ±2, 
represented by the creation (annihilation) operators ). They satisfy usual boson 

commutation relations 

[ai,a]} = 5ij, [ai,aj] = [a], a]} = 0. (27.1) 



The 9 bilinear operators 

Ay = a\ aj (27.2) 

satisfy then the commutation relations 

[A i:j , A M ] = 5 jk A u - 5 u A kj , (27.3) 

which are the standard u(3) commutation relations. The total number of bosons 

N = = a\a + a ] + a + + ata_ (27.4) 

is kept constant. Since we are dealing with a system of bosons, only the totally symmetric irreps 
{N, 0, 0} of u(3) occur. 

In the quantum case one has the u 9 (3) commutation relations given in Table 2 EH3, where are 
the generators of u g (3) and the g-commutator is defined as 

[A, B} q = AB - qBA. (27.5) 



Table 2: u 9 (3) commutation relations 210 , given in the form [A, B] a = C. A is given in the first column, 
B in the first row. C is given at the intersection of the row containing A with the column containing 
B. a, when different from 1, follows C, enclosed in parentheses. 
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In order to obtain a realization of u g (3) in terms of the g-bosons described in sec. 10, one starts 
with EH 

A12 = a\a 2 , A 2 i = a\a u A 23 = a\a 3 , A 32 = a\a 2 . (27.6) 

One can easily verify that the u g (3) commutation relations involving these generators are satisfied. For 
example, one has 

[A12, A 21 ] = [ATi - N 2 ], [A 23 , A 32 ] = [N 2 - N 3 ], (27.7) 
using the identity of eq. (2.6) and the identifications 

N 1 = An, N 2 = A 22 , N 3 = A 33 . (27.8) 



One can now determine the boson realizations of A13 and A31 from other commutation relations, as 
follows 

A 13 = [A 12 , A 23 ] q = a\a 3 q- N \ (27.9) 



A 3 i = [A 32 , A 21 ] q -i 
Using eq. (2.6) once more one can verify that the relation 

[A 13 ,Az 1 ] = [N 1 -N 3 



a 3 aiq . 



(27.10) 



(27.11) 



is fulfilled by the boson images of (27.9), (27.10). It is by now a straightforward task to verify that all 
commutation relations of Table 2 are fulfilled by the boson images obtained above. 

Before turning to the study of the two limits of the model, we give for completeness some additional 
information on su g (3): 

i) The irreps of su g (3) have been studied in 




ii) Clebsch-Gordan coefficients for su g (3) have been given in 

iii) The Casimir operator s of S%(3 ) arid their eigenvalues have been given explicitly in 



the Elliott quantum numbers 
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A = fx - f: 
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V = f: 
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Using 
(27.12) 



where fi represents the number of boxes in the i-th line of the corresponding Young diagram, the irreps 
of su 9 (3) are labelled as (A,/i). The eigenvalues of the second order Casimir operator then read 



+ 



2A fi 1 

h - + 1 

3 3 



+ 



A 2/i 1 
- + — + 1 
3 3 



(27.13) 



while the eigenvalues of the third order Casimir operator are 

C 3 = 2 



2A n 1 

h - + 1 

3 3 



A 2fi i 
- + — + 1 

3 3 



(27.14) 



In the limit q — > 1 these reproduce the ordinary results for su(3): 

2 



C 2 = -(A 2 + // + A// + 3A + 3aO, 



c 3 
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(A-/i)(2A + /i + 3)(A + 2 y u + 3). 



(27.15) 
(27.16) 



2 if 



iv) Irrducible tensor operators under u g (3) have been considered in 

v) A different deformation of sl(3) has been studied in Bl^. 

vi) Coupling (Wigner-Clebsch-Gordan) coefficients for the algebras u g (n) have been considered in 



27.2. The su q (2) limit 

We shall study the su g (2) limit of the model first, since it is technically less demanding. 

So far we have managed to write a boson realization of u g (3) in terms of 3 g-bosons, namely a\, a 2 , 
a 3 . Omitting the generators involving one of the bosons, one is left with an su g (2) subalgebra. Omitting 
the generators involving a 3 , for example, one is left with A 12 , A 2 i, N x , N 2 , which satisfy usual su 3 (2) 
commutation relations if the identifications 

J + = A 12 , J- = A 2U J = i(JVi-JV 2 ) (27.17) 

are made. Jo alone forms an so q (2) subalgebra. Therefore the relevant chain of subalgebras is 



su g (3) D su g (2) D so,(2). 



(27.18) 



The second order Casimir operator of su g (2) has been given in eq. (14.9). Substituting the above 
expressions for the generators one finds 



C 2 (su 9 (2)) 



N± + N 2 



+ 1 



(27.19) 



All of the above equations go to their classical count erpa rts by allowing q — > 1, for which [x 



o-numbers become usual numbers. In the classical case EHSl out of the three bosons ( ) forming 



su(3), one chooses to leave out ao, the boson with zero angular momentum, in order to be left with the 
su(2) subalgebra formed by a + and a_, the two bosons of angular momentum two. The choice of the 
su q (2) subalgebra made above is then consistent with the following correspondence between classical 
bosons and g-bosons 

a + — > a\, a_ — > a 2 , an ~~ > a 3- (27.20) 

(We have opted in using different indices for usual bosons and g-bosons in order to avoid confusion.) 

In the classical case the states of the system are characterized by the quantum numbers characterizing 
the irreducible representations (irreps) of the algebras appearing in the classical counterpart of the chain 
of eq. (27.18). For su(3) the total number of bosons iV is used. For su(2) and so(2) one can use the 
eigenvalues of J 2 and J , or, equivalently, the eigenvalues of a\a + + aLa_ and L 3 = 4J , for which we 
use the symbols n d (the number of bosons with angular momentum 2) and M. Then the basis in the 
classical case can be written as @ 



\N, n d ,M>-- 



J\N-n d I t \n d /2+M/4: 



,t W2-M/4 



(N - n d )\ (n d /2 + M/4)! (n d /2 - M/4)! 



10 > . 



(27.21) 



In the quantum case for each oscillator one defines the basis as in sec. 10. Then the full basis in the 
g-deformed case is 



\N,n d ,M>, 



^n d /2+M/4 



(a t)n d /2-M/4 



[N - n d }\ [ml 2 + M/4]! [n d /2 - M/4]! 



10 >, 



(27.22) 



where N = N% + N 2 + N 3 is the total number of bosons, n d = Ni + N 2 is the number of bosons with 
angular momentum 2, and M is the eigenvalue of L = 4J . n d takes values from up to N, while for a 
given value of n d , M takes the values ±2n d , ±2(n d — 2), . . . , ±2 or 0, depending on whether n d is odd 
or even. In this basis the eigenvalues of the second order Casimir operator of su g (2) are then 



C 2 (su ? (2))|iV, n d ,M > ? 



n d 



+ 1 



\N,n d ,M>„ . 



(27.23) 



In the case of N = 5 one can easily see that the spectrum will be composed by the ground state band, 
consisting of states with M = 0, 2, 4, 6, 8, 10 and n d = M/2, the first excited band with states 
characterized by M = 0, 2, 4, 6 and n d = M/2 + 2, and the second excited band, containing states with 
M = 0, 2 and n d = M/2 + 4. 

In the case that the Hamiltonian has the su g (2) dynamical symmetry, it can be written in terms of 
the Casimir operators of the chain (27.18). Then one has 



H = E + AC 2 (su q (2)) + £?C 2 (so 9 (2)), 
where E , A, B are constants. Its eigenvalues are 



(27.24) 



E = E + A 



n d 



n d 



+ 1 



+ BM 2 



(27.25) 



Realistic nuclear spectra are characterized by strong electric quadrupole transitions among the levels 
of the same band, as well as by interband transitions. In the framework of the present toy model one 
can define, by analogy to the classical case EH^I, quadrupole transition operators 



Q- 



a[a 3 + a\a 2 , 



Q- 



a 2 a 3 + a^ai. 



(27.26) 



In order to calculate transition matrix elements of these operators one only needs eqs. (10.7), (10.8), 
i.e. the action of the g-boson operators on the g-deformed basis. The selection rules, as in the classical 
case, are AM = ±2, An d = ±1, while the corresponding matrix elements are 



<N,n d + l,M± 2\Q±\N, n d , M > c 



[N-n d ][-f±- + l], 



(27.27) 



< N,n d -l,M±2\Q±\N,n d ,M> q = ^ [N - n d + 1][^ t ^] 



(27.28) 



From these equations it is clear that both intraband and interband transitions are possible. 

In order to get a feeling of the qualitative changes in the spectrum and the transition matrix elements 
resulting from the g-deformation of the model, a simple calculation for a system of 20 bosons (N = 20) 
has been performed in E£D. Two cases are distinguished: i) q real (g = e T , with r real), ii) q a phase 
factor (q = e lT , with r real). The main conclusions are: 

i) When q is real the spectrum is increasing more rapidly than in the classical case, while when q is 
a phase the spectrum increases more slowly than in the classical case, in agreement with the findings 
of the g-rotator model. 

ii) The transition matrix elements in the case that q is real increase more rapidly than in the classical 
case, while they increase less rapidly than in the classical case when q is a phase. 

iii) Transition matrix elements are much more sensitive to g-deformation than energy spectra. This 
is an interesting feature, showing that g-deformed algebraic models can be much more flexible in the 
description of transition probabilities than their classical counterparts. 



27.3. The so q (3) limit 

The classical su(3) toy model has, in addition to the above mentioned su(2) chain of subalgebras, an 
so(3) chain. However, the problem of constructing the su g (3)DsOg(3) decomposition is a very difficult 
one. Since this decomposition is needed in constructing the g-deformed versions of several collective 
models, including the Elliott model BUSS, the su(3) limit of the IBM SlH, the Fermion 
Dynamical Symmetry Model (FDSM) El i], the Interacting Vector Boson Model (IVBMlil'Ei'Ei'El, 
the nuclear vibron model for clustering EH, as well as the su(3) limit of the vibron model HI for diatomic 
molecules, we report here the state of the art in this problem: 

i) As far as th e co m pletely s y mmetric irreps of su g (3) are concerned, the problem has been solved 
by Van der Jeugt M&^Wm. This suffices for our needs in the framework of the toy IBM, since 
only completely 

ii) Sciarrino 

However, it was not clear how to impose the Hopf structure on this larger algebra. Trying the other way 
around, he found that by starting from a gig (3) algebra, which already possesses the Hopf structure, 
one loses the Hopf structure of the principal 3-dim subalgebra, which should have been so g (3). 

3 attacked the problem through the use of g-deforming 



iii) Pan E5a and Del Sol Mesa et al. 
functionals of see s 10, 14. 

iv) Quesne e£3 started with g-bosonic operators transforming as vectors under so g (3) and constructed 
a g-deformed u(3) by tensor coupling. 

v) g-deformed subalgebras of several g-deformed algebras have recently been studied by Sciarrino 



/mmetric su g (3) irreps occur in it. 
^3 started from so g (3) and obtained a deformed gl(3) containing so g (3) as a subalgebra. 



vi) A simplified version of the so 9 (3) subalgebra of u 9 (3) has been constructed in Furthermore, 
explicit expressions for the irreducible vector and quadrupole tensor operators under so g (3) have been 
given and the matrix elements of the latter have been calculated. Further details on these developments 
will be given in sections 28 and 29. 



In what follows, it suffices to use the solution given in E3, since in the model under study only 



completely symmetric irreps of u 9 (3) enter. Using the notation 



Ct_|_ > 0\, 



a_ 



0,2, Oq 



0-3, 



(27.29) 



the basis states are of the form 

(aL ) n+ (al) no (aL) n - , 

|n+n n_ >= v + / v 0J y =^\0 >, (27.30) 

r n+]![n ]![n_]! 

with cij|0 >= and Ni\n + riQn- >= ni\n + non- >, where % — +, 0, — . The principal subalgebra so g (3) is 
then generated by 

L = N + -N_, (27.31) 

L + = q N -? No y/q N ++q- N + a ] + a + aja„ q N ^ N ^q N - + q~ N -, (27.32) 



L_ = ala + q N —? No y/q N + + q- N + + q N +-^ N °^q N - + q- N -ala , (27.33) 
satisfying the commutation relations 

[Lo, L ± ] = ±L±, [L+, L_] = [2L ]. (27.34) 

L alone generates then the so g (2) subalgebra. Therefore the relevant chain of subalgebras is 

su g (3) D so 9 (3) D bo,(2). (27.35) 

The so g (3) basis vectors can be written in terms of the vectors of eq. (27.30) as 



\v(N,L,M) >= q -l(L+M)(L+M-W\ 



[N + L]\\[2L + 1][L + M]\[L - M]\ 
\ [N - L]\\[N + L + 1]\ 



y^ q (2L-i)x/2 s (n-l)/2 \x, L + M — 2x, x — M > (27.36) 
^' / [2a;]!![L + M - 2x]![2x - 2M}\\ ' 



where 



[alfq N++N - +l 



[2iV + ][2iV_j t t „„ 



a\aLq (27.37) 



x takes values from max(0,M) to [(L + M) /2] in steps of 1, L = N, N- 2, . . . , 1 or 0, M = -L, -L + 1, 
. . . , +L, and [2x]\\ = [2x][2x — 2} . . . [2]. The action of the generators of so q (3) on these states is given 
by 

L \v(N, L, M) >= M \v(N, L, M) >, (27.38) 



L± \v(N, L, M) >= sJ[L =f M] [L ± M + 1] \v(N,L,M ± 1) > . (27.39) 
The second order Casimir operator of so g (3) has the form 

C 2 (so 9 (3)) = L 2 = L_L + + [L ][L + 1]. (27.40) 

Its eigenvalues in the above basis are given by 

C 2 (so 9 (3)) \v(N, L, M) >= [L] [L + 1] \v(N, L,M)> . (27.41) 

All of the above equations go to their classical count erp arts by allowing q — > 1, for which [x] — > x, i.e. 
g-numbers become usual numbers. In the classical case c3B the states of the system are characterized by 
the quantum numbers characterizing the irreducible representations (irreps) of the algebras appearing 
in the classical counterpart of the chain of eq. (27.35). For sn(3 ) the total number of bosons iV is used. 
For so (3) and so (2) one can use L and M, respectively. In cHa, however, the eigenvalue of L' = 2Lq is 
used, which is M' = 2M. 

Since the rules for the decomposition of the totally symmetric u g (3) irreps into so g (3) irreps are the 
same as in the classical case, it is easy to verify that for a system with N = 6 the spectrum will be 
composed by the ground state band, consisting of states with M' = 0, 2, 4, 6, 8, 10, 12 and L = N, the 
first excited band with states characterized by M' = 0, 2, 4, 6, 8 and L = N — 2, the second excited 



band, containing states with M = 0, 2, 4 and L = N — 4, and the third excited band, containing a state 
with M' = and L = N - 6. 

In the case that the Hamiltonian has the so g (3) dynamical symmetry, it can be written in terms of 
the Casimir operators of the chain (27.35). Then one has 



H = E + AC 2 (so 9 (3)) + BC 2 {so q {2)), 
where E , A, B are constants. Its eigenvalues are 

E = E + A[L][L + 1] + BM 2 . 



(27.42) 



(27.43) 



It is then clear that in this simple model the internal structure of the rotational bands is not influenced 
by g-deformation. What is changed is the position of the bandheads. 

We turn now to the study of electromagnetic transitions. In the present limit one can define, by 
analogy to the classical case, quadrupole transition operators Q± proportional to the so g (3) generators 
L ± 



(27.44) 



In order to calculate transition matrix elements of these operators one only needs eq. (27.39). The 
selection rules, as in the classical case, are AM' = ±2, AL = 0, i.e. only intraband transitions are 
allowed. The relevant matrix elements are 



< v(N, L, M' + 2)\Q + \v{N, L, M') > 



< v(N,L } M' -2)\Q-\v(N,L,M') >- 




L + 



M' 



M' 



+ 1 



(27.45) 



(27.46) 



In order to get a feeling of the qualitative changes in the spectrum and the transition matrix elements 
resulting from the g-deformation of this limit of the model, a simple calculation has been done in . 



Again the cases of q being real or q being a phase factor have been considered. The main conclusions 
are: 

i) g-deformation influences only the position of bandheads, while it leaves the internal structure of 
the bands intact. 

ii) When q is real the bandheads are increasing more rapidly than in the classical case, while when 
q is a phase the bandheads increase more slowly than in the classical case. This result is in qualitative 
agreement with the findings of the su g (2) model. 

iii) Transition matrix elements in the case that q is real have values higher than in the classical case, 
while they have values lower than in the classical case when q is a phase. 



28. The 3-dimensional g-deformed harmonic oscillator and the nuclear shell model 

Having developed in the previous section several techniques related to the u 9 (3) algebra and its 
so g (3) subalgebra, it is instructive to see how a g-deformed version of the 3-dimensional harmonic 
oscillator possessing this symmetry can be built and how this oscillator can be used for reproducing the 
predictions of the Modified Harmonic Oscillator, first suggested by Nilsson £53 and studied in detail in 
2,1 -2, i 272| rj^g cons ^ ruc ^i on Q f t;h e Hamiltonian of the 3-dimensional g-deformed harmonic oscillator is 



a non-trivial problem, since one has to construct first the square of the g-deformed angular momentum 
operator. 



28.1. Simplified so q (3) subalgebra and so q (3) basis states 

The so g (3) subalgebra of u g (3) has already been mentioned in subsec. 27.3. Here we shall derive a 
simplified form of it. 

We shall use three independent g-deformed boson operators 6j and b\ (i = +,0,—), which satisfy 
the commutation relations 



[NiM] [N h h] 



hb\ - q ±x b\b, = q* Ni 



(28.1) 



where iVj are the corresponding number operators. (These bosons are the same as the bosons a and a', 
with i — +, 0, — , used in the previous section.) 

The elements of the so 9 (3) algebra (i.e. the angular momentum operators) acting in the Fock space 
of the totally symmetric representations [N,0,0] of u 9 (3) (which is in fact the space of the 3-dimensional 
g-deformed harmonic oscillator) have been derived by Van der Jeugt HHHl^li^li^ll. As it has already 
been mentioned in the previous section^ a simplified form of the subalgebra so g (3) C u g (3), can be 
obtained by introducing the operators 



2(57 



B = q-* No b Q , 



Bl 



'69 



jiVo 



Bi 



B 



\ [Ni] 



+,-, 



(28.2) 



where the g-numbers [x] are defined as in eq. (2.1). We shall consider only real values for the deformation 
parameter q (i.e. q = e T with r being real). These new operators satisfy the usual commutation relations 



\^B\ 



B] 



[Ni,By 



-B, 



(28.3) 



while in the Fock space, spanned on the normalized eigenvectors of the excitation number operators 
N + , N , N_, they satisfy the relations 



BlB = q- N ° +1 [N ] 



b\b 1 



,2JVi-lr 



q-' l [2Ni] 
and hence, the commutation relations 

[B ,Bl] = q 



B l B i 



B Bl = q- No [N + l}, 
q 2N * +l [2N t + 2\, i = +, 



2N 



[Bi, Bj] = [2}q 



+, 



(28.4) 
(28.5) 

(28.6) 



It was shown in that the angular momentum operators defined in and reported in eqs (27.31)- 
(27.33), when expressed in terms of the modified operators of eq. (28.2) , take the simplified form 

L = N + - AL, L + = q- Lo+ ^B j + B + q L °-^B j B_, L_ = q~ L °-^B j B + + q Lo+ ^B j _B , (28.7) 

and satisfy the standard so q (3) commutation relations 

[Lq, L ± ] = ±L±, [L+, L_] = [2L ). (28.8) 

As it has been discussed in detail in the algebra of eq. (28.7) is not a subalgebra of u g (3) when these 
algebras are considered as Hopf algebras, but it is a subalgebra when these algebras are considered as 
g-deformed enveloping algebras. However, in our context, working in the Fock space (constructed from 
three independent g-bosons), we restrict ourselves to the symmetric representation of u 9 (3), where the 
embedding so g (3) C u g (3) is valid. It should also be noted that the operators L+, L__, L given by eq. 
(28.7) are not invariant with respect to the replacement q <-> q~ l a fact that restricts us to real q. 
The Casimir operator of so g (3) can be written in the form lEiME 6 *" 1 
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= l - {L + L_ + + [2][L ] 2 } = + [L ][L + 1] = + [L ][L - 1]. 

One can also define g-deformed states \n£m) q satisfying the eigenvalue equations 

q ~L 2 \n£m) q = [£}[£+ l]\n£m) q , L \n£m) q = m\n£m) q , N\n£m) q = n\n£m) q . 



(28.9) 



(28.10) 



Here q L 2 = f C ( * ] is the "square" of the g-deformed angular momentum and = N < 4 - A'q + N_ is the 
total number operator for the g-deformed bosons. These states have the form (see e£3 for a derivation) 



\n£m) t 



\{n~£){n+t+l)-\m 2 



\ 



n 



+ m}\[£ - m}\[2£ + 1] 
[n + e + l]U 



x 



(n-t)/2 \_(l+m)/2\ (_-\Xt Q -(n+e+l)t tr>\ \p+t / r>Un+m-2p-2t f r>] \p+t-m 

h „ { o, m) {n ] -[n-i-2t}\\ [2p}\\ [n + m-2p]\[2p-2m}\\ lh 1 ' > 

where 10) is the vacum state, [n]\ = [n] [n — 1] . . . [1], [n}\\ = [n][n — 2] ... [2] or [1], and, as shown 
in r 57 lr 59 l, form a basis for the most symmetric representation [n, 0, 0] of u g (3), corrresponding to the 
u 9 (3) D so g (3) chain. 

28.2. Definitions of so q (3) irreducible tensor operators, tensor and scalar products 

We shall recall now some definitions about g-deformed tensor operators within the framework of the 
algebra so g (3). 

An irreducible tensor operator of rank j with parameter q according to the algebra so g (3) is a set 
of 2j + 1 operators Tj™, satisfying the relations 



[L 0l T^\ =mT$, [L ± ,T$} gm q L ° = ^J[j T m}[j ± m + 1] T$ ±1 , (28.12) 

where, in order to express the adjoint action of the generators Lq,L± of so g (3) on the components of 
the tensor operator Tj™, we use the usual notation of the q- commutator 

[A, B] qa = AB- q a BA. (28. 13) 
By Tj£ we denote the conjugate irreducible q-tensor operator 

f(M) = (-l)i-™,-^ m , (28.14) 

which satisfies the relations 



[T$, L ] = m 7$, q L °[f^, L ± } qm = yj\j ± m][j T m + 1] f^ Tl • (28.15) 
Then the operator 

Vfi = (Tgf)* = {-iy- m q- m {T^ m )\ (28.16) 

where t denotes hermitian conjugation, transforms in the same way (given in eq. (28.12)) as the tensor 

Tjm, i.e. also is an irreducible so q (3) tensor operator of rank j. 

Let and B^ 2 be two irreducible tensor operators. We shall define the tensor and scalar 

oduct of these tensor operators following some of the prescriptions, summarized, for example, in 
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For the tensor product one can introduce the following operator 



^31 A °72 



(93) _ q 3r ,jm Aqi) »(? 2 ) /r,o i y\ 



m\,m,2 



where 93 Cj^ n j 2 OT 2 are ^ ne Clebsch-Gordan coefficients corresponding to the deformation parameter q 3 . 
In general the deformation parameters qi,q% and q$ can be arbitrary. It turns out, however, if one 
imposes the condition that the left hand side of eq. (28.17) transforms as an irreducible g-tensor of 
rank j in a way that the Wigner-Eckart theorem can be applied to eq. (28.17) as a whole, not all of 
the combinations of q\, q2 and q$ are allowed. 

If the tensors A-i 91 and depend on one and the same variable and act on a single vector, which 
depends on the same variable, the mentioned requirement will be satisfied only if q\ = q2 = q and 
g 3 = 1/q, i.e. the operator 



a(q) x »(*' = v 1 /«rtf m A (q) B {q) f28 18) 



(V<z) _ l/qni™ a(i) »(<?) 



mi ,r?i2 



transforms as an irreducible g-tensor of rank j according to the algebra so q (3) . Then, the definition of 
eq. (28.18) is in agreement with the property 

(a',£'\\[A^xB^}f^\\a,£) = 



[2j + 1] E(-i 



J2 



{a.\l'\\Af\\a",f){a",f\\Bf\ 



a,£), 



(28.19) 



which is a g-analogue o f the w ell known classical identity for the matrix elements of the tensor product 
of two tensor operators 



271] 
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In what follows we shall consider also the scalar product of two tensor operators depending on two 
different variables (1) and (2) and acting on different vectors, depending on these different variables. 



In this particular case the scalar product of irreducible q and q 1 - tensor operators Aj (1 
with the same rank, acting on different vectors (1) and (2), is given by means of the following definition 



Ef /q) {2) 



{Af{l) ■ Bf' q \2)Y = (-i)vpy + 1] 



?(!/<?), 



i (<?) i 



1) x B) 



(1/9) 



(2) 



(?) 

00 



E(-<7)"V^U) B?J%(2) . (28.20) 



(!/<?), 



28.3. Construction of so q (3) vector operators and spherical vector operators 

It should be noticed that the simlified angular momentum operators of eq. (28.7) satisfy the com- 
mutation relations of eq. (28.8), but are not tensor operators in the sense of definition of eq. (28.12). 
One can form, however, the angular momentum operators 

- =c 1 n~ L ° T i 
*-±i r i — 1 lj ±y 



'[2] 

4 9) = ij {q[2L } + (q- q- 1 )L.L+} = ± {q[2L ] + (q- q' 1 ) (<# - [L ][L + 1])} , (28.21) 

which are tensors of first rank, i.e. so g (3) vectors. 

One can easily check that the operators Bq, B± and B , B± do not form a spherical vector (see 
for example 00IM21). The situation is the same as in the "standard" theory of angular momentum. 
Indeed, in the "standard" theory, the operators L = L 3 , L± — L\ ± iL 2 , do not form a spherical 
tensor. However, in the "standard" theory one constructs the operators J± = T^C^i ± iL 2 ), Jo = L , 
which are indeed the components of a spherical tensor of first rank according the standard so(3) algebra. 

Here we proceed in an analogous way. As shown in c53, one can define a vector operator Ti\ of the 
form 



T\ = -L {bU 2N+ ~ N ^ - (q - q-^B+iBlfq-^+^+l] . (28.22) 
J\9] I J 



'[2] 

The corresponding expressions for the conjugate operators T m = (— l) m g~ m (Tl fn )"'" are 
f +1 = - (q - q- l )q- 2N + +N ^Bi{B f\ , 



T = q~ 2N++N B , T_i = L — q -™++ N +\ B+ . (28.23) 



One can easily check, that the vector operators and T satisfy the commutation relation 
[T-i, 7+i] g -2 = -q 2N+ \ [f ,T t] = q™ + q~\q 2 - g^T^, 

[T+i, 2li] r > = -g^" 1 + q~\q 2 ~ q~ 2 ){T ] f Q + (q - g^T^-i} , (28.24) 



and 

[To, Tl] = 0, [f +1 , Tl] q i = 0, [T +1 , Tj] = (g 2 - g" 2 )^, 

—2\nr$r 



[T_ 1; T T ] = 0, [T-uT^ = 0, [TcTlJ = (g 2 - fT 3 )!^. (28.25) 



Unlike the operators of eq. (28.2), the commutators [T m ,T n ] and [7^,7^] do not vanish, but are equal 
to 



[7^,7^=0, [TlTl x 



t m\ 



Ug 2 



0, [T +1 ,T ] 



0. 
0. 



= ( g -g" 1 ) (f ) 2 , 



(28.26) 



which is in agreement with the results obtained in . 



28.4- A choice for the physical angular momentum 

It should be noted that the angular momentum operator C$ (considered as a vector according 
so q (3)) can be represented in the form 



r (q) 



(V?) 

XM 



T f x T 

Its "square" differs from C$ and equals to EUSEi 



111 1/qrilM rp] fp 
r i °lm,ln i m i n- 



\ [2] 



(28.27) 
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(£{g)) 2 = . = YX-qY M ^M^u = Ac { 2 q) + ( q q 



[2] 



[2] 



(C : 



(</)\2 



(28.28) 



This difference, however, is not essential. Let us consider, for instance, the expectation values of the 
scalar operator of eq. (28.28). It has the form 



q {lm\£«> ■ C M \£m) q = = [e]q2[£ + l]q2 . 



(28.29) 



In this sense the replacement of 9 L 2 = with ■ is equivalent to the replacement q — > q 2 
and leads to the renormalization of some constant. For this reason we shall accept as the square of the 
physical angular momentum the quantity 9 L 2 = whose eigenvalues are [£][£ + 1]. 



28. 5. A choice for the Hamiltonian 

Now let us consider the scalar operator, constructed in terms of and 7m. We have 

X ( q) = -y/\$[T* x f\§/<> 
We define the Hamiltonian of the three dimensional g-deformed oscillator (3-dim g-HO) as 



-q~ x T\. x T_x + TqT — qTl x T x . 



q H 3dim = Tluj X = hu (-q l T ] +1 T_ x + tJTq - qTl x T x 



(28.30) 



(28.31) 



The motivations for such an ansatz are: 

a) The operator so defined is an so g (3) scalar, i.e. it is simultaneously measurable with the physical 
g-deformed angular momentum square 9 L 2 and the z-projection L . 

b) Only this so g (3) scalar has the property of being the sum of terms conserving the number of 
bosons. (Each term contains an equal number of creation and annihilation operators.) 

c) In the limit q — > 1 the Hamiltonian of eq. (28.31) goes into 



lim q H- 



3dim 



huQ (a) +x a +x + ajao + a)_ x a_ x 



(28.32) 



where [a m ,ajj = S mn , i.e. in this limit the Hamiltonian of eq. (28.31) coincides with the Hamiltonian 
of the 3-dimensional (spherically symmetric) harmonic oscillator up to an additive constant. 

What we have done so far can be summarized in the following way. We have shown that in the en- 
veloping algebra of so g (3) there exists an so g (3) scalar operator, built out of irreducible vector operators 
according to the reduction u g (3)Dso g (3), for the case of the most symmetric irreducible representations 
of u g (3), which are the space of the 3-dim g-HO. In the limit q — > 1 this operator coincides with the 



Hamiltonian of the 3-dim harmonic oscillator. It is therefore reasonable to consider this operator as 
the Hamiltonian of the 3-dim g-HO. (In the next subsections it will be shown that this Hamiltonian is 
related to the Hamiltonian of the Modified Harmonic Oscillator of Nilsson 26^, 270 ,^7i ^72 _ \ 



It should be underlined that the Hamiltonian of eq. (28.31) does not commute with the square 
of the "classical" (or "standard") angular momentum L 2 . This is due to the fact that the g-deformed 
oscillator is really space deformed and the "standard" quantum number / is not a good quantum number 
for this system. On the other hand 9 L 2 commutes with q H^dim and the "quantum angular momentum" 
£ can be used for the classification of the states of the g-oscillator. The projection of the standard 
angular momentum on the z axis, l z , coincides however with the quantum projection, £ z , i.e. it is a 
good quantum number. 

Eq. (28.31) can be cast in a simpler and physically more transparent form. Indeed, making use of 
the third component of the g-deformed angular momentum (considered as an so q (3) vector) 



C 



(?) 



[4]/[2] 



T f x T 



(1/5) 
10 



-Tl^ + iq-q-^fo + T^n, 



we obtain 
Since 



+ qCf 



[2] 



we get upon summation 
and, after some calculations, 



-q 



2JV++27V+1 



[2]ri a r_! + q 2 T$T . 

[N ]q- m +- N ° +2N - 1 , 

qC<8j + [N + L ]q N - L ° + \ 



q H- 



3dim 



HuJoX^ = hu { [N]q 



N+l 



q(q-q 



-ci q) 



The eigenvalues of such a g-deformed Hamiltonian are 

g(g-g- x ) 



q E 3dim = hu [n]q n+L - 



[2] 



-[£}[£ +1} 



[2] 



n, n — 2, 



(28.33) 
(28.34) 
(28.35) 

(28.36) 
(28.37) 



.0 or 1. 



(28.38) 



In the g — ► 1 limit we have lim^! q E^i m = Huqu, which coincides with the classical result. 
One should note that the expression of eq. (28.38) can also be put in the form 



q E^ 



3dim 



HuJq 



n+l 



[2] 



{q e+1 [n-£}+q- e - 1 [n + £}} 



(28.39) 



which perhaps looks more elegant. However, the form of eq. (28.38) has the advantage that the 
vibrational ([n]) and the rotational ([£][£ + 1]) degrees of freedom are clearly separated. 

For small values of the deformation parameter r, where q = e T , one can expand the rhs of eq. (28.38) 
in powers of r obtaining 

q E 3dim = Kujq n — KujqT (£(£ + 1) — n(n + 1)) - hu r 2 (V + 1) - \<n + l)(2n + 1)) + C(r 3 ). (28.40) 

To leading order in r the expression of eq. (28.40) closely resembles the one giving the energy eigenvalues 
of the Modified Harmonic Oscillator, as we shall see in the next subsections. 



?. 6. The Modified Harmonic Oscillator of Nilsson 



Before examining the connection he tweeri the Harniltoniari of the last subsection and the Modified 

3 



Harmonic Oscillator (MHO) of Nilsson 



2(i!) 



27i 17 1 . 272 



a few words about the MHO are in place. 



The basis of the shell model in finite nuclei is the assumption of an independent particle motion 
within a mean field. The average field in which the particles move is correctly described by the Woods- 
Saxon (WS) potential Vws(t). However, the corresponding eigenvalue equation can be solved only 



numerically. It is therefore desirable for many purposes to have an exactly soluble model producing 
approximately the same spectrum as the Woods-Saxon potential. The Modified Harmonic Oscillator 
(MHO) Hamiltonian, first suggested by Nilsson EH and studied in detail in 5Z2clillZa, indeed reproduces 



approximately the energy spectrum of the WS potential and at the same time has a simple analytical 
solution. The Hamiltonian has the form 



V MHO = lhup 2 -hu;k{ f i(L 2 -(L 2 ) N )+2L-S}, P=] j^r, (28.41) 

where L 2 is the square of the angular momentum and L • S is the spin-orbit interaction. The subtraction 
of the average value of (L 2 ), taken over each TV-shell, 

L% - (28.42) 

results in avoiding the shell compression induced by the L 2 term, leaving the "center of gravity" of each 
iV-shell unchanged. The Vmho (without the spin-orbit term L • S) indeed reproduces effectively the 
Woods-Saxon radial potential. A better agreement with the experimental data can be achieved if the 
parameters of the potential are let to vary smoothly from shell to shell. 

It should be noted, however, that from a mathematical point of view the introduction of an In- 
dependent term (and of the corresponding correction (L 2 )^) in the MHO potential is not so innocent, 
as it may look at first sight, because this term depends on the state in which the particle occurs, and 
it this sense the potential of eq. (28.41) is "non-local" and "deformable" (with variable deformation). 
This effect is amplified by the fact that the constant p! = kp is allowed to vary from shell to shell, i.e. 
depends on N. 

The eigenvalues of the MHO (if one neglects the spin-orbit term) are 

E nl = huon - hoop! + 1) - -n(n + 3)) , (28.43) 
where p! is allowed to vary from shell to shell. 



28.7. Connection between the 3 -dimensional q-deformed harmonic oscillator and the Modified Harmonic 
Oscillator of Nilsson 

Comparing eqs (28.40) and (28.43) we see that eq. (28.40) closely resembles, to leading order in r, 
eq. (28.43), thus establishing a connection between the 3-dimensional g-deformed harmonic oscillator 
and the Modified Harmonic Oscillator of Nilsson, with the parameter r of the 3-dim g-HO playing the 
role of the parameter p' = kp of the MHO. (Remember that the spin-orbit term has been omitted in 
both cases.) 

Quantitatively, one obtains a good fit of the energy spectrum produced by the MHO potential of 
eq. (28.41) by choosing for Tioj® and r of the 3-dimensional g-deformed harmonic oscillator the values 
hu = 0.94109 and r = 0.021948 El. (It has been assumed that fhu = 1 for the MHO and the spin-orbit 
term has been neglected.) The value of r is close to the values usually adopted for p! = kp. Indeed 
the MHO fit in the 2 ° 8 Pb region yields 0E1 p' = for iV = 2, p! = 0.0263 for N = 3, and p! = 0.024 
for N > 4, while for k, fixed by the condition that the observed order of sub-shells be reproduced, the 
relevant values are k = 0.08 for N = 2, k = 0.075 for iV = 3 and k = 0.07 for iV > 4. 

It appears surprising at first that Tlujq differs slightly from 1. Indeed the correction term n(n + 1) 
in the 3-dim g-HO spectrum is slightly different from the corresponding piece n(n + 3)/2 in the MHO. 
Some compression of the 3-dim g-HO spectrum is therefore inevitable. From a physical point of view, 
one may say that the mean radius of the deformed oscillator is slightly larger than the radius of the 
classical isotropic HO. 

When comparing the MHO and 3-dim g-HO spectra it should be remembered that the constant p! 
in the MHO potential takes different values for shells with different values of N, while in the 3-dim 
g-HO model the parameters r and Tiojq have the same values for all shells. The comparison shows that 
the g-deformation of the 3-dimensional harmonic oscillator effectively reproduces the non-locality and 
the "deformations" induced in the MHO model through the terms L 2 , (L 2 )^ and the variability of pi 
with the shell number N. 



28.8. 3- dimensional q- deformed harmonic oscillator with q-deformed spin-orbit term 

For a full comparison with the MHO and a more realistic description of the single-particle spectrum, 
one needs to include a g-deformed spin-orbit term ■ S^- 1 ^ 9 ' in the g-deformed harmonic oscillator 
Hamiltonian q H 3 dim- To this purpose one introduces the spin operators S + , So, SL, which are elements of 
another (independent) su 9 (2) algebra. These operators satisfy the commutation relations of eq. (28.8), 
i.e. 

[5-0, S±] = ±S ± , [S+, S-] = [2S ], (28.44) 

and act in the two-dimensional representation space of this algebra. The orthonormalized basis vectors 
of this space will be denoted by \^m s ) q , in analogy with the usual (non-deformed) case. 
One can now define the g-deformed total angular momentum as 1 
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Jo = L + S , J± = L ± q So + S ± q- Lo , (28.45) 

where the g-deformed orbital angular momentum is given by eq. (28.7). The operators of eq. (28.45) 
satisfy the same commutation relations as the operators of the g-deformed orbital angular momentum 
(eq. (28.8)) and spin (eq. (28.44)), and the corresponding expression for the Casimir operator of the 
algebra of eq. (28.45) can be written in the form 

C^j = J-J+ + [J Q ][J + l]. (28.46) 

As in the case of the g-deformed orbital momentum we shall consider q J 2 = C^j as the total angular 
momentum "square". 

The common eig envectors of C^j, Cf L and Cf s can be written in the usual form 

\n{t\)jm) q = q ° jm 1 \n£m e ) q \±m s ) q , (28.47) 

m,i,m s 2 

where q C J e ™ srrls are the g-deformed Clebsch-Gordan coefficients 0000. 

One can define the spin-orbit term as a scalar product, according to the definition of eq. (28.20) 



■ = £ {-q)»ms£lt> = i jcg - C% - C$ - ^^Cgcgj , (28.48) 

where S^ 1 /^ is a vector operator according to the algebra of eq. (28.44), constructed from the g-deformed 
spin operators taking into account the rule of eq. (28.21), but for deformation parameter 1/q. One can 
then suggest the following Hamiltonian for the 3-dim g-deformed harmonic oscillator with spin-orbit 
interaction 

q H = hu [X^ - k[2}£^ . ^(i/?)! = 



-1\2 



»■* I [N]q N+1 - "-^P-C^l - k let] - - - (LZ|3! C « C g ) \ . ( 28 . 49) 



In eq. (28.49) the various factors have been chosen in accordance with the usual convention used in 
the classical (spherically symmetric) shell model with spin-orbit coupling. The eigenvectors of this 
Hamiltonian are given by eq. (28.47), while the corresponding eigenvalues are 

q E nl0 = hco I [n]q n+1 - g(g ~ ^ [I] [I + 1] 

-« (m + 1] - iw+i] - im - {q ~^ 1)2 m+mii]))- (28.50) 



In order to compare the expression for the g-deformed L-S interaction (the term proportional to 
k) with the classical results one needs the expansion of the expectation value of the term given in eq. 



(28.48) in powers of r (g = e T , with r real) 



li(l + Tl(A£ 2 -7) + ...) if j =£+1/2 




{V - -3) + ...l if./ i • 1/2 



(28.51) 

It can be easily seen that the expression of eq. (28.48) introduces some corrections to the classical 
expressions for the spin-orbit interaction, which are proportional to r 2 . These corrections are small for 
light nuclei, where the values of I are small, but they are non-negligible for heavier nuclei, where shells 
with higher values of £ are important. 

One can easily see that the level scheme produced by the 3-dim g-HO withspin-orbit interaction 
reproduces very well the level scheme of the MHO with spin-orbit interaction £13. 

29. Further development of u 9 (3)Dso g (3) models 

The u g (3)Dso g (3) mathematical framework developed in the last section for the case of the 3- 
dimensional g-deformed harmonic oscillator is of more general interest, since one can use it for studying 
-deformed versions of various models having the u(3)Dso(3) symmetry, such as the Elliott model 
and the su(3) limit of the Interacting Boson Model. Since the quadrupole degree of freedom 
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plays a central role in nuclear collectivity, one needs to be able to calculate matrix elements of the 
quadrupole operator in the u ? (3)Dso 9 (3) basis. For the quadrupole operator itself two choices appear: 

a) The quadrupole operator can be a g-deformed second rank tensor operator, with its matrix 
elements being calculated in the g-deformed basis. In this case one assumes that observables are 
described by g-deformed operators and physical states are described by g-deformed wavef unctions. 

b) The quadrupole operator can be the usual (non-deformed) operator, with its matrix elements 
being calculated in the g-deformed basis. In this (less extreme) case one assumes that observables are 
described by the usual physical operators, while the states are described by g-deformed wavef unctions, 
corresponding to highly correlated eigenfunctions of the undeformed Hamiltonian. 

In what follows both cases will be examined. 

29.1. Matrix elements of a q-deformed quadrupole operator in the u q (3)Z)so q (3) basis 

We are going to compute the reduced matrix elements of an so 9 (3) quadrupole operator (second rank 
tensor operator under so 9 (3) ) in the u g (3)Dso ? (3) basis (for the most symmetric u 9 (3) representation 
and for q being real). 

The generalization of the Wigner-Eckart theorem in the case of the algebra so g (3) is 
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riL M' 

(a', L'M'\Ti\a, LM) = (-l) 2j 9 . LM ' jm (a', L'\\Ti\\a, L), (29.1) 

+ 1] 

where \a,LM) are orthonormalized basis vectors of the irreducible representation q D of the algebra 
so g (3) and qC^M^ LaM2 are the Clebsch-Gordan coefficients of the same algebra. 

Our purpose is to calculate the reduced matrix elements of the g-deformed quadrupole operator Q 2 , 
namely the quantities (A, L + 2||<5 2 || A, L) and (A, L||<5 2 || A, L), where the operator Q 2 is 



Q 



2 

M 



A-ki A 2 M — [T" 1 " (8) T] M — Q-iC^^T^Tn. (29.2) 



\ [2] 



m,n 



It is clear that this operator is a second rank tensor under so g (3), according to the definitions given in 
the previous section, while in eq. (29.2) the overall constant factor guarantees the agreement between 
the present expression and the usual classical operator in the limit q — *■ 1. Other reduced matrix 
elements do not occur, since in the most symmetric representation of u g (3) only states with equal parity 
of A and L occur. 



Using the Wigner-Eckart theorem (eq. (29.1)) one obtains for the reduced matrix elements of the 
quadrupole operator Q 2 of eq. (29.2) the following results (see EiZl for a proof) 



(A,L + 2||Q 2 ||A,L> 



[3] [4] 
[2] \ [2] \ 



[A - L] [A + L + 3] [2L + 4] [2L + 2] 
[2L + 3] ' 



(A,L||Q 2 ||A,L> 



[2L] [2L + 1] [2L + 2] 
PTxl [2L-l][2L + 3] 



{g L -5[A-L]+g- L+ 5[A + L + 3]} 



(29.3) 



(29.4) 



which are in agreement with the classical case in the limit q — > 1. Taking into account the Wigner- 
Eckart theorem and the symmetry properties of the g-deformed Clebsch-Gordan coefficients it can also 
be shown that the reduced matrix element given in eq. (29.3) has the following symmetry property 



(A,L + 2||Q 2 ||A,L> = (A,L||Q 2 ||A,L + 2>. 



(29.5) 



For small values of the deformation parameter r (q = e T , with r being real) the following Taylor 
expansions for the reduced matrix elements of eqs (29.3) and (29.4) are obtained 



<A,L + 2||Q 2 ||A,L> 

xh + (x- 



'6(A-L)(A + L + 3)(L + 2)(L + 1) 



(A,L||Q 2 ||A,L> = -(2A + 3) 



T + 



V 2L + 3 

2^ 2 l r2 3 r 65\ 

-A 2 + -L + -L H 

3 2 2 24/ 



r 2 + 0(r 3 )} 



(29.6) 



L(L + 1)(2L + 1) 



1 + 



2A + 3 



{A(A + l)-L(L + l)}r 



+ 



\ (2L- l)(2L + 3) 

^— {(2A + 15)L(L + 1) + (2A + 1)(2A 2 + 2A + 3)}r 2 + 0(r 3 ) 1. 
i + 3) J 



3(2A + 3) —J" (297) 

In the limit r = one obtains the classical expressions for the corresponding reduced matrix elements. 

From a physical point of view of great interest are the E2 transition probabilities (B(E2) factors), 
which in the classical case are expressed by means of the reduced matrix elements of the u(3) quadrupole 
operator. By analogy the B(E2) factors corresponding to the chain u q (3)Dso g (3) can be assumed to be 
of the form 



B(£2;(A,L + 2)-(A,L)) 4 



1 



-|(A,L + 2||Q 2 ||A,L>r 



(29.* 



[2L + 5] 

where Q 2 is the g-deformed quadrupole operator of eq. (29.2). Again in analogy with the classical case, 
the reduced matrix element of eq. (29.4) is related to the deformation of the physical system in the 
state with angular momentum L. 

It should be remembered that the results obtained here are valid only for real values of the deforma- 
tion parameter q. On the other hand the comparison of the predictions of various models based on the 
su g (2) algebra with the experimental data (sections 19-23) shows that one can achieve good agreement 
between theory and experiment in the case of q being a phase factor (q = e"~, with r real). It is therefore 
desirable to repeat the steps followed in this section for the case of q being a phase factor. 



29.2. Matrix elements of the usual quadrupole operator in the u q (3) D so q (3) basis 

In order to be able to calculate matrix elements of the usual quadrupole operator in the g-deformed 
basis of eq. (28.11) we should first express the basis vectors in terms of usual (undeformed) bosons a\, 
d, which are connected to the g-deformed bosons b\, hi through HS^l 



*t = V^U fc = a *vW' * = +' '-' ( 29 - 9 ) 



where Ni is the number operator of the undeformed bosons. In this way the operators B\, B t of eq. 
(28.2) can be written in the form 



B = q 



q 1 2 di 



'[2Ni\ 



[N t 



o 



No 
B\ 



-a q 



z = +,-. 



(29.10) 



(29.11) 



N> 1 V Ni 

As a next step we can replace B\ in eq. (28.11) with their expressions from eq. (29.11). Then, using 
the following identity, which holds in the Fock space, 



{Bjy (Bpy (Bly = k^-i^d (4r (4)!(aLr |0) 

[2s]!! J\y\\ J\2z\U Vx\ \fz\ ' 



(29.12) 



we obtain (in a slightly different notation) EHEE 



A 

L,M 



■\{(2X+M)(M- 



(A-L)/2 l(L+M)/2} 



1)+L[L+l)] N^ M E " E 



t„(A-i)x-(L+|)t 



t=0 x=max(0,Af) 



[2t]!![A-L-2t]!! 



x 



[2x + 2t]!![A + M — 2x — 2t]\[2x + 2t- 2M]\\ ( a \) x+t 



a, 



]\X+M-2x-2t 

oJ 



[a 



t \x+t~M 



[2x\\\[L + M- 2x]![2x - 2M]H 



= |0>, 



where the normalization factor N^ M is 



^(x + t)\ + M - 2x - 2t)\ yj(x + t- M)\ 

(29.13) 



N 



(?) 



XLM 



[A - L]!![L + M]![L - M]![2L + 1] 
[A + L + l]!! 



(29.14) 



The expression of eq. (29.13) can be used for the calculation of the matrix elements of "standard" 
(or classical) tensor operators (which are built up from standard boson operators) between the states 
with given g-deformed angular momentum. 

Let us consider now the classical quadrupole operator Q$. Its zero component is of the form 



Qo = V&[a^ <S> a] 2 o = a+a_ + 2aJao + a!_a + = — a+a + + 2aJao — ala_ = 3N — N, 



(29.15) 



where a^, a m 



m 



-, 0, — ) are classical boson operators and a m = (— l) m a_ 



The explicit expression for the matrix elements of Q between the g-deformed states of eq. (29.13) 
then becomes 



A 

1M' 



Q 



A 

L,M 



6* 



>M,M' \ 
Q <1 



A 

J,M 



Q 



A 

L,M 



(29.16) 



and 



A 



A 

L,M 



(2A + 3M)5 L ,j - 6g -3(^)(w-i)-i{L(i+i)+J(J+i)}^) MiV W M 

[2x]!![A + M-2x]![2x-2M]!! 



(A-L)/2 (A-J)/2 f _.s t+r / X +L+l)t-(X+J+l)r 

x V V 1 > q V 

h> r=o l 2t } ]] l X ~ L ~ 2t]!![2r]!![A - J - 2r]!! V [2z - 2t]\\[L + M - 2x + 2t]\[2x - 2t - 2M]!! 



where: 



[2x - 2r]!![J + M - 2x + 2r]![2x - 2r - 2M]H ' 
max{t, r, M + t, M + r} < x < min{ [(L + M) /2J + t, [(J + M) /2J + r}. 



(29.17) 
(29.18) 



Clearly the recipe described above can be applied for the calculation of matrix elements between 
the g-deformed angular momentum states of any classical operator. The disadvanage of this method is 
that one must perform special calculations for any specific operator. This can be avoided by the use of 
the transformation matrix between the "standard" and g-deformed states, which will be described in 
the next subsection. 



29.3. Transformation between so q (3) and so(3) basis states 

For fixed A one can expand the g-deformed basis states of eq. (29.13) in terms of classical ones 
(with q — 1), since both sets form orthonormal bases in the eigenspace 7i\ of the number operator 
N = N + + iVo + iV_, corresponding to the eigenvalue A. In particular, dim7^A = ^(A + 1)(A + 2) and 
the operator L has the same form L = N + — N_ in both classical and g-deformed cases. In this way 



(29.19) 



where J = A, A — 2, ... , \M\ or \M\ + 1 and the subscript c denotes the classical (undeformed) basis 
states. Taking the scalar product between g-deformed basis states of eq. (29.13) and classical basis 
states, for the transformation matrix in eq. (29.19) one obtains 



OM,M' y 4 iv AJM ly XLM 










(A-L)/2 (A-J)/2 

x E E 



-iy+r q-(\+L+l)t 



E 



q^-^ x y/[2x]\\[X + M - 2x}\[2x - 2M]H 



t=0 



where 



r=0 



x 



[2t]!![A — L — 2t]!!(2r)!!(A -J - 2r)!! V [2x - 2t]\\[L + M - 2x + 2t]\[2x -2t- 2M}\\ 

(29.20) 

(29.21) 



(2x)!!(A + M- 2x)\(2x - 2M)\\ 
(2x - 2r)!!(J + M - 2x + 2r)!(2a; - 2r - 2M)\\ ' 



max{t, r,M + t,M + r} < x < min{ [(L + M) /2J + t, [(J + M) /2J + r}. 
It should be noted that the scalar product given in eq. (29.20) is real and has the properties 
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M,M' 







(29.22) 



which have been used in eq. (29.19). Though looking very involved, eq. (29.20) is quite handy for 
computational purposes. Indeed, only few terms are contained in the sum. The same equation shows 
that the transformation matrix is diagonal in the third projections of J and L, and is nonvanishing even 
for high values of the difference A (J) = \L — J\. However, the values are peaked around the classical 
value A (J) =0, if the parameter of deformation q tends to unity. 

Let us consider now a classical tensor operator of rank j according to the algebra so(3), which 
conserves the number of particles. Using the Wigner-Eckart theorem in the classical case we can express 
the matrix elements of between the classical so(3) basis states in the form 




A 



(c) 
jm 




fi.JM' 

= ^ ^2.7 LM,jm n 7 || A (c) 



V2J + 1 



{XJ\\A \\\L), 



(29.23) 



where Cj^jm are ^ e classical Clebsch-Gordan coefficients of the so(3) algebra, and we assume that 

the classical reduced matrix elements of A^ in eq. (29.23) are known. In this way, from eq. (29.19) 
follows that 

Ac) 




A 




A 


R,M 




L,M 



= EE 

R P,Mp 

and using eq. (29.23) we have the expansion 




Ac) 

jm 






(29.24) 
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EE( 

R P 



(XPWA^WXR) 



fiPM' 



V2P+1 




A 
L. 



M 



(29.25) 



In particular, eq. (29.25) shows that the matrix elements of the classical tensor operator between 
g-deformed so g (3) states are nonvanishing if M' = M + m, as in the classical basis. It is also worth to 
remark that the summation expression of eq. (29.25) gives an universal way for (numerical) computation 
of matrix elements of different types of classical operators between the deformed basis states of eq. 
(29.13). For instance, it is straightforward to show that the numerical values of the quadrupole operator 
calculated by eq. (29.17) of the previous subsection and the method used in this subsection (involving 
transformation matrices) coincide. However, the method descibed in this subsection is universal and 
can be applied to any classical operator. 

Examples of transformation brackets can be found in 
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30. The question of complete breaking of symmetries and some applications 

In the cases examined so far the g-deformed symmetries considered were close to their classical 
counterparts, to which they reduce for r = (g = 1), since the values of r were relatively small. 
One can then argue that results similar to the ones provided by the quantum symmetries can also 
be obtained from the usual Lie symmetries through the addition of suitable perturbations. What can 
be very useful is to start with one limiting symmetry and, through large deformations, reach another 
limiting symmetry. We shall refer to this as the complete breaking of the symmetry. In this way one 
could hope to "bridge" different Lie symmetries through the use of q- deformations, providing in addition 
new symmetries in the regions intermediate between the existing Lie ones. 

The question of complete breaking of the symmetry in the framework of the toy IBM of sec. 27 has 
been studied by Cseh E3, Gupta El, and Del Sol Mesa et al. EE. 



Cseh till started with the su g (2) (vibrational) limit (in a form different from the one used in sec. 
27) and tried to reach the so g (3) (rotational) limit. He noted that for q being a phase factor this is not 
possible, while for real q some rotational features are obtained, but without all of the requirements for 
rotational behaviour being satisfied simultaneously. 

Gupta tla started with the su g (2) limit, in the form given in sec. 27. He noted that for q being a 
phase factor a recovery of the su(3) symmetry occurs (see also the next paragraph), while for real q, 
and even better for q complex (q = e s with s = a + ib), the so g (3) limit is indeed reached. 

Del Sol Mesa et al. c£-3 considered the o g (3) limit of the model of sec. 27, since it corresponds to the 
symmetry of a g-deformed version of the spherical Nilsson Hamiltonian with spin-orbit coupling term. 
They found that for q being a phase factor (q = e lT ) and for r obtaining values in the region 0.5 < r < 2 
the u(3) symmetry, which is broken in the initial model because of the presence of the spin-orbit term, 
is recovered. This offers a way of recovering the u(3) symmetry a lterna tive to the one developed for the 



spherical Nilsson model P 79 IP 80 I and the deformed Nilsson model EBt^l through the use of appropriate 
unitary operators. 



Complex deformations have also been used in in the framework of a deformed u(2) model 



possessing the u(2)Du(l) and u(2)Do(2) chains. Again it has been found that complex deformations 
can bridge the two limiting symmetries. 

Possible complete breaking of the symmetry has also been studied in the framework of the g-deformed 
version of the full Interacting Boson Model (see sec. 31). 

A problem associated with complex deformations as the ones considered above is that the energy 
eigenvalues become complex as w ell. A way to avoid this problem has been introduced recently by 



Jannussis and collaborators 



Finally, the o q (3) limit of the model of sec. 27 has been used for describing the 16 + a cluster 
states in 20 Ne iH. It turns out that an improved description of the energy spectrum and the a-particle 
spectroscopic factors occurs for q = e 
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31. g-deformation of the Interacting Boson Model (IBM) 



The Interacting Boson Model l 179 H 215 lr 16 l (see li£§E^j for recent overviews) is the most popular algebraic 
model of nuclear structure. It describes the collective properties of medium-mass and heavy nuclei 
away from closed shells in terms of bosons, which correspond to correlated valence fermion pairs. In its 
simplest form, called IBM-1, only s ( J — 0) and d (J = 2) bosons are used. The overall symmetry of 
the model is u(6), possessing three limiting symmetries: the u(5) (vibrational) limit, corresponding to 
the chain of subalgebras 

u(6) D u(5) D o(5) D o(3), (31.1) 
the su(3) (rotational) limit, characterized by the chain 



u(6) D su(3) D o(3), 
and the o(6) (7-unstable) limit, for which the relevant chain is 

u(6) D 0(6) D o(5) D o(3). 



(31.2) 



(31.3) 



If one of these dynamical symmetries is present, the Hamiltonian can be written in terms of the 
Casimir operators of the algebras appearing in the relevant chain. Thus the Hamiltonian can be ana- 
lytically diagonalized in the corresponding basis. This is a great advantage of IBM and its numerous 
generalizations: they provide us with a large number of exactly soluble models, the predictions of which 
can be directly compared to experiment, without any need for lengthy numerical calculations. 

From what we have already seen in sec. 27, it is worth examining if a g-deformed version of the 
IBM has any advantages in comparison to the standard version. In order to accomplish this, one has 
to construct the g-analogues of the three chains mentioned above. The difficulties associated with the 
su(3) chain have already been discussed in subsec. 27.3. In what follows we are going to focus attention 
on the 0(6) chain, for which the relevant construction has been carried out Hll. The technique used is 
based on the notion of complementary subalgebras, which is explained in detail in EH, while here only 
final results will be reported. We only mention here that the notion of complementary subalgebras was 
introduced by Moshinsky and Quesne EH-SEEH. Two subalgebras Ai and A2 of a larger algebra A are 
complementary within a definite irrep of A, if there is an one-to-one correspondence between all the 
irreps of Ai and A2 contained in this irrep of A. 



In the 0(6) limit of IBM the Hamiltonian is 

H = E Q + /3C 2 (o(5)) + 7 C 2 (o(3)) + »A(o(6)). 



(31.4) 



The eigenvalues of the energy in the relevant basis have already been given in subsec. 26.1. 

Using the notion of complementarity it turns out that, instead of using the 0(6) chain mentioned 
above, it suffices to study the chain 



su 



sd 1 



1,1) 



so(6) D su d (l, 1) ® so(5) D so(3), 



(31.5) 



where su sd (l,l) is the algebra closed by the pair operators formed out of the s and d bosons, while 
su d (l,l) is the algebra closed by the pair operators fo rme r! out of d bosons alone. (Details on the basis 
for symmetric irreps of su(l,l)®o(5) can be found in E23.) The irreps of su sd (l,l) are characterized by 
the same quantum numbers as the irreps of 0(6) in the 0(6) chain of the IBM, while the irreps of su d (l,l) 
are characterized by the same quantum numbers as the irreps of o(5) in the 0(6) limit of IBM. Therefore 
in the Hamiltonian one can use the Casimir operators of the su sd (l,l), su d (l,l) and su(2) subalgebras 
(the deformed versions of which are well known, as seen in sees 14-16) instead of the Casimir operators 
of 0(6), o(5), o(3) respectively. Keeping the same notation as in eq. (26.3) the final result reads 



E(N, a, r, p A , J, Mj) = E + P8 
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(31.6) 



where the free parameters have been chosen so that the present equation reduces to its classical coun- 
terpart for q — > 1. For the ground state band then the analog of eq. (26.4) is 



E(J) = E + {3'8[J} ql/ 4J + 6] gl/ 4 + -y2[J] q [J + 1], + t/8[iV] 9 i /a [N 



4] 9 i/2, 



(31.7) 



where the identities 



i<l 1/2 + q- 1/2 ) 



X 



W ? v 4 (g 1/2 + ?- 1/2 )- 1 (? 1/4 + g- 1/4 )- 1 , 



(31.8) 



(31.9) 



have been used and /3', rf are related to (3, rj and q in an obvious way. We remark that the Casimir 
operator of sUg(l,l), which is complementary to o(5) in the undeformed case, leads to a term of the 
form [J] g '[« Z + 6b / with q' = g 1 / 4 . 

In refs | 293 lr 94 | the question has been studied if large values of the deformation parameter can lead 
us from the o(6) limit to the su(3) (rotational) or u(5) (vibrational) limits, so that complete breaking 
of the symmetry, in the sense of sec. 30, could be obtained. It turns out that for q real the spectrum 
of the ground state band goes towards the rotational limit, while q being a phase factor leads towards 
the vibrational limit. Many more detailed studies, of both spectra and electromagnetic transition 
probabilities, are required before such a claim can be made. 

A different method fox constructing the g-deformed versions of the u(5) and o(6) limits of the IBM 
has been used by Pan The method isbased on the use of g-deforming functionals (see sees 10, 
14). The same method has been used in for studying the g-deformed version of the su(3)Dso(3) 
decomposition. The final result for the energy eigenvalues in the o(6) case is similar to the one reported 
in eq. (31.6), the main difference being that different deformation parameters are allowed in each of 
the three deformed terms in the rhs. Some comparisons of the model predictions for spectra and B(E2) 
values to the experimental data have been performed 

It is clear that several deformed versions of the IBM can be constructed, providing us with a large 
number of exactly soluble models. In order to demonstrate their usefulness, one has to show that by 
deforming the model one gets some advantages over the classical (non-deformed) version. One way to 
achieve this is the use of parameter-independent tests based on systematics of the data, like the ones 
used in r 96 ^ 97 !! 298 ! for the usual IBM. It is also desirable for the deformation parameter to be associated 
with some physical quantity, as in the case of the su g (2) model. Much work is still required in these 
directions. Some mathematical results which can be useful in these efforts are reported below: 

i) Casimir operators for su,j(n) have been given ini 9 !, while the quadratic Casimir of so g (5) can be 
found in E22 . 



ii) Raising and lowering operators for u g (n) have been given by Quesne E2il. 

iii) Irreps of u 9 (m+n) in the u 9 (m)©u 9 (n) basis have been constructed in E23, while generalized 
g-bosonic operators acting in a tensor product of m Fock spaces have been constructed as double 
irreducible tensors with respect to u g (m)©u 5 (n) in c2i 



32. Deformed versions of other collective models 



The Moszkowski model El3 is a schematic two-level model which provides a description of the phase 
transition from the vibrational regime to the rotational one. A g-deformed version of the model has been 
developed and the RPA modes in it have been discussed IM In addition the validity of the coherent 
states variational method has been checked in the framework of this model and the time evolution of 
the system has been studied by means of a time dependent variational principle 0. Furthermore, the 
g-deformed Moszkowski model with cranking has been studied in thejuean field approximation and the 
relation between g-deformation and temperature has been discussed It should be noticed here that 
quantum algebraic techniques have also been found usefu l in describing thermal effects in the framework 
of the g-deformed Thouless model for superconductivity 

The Lipkin-Meshkov-Glick (LMG) model ^ is an exactly solubl e sch e m a tic s hell model, g-deformed 

and of the 3-level LMG 



versions of the 2-level LMG model (in 
model (in terms of an su g (3) algebra) 



2rms of an su g (2) algebra) E13't 
^13 have been developed. 



33. Fermion pairs as deformed bosons: approximate mapping 



We have seen so far that several quantum algebraic phenomenological models have been proposed 
for the description of nuclear collective properties. These models make use of g-deformed bosons, which 
satisfy commutation relations differing from the standard boson commutation relations, to which they 
reduce in the limit q — > 1. 

On the other hand, it is known that vibrational nuclear spectra, which are described in the simplest 
way by a pairing Hamiltonian, show anharmonicities (see also sec. 26), described, for example, by the 
Anharmonic Vibrator Model (AVM) 111 

E(J) =aJ + bJ(J -2). (33.1) 
r ii,Eii,5i3^2]|J322| 



In the framework of the single-j shell model Eia'Lij'Eijf W I, which can be extended to several non 



degenerate j-shells rf f i these anharmonicities are related to the fact that correlated fermion pairs 



satisfy commutation relations which resemble boson commutation relations but in addition include 
corrections due to the presence of the Pauli principle. This fact has been the cause for the devel opm ent 
of boson mapping techniques (see the recent reviews by Klein and Marshalek e3 and Hecht E2a and 
references therein), by which the description of systems of fermions in terms of bosons is achieved. In 
recent years boson mappings have attracted additional attention in nuclear physics as a necessary tool 
in providing a theoretical justification for the success of the phenomenological Interacting Boson Model 
Ei! and its various extentions, in which low lying collective states of medium and heavy mass nuclei are 
described in terms of bosons. 

From the above observations it is clear that both q-bosons and correlated fermion pairs satisfy 
commutation relations which resemble the standard boson commutation relations but they deviate 
from them, due to the g-deformation in the former case and to the Pauli principle in the latter. A 
question is thus created: Are g-bosons suitable for the approximate description of correlated fermion 
pairs? In particular, is it possible to construct a boson mapping in which correlated fermion pairs are 
mapped onto g-bosons, in a way that the g-boson operators approximately satisfy the same commutation 
relations as the correlated fermion pair operators? In this section we show for the simple case of su(2) 
that such a mapping is indeed possible. 

33.1. The single-j shell model 

Let us consider the single-j shell model ^.^.^^.ill. One can define fermion pair and multipole 
operators as 

A \u = -4 E (jmjm'\JM)a] m a] m ,, (33.2) 

v ^ mm' 

B JM = J— £ (jmj - m'\JM)(-iy- m 'a] m a jm ,, (33.3) 

V ZJ + 1 mm' 

with the following definitions 

Ajm = [A ] jm }\ B\ M =[B JM t (33.4) 

In the above aj m (a>jm) are fermion creation (annihilation) operators and (jmjm'\JM) are the usual 
Clebsch-Gordan coefficients. 

The pair and multipole operators given above satisfy the following commutation relations: 

[AjMi A\i M i] = 0, (33.5) 

[Ajm,A\, m ,\ = 5jj,5mm> - 2 £(-l) 2 ^ +M ^(2J + 1)(2J' + 1)(2J" + 1) 

j" 

(J - MJ'M'\J"M' - M) | J J J . j^M'-M, (33.6) 

[Bjm> Aj J'M'} = E 2v / 2J / + l(-l) 2j '- A/ (J - MJ'M'\J"M' - M) 
j" 



J J ' J "\A^ 1 + f33 7 ) 

j j j j /± .J",M'-M 2 ' ^ > 

[Bj M ,Bj /M/ ] = E(-1) 2J - J "[1 " (-l) J+J ' +J "]v / 2J^TT 
j" 

{JMJ'M'\J"M + M')| J J J |fi,/», M +M', (33.8) 

in which the curly brackets are the usual 6-j symbols. These are the commutation relations of the 
so(2(2j + l)) algebra. 

33.2. Fermion pairs of zero angular momentum 

In the present subsection we will restrict ourselves to fermion pairs coupled to angular momentum 
zero. The relevant commutation relations take the form 

[4>,4l = l-7p [^,41=4, [^,A ] = -A , (33.9) 
where Np is the number of fermions, 2Q = 2 j + 1 is the size of the shell, and 

B = N F /V2tt. (33.10) 

With the identifications 

J+ = VUAl J_ = VnA , Jo = NF ~ n , (33.11) 
eqs (33.9) take the form of the usual su(2) commutation relations 

[J+,J_]=2J , [J ,J + ] = J+, [J ,J-] = -J-. (33.12) 
An exact boson mapping of the su(2) algebra is given in 
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4 = aSyi-^, A = ^1-^, N F = 2n , (33.13) 

where a\ (ao) are boson creation (annihilation) operators carrying angular momentum zero and no is 
the number of these bosons. 

The simplest pairing Hamiltonian one can consider has the form 

H = -GVLA\A Q . (33.14) 

The Casimir operator of su(2) can be written as 

r ,t . -. ^ ( n f\ 2 

{Al,A } + - =17 + 1 > ( 33 - 15 ) 



2 V 

while the pairing energy takes the form 

E N F Np N F 

^GQ)~^r~^n + 2n- (33 - 16) 

Our aim is to check if there is a boson mapping for the operators A\, A q and iV^ in terms of 
g-deformed bosons, having the following properties: 

i) The mapping is simpler than the one of eq. (33.13), i.e. to each fermion pair operator Aq, A 
corresponds a bare g-boson operator and not a boson operator accompanied by a square root (the Pauli 
reduction factor). 

ii) The commutation relations (33.9) are satisfied up to a certain order. 

ii) The pairing energies of eq. (33.16) are reproduced up to the same order. 



33.3. Mapping using the q- deformed oscillator 

In the case of the g-deformed harmonic oscillator (sec. 10), the commutation relation 

[a, a)} = [N + 1] - [N] (33.17) 

for q being a phase can be written as 

[a, « f ] = ° • (33.18) 

In physical situations r is expected to be small (i.e. of the order of 0.01). Therefore in eq. (33.18) 
one can take Taylor expansions of the functions appearing there and thus find an expansion of the form 

r 2 r 4 

[a, a f ] = 1 - — (N 2 + N) + — (iV 4 + 2iV 3 - N) - . . . . (33.19) 

We remark that the first order corrections contain not only a term proportional to N, but in addition 
a term proportional to iV 2 , which is larger than N. Thus one cannot make the simple mapping 

A -»• a, a\ N F -> 2N, (33.20) 

because then one cannot get the first of the commutation relations (33.9) correctly up to the first order 
of the corrections. The same problem appears in the case that q is real as well. In addition, by making 
the simple mapping of eq. (33.20) the pairing Hamiltonian can be written as 

a^a=[N}. (33.21) 



-GQ 



In the case of small r, one can take Taylor expansions of the functions appearing in the definition of 
the g-numbers (eq. (2.2) or (2.3)) and thus obtain the following expansion 

r 2 t 4 r 6 

\N] = N±—(N-N 3 ) + (7N — 10iV 3 + 3iV 5 ) ± (SIN - 49iV 3 + 21iV 5 - 3iV 7 ) + . . . , (33.22) 

1 J 6 V ; 360 V ; 15120 V ; V ; 

where the upper (lower) sign corresponds to q being a phase factor (real). We remark that while the 
first order corrections in eq. (33.16) are proportional to Np and N F , here the first order corrections 
are proportional to iV and N 3 . Thus neither the pairing energies can be reproduced correctly by this 
mapping. 

33.4- Mapping using the Q -deformed oscillator 

In the case of the Q-oscillator of sec. 11, however, the commutation relation among the bosons is 

[b,tf] = Q N . (33.23) 

Defining Q = e T this can be written as 

T 2 N 2 T 3 N 3 

[h, &t] = 1 + TN + -— + — — + . . . . (33.24) 

We remark that the first order correction is proportional to N. Thus, by making the boson mapping 

A\ -> b\ A -> 6, N F -> 2N, (33.25) 

one can satisfy the first commutation relation of eq. (33.9) up to the first order of the corrections by 
determining T = —2/Q. 



We should now check if the pairing energies (eq. (33.16)) can be found correctly up to the same 
order of approximation when this mapping is employed. The pairing Hamiltonian in this case takes the 
form 

JL = &t 6 = [N]q . (33.26) 

Defining Q = e T it is instructive to construct the expansion of the Q-number of eq. (6.1) in powers 
of T. Assuming that T is small and taking Taylor expansions in eq. (6.1) one finally has 

[N]q = N + ^(N 2 - N) + ^-(2iV 3 - 3iV 2 + 1) + ^(N 4 - 2N 3 + N 2 ) + . . . (33.27) 

Using the value of the deformation parameter T = —2/Q, determined above from the requirement that 
the commutation relations are satisfied up to first order corrections, the pairing energies become 

E AT N 2 -N 2N 3 -3N 2 + 1 N 4 - 2N 3 + N 2 

N + — — — — + .... (33.28) 



-GQ Q 3fi 2 3ft 3 

The first two terms in the rhs of eq. (33.28), which correspond to the leading term plus the first 
order corrections, are exactly equal to the pairing energies of eq. (33.16), since N F — > 2N. We 
therefore conclude that through the boson mapping of eq. (33.25) one can both satisfy the fermion pair 
commutation relations of eq. (33.9) and reproduce the pairing energies of eq. (33.16) up to the first 
order corrections. 

The following comments are also in place: 

i) By studying the spectra of the two versions of the g-deformed harmonic oscillator, given in eqs. 
(10.10) and (11.10), one can easily draw the following conclusions: when compared to the usual oscillator 
spectrum, which is equidistant, the spectrum of the g-oscillator is getting shrunk for q being a phase, 
while the spectrum of the Q-oscillator gets shrunk when T < 0. In a similar way, the spectrum of the 
(/-oscillator gets expanded for q real, while the spectrum of the Q-oscillator gets expanded for T > 0. In 
physical situations (sees 19-23) it has been found that the physically interesting results are gotten with 
q being a phase. Thus in the case of the Q-oscillator it is the case T < the one which corresponds 
to the physically interesting case. As we have already seen, it is exactly for T = —2/Q < that the 
present mapping gives the fermion pair results. 

ii) It should be recalled that the pairing model under discussion is studied under the assumptions 
that the degeneracy of the shell is large (Q » 1), that the number of particles is large (N >> 1), and 
that one stays away from the center of the shell (Q — N = 0(N)). The accuracy of the present mapping 



in reproducing the pairing energies has been checked in where results for Q — 11 (the size of the 
nuclear fp major shell), Q = 16 (the size of the nuclear sdg major shell) and Q = 22 (the size of the 
nuclear pfh major shell) are reported, along with results for the case Q = 50 (as an example of a large 
shell). In all cases good agreement between the classical pairing model results and the Q-Hamiltonian 
of eq. (33.26) is obtained up to the point at which about 1/4 of the shell is filled. The deviations 
observed near the middle of the shell are expected, since there the expansion used breaks down. 

We have thus shown that an approximate mapping of the fermion pairs coupled to angular mo- 
mentum zero in a single-j shell onto suitably defined g-bosons (the Q-bosons) is possible. The su(2) 
commutation relations are satisfied up to the first order corrections, while at the same time the eigenval- 
ues of a simple pairing Hamiltonian are correctly reproduced up to the same order. The small parameter 
of the expansion, which is T (where Q = e T ), turns out to be negative and inversely proportional to 
the size of the shell. 

The present results are an indication that suitably defined g-bosons could be used for approximately 
describing systems of correlated fermions under certain conditions in a simplified way. The construction 
of g-bosons which would exactly satisfy the fermion pair su(2) commutation relations will be undertaken 
in the following section. 



34. Fermion pairs as deformed bosons: exact mapping 



From the contents of the previous section, the following question is created: Is it possible to construct 
a generalized deformed oscillator (as in sec. 12) using deformed bosons in such a way that the spectrum 



of the oscillator will exactly correspond to the pairing energy in the single-j shell model, while the 
commutation relations of the deformed bosons will exactly correspond to the commutation relations of 
the correlated fermion pairs in the single-j shell under discussion? In this section we show that such an 
oscillator can indeed be constructed EHa by using the method of sec. 12. 



34-1- An appropriate generalized deformed oscillator 

We apply the procedure of sec. 12 in the case of the pairing in a single-j shell mentioned before. 
The boson number is half the fermion number, i.e. A" = Np/2. Then eq. (33.16) can be written as 

E N 2 N 

^ N - ^ + (34.1) 



-Gn n q 

One can then use a generalized deformed oscillator with structure function 

N 2 N 

F(N)=<Ja = N- — + -. (34.2) 

In addition one has 

F(N + 1) = aa < = N+1 -(Ji+}l + E±l. (34.3) 

What we have constructed is a boson mapping for the operators A , Aq, N f : 

A ^a, Al^a\ N F ^2N. (34.4) 

From eq. (34.2) it is clear that this mapping gives the correct pairing energy. In addition one has 

2N 

[a, at] = F(N + 1) - F(N) = 1 " -jp (34.5) 

in agreement to the first commutation relation of eq. (33.9). Thus the correct commutation relations 
are also obeyed. (The last two commutation relations of eq. (33.9) are satisfied because of (12.1).) 

As we have already seen in the previous section, an exact hermitian boson mapping for the su(2) 
algebra is known to have the form of eq. (33.13). In this mapping the Pauli principle effects are carried 
by the square roots accompanying the ordinary boson operators, while in the mapping of eq. (34.4) the 
Pauli principle effects are "built in" the deformed bosons. 

The generalized oscillator obtained here has energy spectrum 

E N = ±(F(N) + F(N + 1)) = N + l - - ^, (34.6) 
which is the spectrum of an anharmonic oscillator. 

34-2. Related potentials 

The classical potential giving the same spectrum, up to first order perturbation theory, can be easily 
determined (see also subsecs 13.1, 13.2). The potential 

V(x) = kx 2 + Xx\ (34.7) 

is known to give in first order perturbation theory the spectrum 

E n = n(2n + 1) + \(<6n 2 + 6n + 3) = (2k + 6A)(n + -) + 6An 2 . (34.8) 

Comparing eqs. (34.6) and (34.8) one finds 

' t = 5 (1 + S»' A = -6?f < 34 ' 9 > 

Then the classical potential giving the same spectrum, up to first order perturbation theory, as the 
generalized oscillator determined here, is 

v v = 1 -P + b*-wf- (34 - 10) 



It is therefore demonstrated that the Pauli principle effects in a single-j shell with pairing interaction 
are equivalent to an x 4 anharmonicity. 

The generalization of the results obtained in this section for the pairing hamiltonian to any anhar- 
monic oscillator is straightforward. For example, the potential 

V(x) = kx 2 + Ax 4 + /ix 6 + £x 8 , (34.11) 

is known to give up to first order perturbation theory the spectrum 

E n = /t(2n+l)+A(6n 2 +6n+3)+ / u(20n 3 +30n 2 +40n+15)+^(70n 4 +140n 3 +350n 2 +280n+105), (34.12) 
which can be rewritten in the form 

E n = (n + (n + l))(/t + 5 / u) + (n 2 + (n + l) 2 )(3A + 70O + (ri 3 + (n+l) 3 )(10 / u) + (n 4 + (n + l) 4 )(35O. (34.13) 
Taking into account eq. (12.11), from eq. (34.13) one gets 

E^Q. = (« + 5^ + (3A + 70£)n 2 + {I0fi)n 3 + (350^ 4 . (34.14) 

For \i — £ — and n, A given from eq. (34.9), the results for the pairing problem are regained. 

It is worth mentioning at this point that the energy spectrum of the generalized oscillator corre- 
sponding to the pairing correlations (eq. (34.6)) can be rewritten as 

E N= *(4 + ^±( N+ l)l( N+ l) 2 ). (34.15) 



Q \ 8 2 V 2/ 2 V 2 

On the other hand, it is known that for the modified Poschl-Teller potential (see also subsec. 13.2) 

V(x) = Dtanh\x/R), (34.16) 

the energy spectrum is given by 
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E "-^p{ j s + y smDR2 ' h2+i o^K^+i) 2 )- (34 - i7 » 

It is thus clear that the energy spectrum of the generalized oscillator studied here can be obtained from 
the modified Poschl-Teller potential for special values of the potential depth D. 

It is also worth remarking that the "structure function" F(N) of the generalized oscillator obtained 
here (eq. (34.2)) can be written as 

N 

F(N) = —(n + l-N), (34.18) 

which is similar to the one of the para-fermionic oscillator of Ohnuki and Kamefuchi (see also sees 
12,18). 

In summary, we have constructed a generalized deformed oscillator which satisfies the same commu- 
tation relations as fermion pair and multipole operators of zero angular momentum in a single-j shell, 
and, in addition, reproduces the pairing energy of this shell exactly. We have thus demonstrated that 
an exact hermitian boson mapping of a system of angular-momentum-zero fermion pairs in terms of 
bare deformed bosons can be constructed, while in the usual case the ordinary bosons are accompanied 
by square roots due to the Pauli principle effects. The oscillator corresponding to the pairing problem 
has a spectrum which can be reproduced up to first order perturbation theory by a harmonic oscillator 
with an x 4 anharmonicity. The construction of a generalized deformed oscillator corresponding to any 
anharmonic oscillator has also been achieved. 

The results obtained in this section indicate that deformed bosons might be a convenient tool for 
describing systems of fermion pairs under certain conditions. The generalisation of the results obtained 



here to fermion pairs of nonzero angular momentum, which will allow for a fuller treatment of the 
single-j shell in terms of deformed bosons, is a very interesting problem. 



35. The seniority scheme 

In the previous two sections we have seen how correlated fermion pairs of zero angular momentum 
can be mapped onto deformed bosons. It is however known that pairs of non-zero angular momentum 
play an important role in the formation of nuclear properties. In the present section a first step in the 
direction of describing the J ^ pairs in terms of deformed bosons is taken. 



35.1. Uncovering a dynamical symmetry 



In the usual formulation of the theory of pairing in a single-j shell fermion pairs of angular 
momentum J = are created by the pair creation operators 

V" m>0 

where Oj m are fermion creation operators and 2fl = 2j + 1 is the degeneracy of the shell. In addition, 
pairs of nonzero angular momentum are created by the Q — 1 operators 

B\ = £ (-l) j+m Umj - m\J0)a] m a]_ m , (35.2) 

m>0 



where (jmj — m\ JO) are the usual Clebsch Gordan coefficients. The fermion number operator is defined 
as 

Np = ^""^ O'jrrfljm = { a jm a jm 4" a >j— m a 'j~iTi) ■ (35.3) 



m m>0 

As we have already seen, the J = pair creation and annihilation operators satisfy the commutation 
relation 

[5,5t] = l-^, (35.4) 

while the pairing Hamiltonian is 

H = -GQS^S. (35.5) 

The seniority Vp is defined as the number of fermions not coupled to J = 0. If only pairs of J = are 
present (i.e. Vp = 0), the eigenvalues of the Hamiltonian are (as already seen in eq. (33.16)) 

E{NF ,v F = 0) = -an{f^-§). (35,, 

For non-zero seniority the eigenvalues of the Hamiltonian are 

E(Np, V F ) =-j(N F - V F )(2n -Np-Vp + 2). (35.7) 

We denote the operators Np, Vp and their eigenvalues by the same symbol for simplicity. 

In subsec. 33.4 it has been proved that the behaviour of the J = pairs can be described, up to 
first order corrections, in terms of Q-bosons. In particular, making the mapping 

S j S -> 6, Np -> 2iV, (35.8) 

the relevant pairing Hamiltonian of eq. (35.5) becomes 

H{N, V = 0) = -GQtfb = -GQ[N] Q , (35.9) 

which coincides with eq. (35.6) up to first order corrections in the small parameter, which is identified as 
T = — 2/0. Furthermore, the Q-bosons satisfy the commutation relation of eq. (33.24), which coincides 



with eq. (35.4) up to first order corrections in the small parameter, which is, consistently with the above 
finding, identified as T = —2/Q. Therefore the fermion pairs of J = can be approximately described 
as Q-bosons, which correctly reproduce both the pairing energies and the commutation relations up to 
first order corrections in the small parameter. 

For the case of nonzero seniority, one observes that eq. (35.7) can be written as 

i.e. it can be separated into two parts, formally identical to each other. Since the second part (which 
corresponds to the J = pairs) can be adequately described by the Q-bosons b, b\ and their number 
operator N, as we have already seen, it is reasonable to assume that the first part can also be described 
in terms of some Q-bosons d, d\ and their number operator V (with Vp — > 2V), satisfying commutation 
relations similar to eqs (11.2) and (11.3): 

[V, S] = d\ [V, d] = -d, dS - Q$d = 1. (35.11) 

From the physical point of view this description means that a set of Q-bosons is used for the J = 
pairs and another set for the J ^ pairs. The latter is reasonable, since in the context of this theory 
the angular momentum value of the J ^ pairs is not used explicitly. The J ^ pairs are just counted 
separately from the J = pairs. A Hamiltonian giving the same spectrum as in eq. (35.10), up to first 
order corrections in the small parameter, can then be written as 

H(N,V)=Gn([V] Q -[N} Q ). (35.12) 

Using eq. (33.27) it is easy to see that this expression agrees to eq. (35.10) up to first order corrections 
in the small parameter T = —2/Q. 

Two comments concerning eq. (35.12) are in place: 

i) In the classical theory states of maximum seniority (i.e. states with N = V) have zero energy. 
This is also holding for the Hamiltonian of eq. (35.12) to all orders in the deformation parameter. 

ii) A landmark of the classical theory is that E(N, V) — E(N, V = 0) is independent of N. This 
also holds for eq. (35.12) to all orders in the deformation parameter. 

Knowing the Schwinger realization of the su g (2) algebra in terms of g-bosons (sec. 15), one may 
wonder if the operators used here close an algebra. It is easy to see that_the operators b + d, d + b and 
N — V do not close an algebra. Considering, however, the operators (see S with p — 1) 

J + = b^Q' v/2 d, J_ = d^Q- v/2 b, J = -(N - V), (35.13) 

one can easily see that they satisfy the commutation relations 

[Jo, J±] = ±J±, J+J- - Q- 1 J_ J + = [2J ] Q . (35.14) 

Using the transformation 

J = J , J + = Q(l/2)(Jo-l/2)j +j J = J_Q(l/2)(J -l/2) ) (35 15) 

one goes to the usual su g (2) commutation relations 

[J , J±] = ±4, [J+, J-\ = [2 J ], (35.16) 

where q 2 = Q. One can thus consider eq. (35.14) as a rewriting of the algebra su ? (2), suitable for boson 
realization in terms of Q-bosons. 

It is clear that N + V is the first order Casimir operator of the uq(2) algebra formed above (since it 
commutes with all the generators given in eq. (35.13)), while N — V is the first order Casimir operator 
of its uq(1) subalgebra, which is generated by Jo alone. Therefore the Hamiltonian of eq. (35.12) can 
be expressed in terms of the Casimir operators of the algebras appearing in the chain uq(2) D uq(1) as 



E(N, V) = GQ( 



C 1 (u Q (2))-C 1 (u Q (l)) 



Q 



C 1 (u Q (2)) + C 1 (u Q (l)) 



), (35.17) 



i.e. the Hamiltonian has a uq(2) D uq(1) dynamical symmetry. 



35.2. Comparison to experiment 

In the construction given above we have shown that Q-bosons can be used for the approximate 
description of correlated fermion pairs in a single-j shell. The results obtained in the Q-formalism agree 
to the classical (non-deformed) results up to first order corrections in the small parameter. However, 
the Q-formalism contains in addition higher order terms. The question is then born if these additional 
terms are useful or not. For answering this question, the simplest comparison with experimental data 
which can be made concerns the classic example of the neutron pair separation energies of the Sn 
isotopes, used by Talmi EEUM 

In Talmi's formulation of the pairing theory, the energy of the states with zero seniority is given by 

E{N) d = NV + N ( N -^ A, (35.18) 

where N is the number of fermion pairs and Vq, A are constants. We remark that this expression is the 
same as the one in eq. (35.6), with the identifications 

A/(2Vb) = A = 2G, N F = 2N. (35.19) 

The neutron pair separation energies are given by 

AE(N + l) cl = E(N + l) d - E(N) c i = V (l + ^-iv) . (35.20) 

Thus the neutron pair separation energies are expected to decrease linearly with increasing N. (Notice 
from eq. (35.19) that A/Vn < 0, since Q > 0.) A similar linear decrease is predicted also by the 
Interacting Boson Model £53. 

In our formalism the neutron pair separation energies are given by 

AE(N + 1) Q = -GQ{[N + 1} Q - [N] Q ) = -GQQ N = -GQe TN . (35.21) 

Since, as we have seen, T is expected to be —2/ ft, i.e. negative and small, the neutron pair separation 
energies are expected to fall exponentially with increasing N, but the small value of T can bring this 
exponential fall very close to a linear one. 

The neutron pair separation energies of the even Sn isotopes from 104 Sn to 130 Sri (i.e. across the whole 
sdg neutron shell) have been fitted in E 3 ^ using both theories. Furthermore, in a fit of the logarithms 
of the energies has been performed, since eq. (35.21) predicts a linear decrease of the logarithm of 
the energies with increasing N. Both fits give almost identical results. Eq. (35.21) (in which the free 
parameters are GQ and T), gives a better result than eq. (35.20) (in which the free parameters are Vq 
and A/Vo) for every single isotope, without introducing any additional parameter, indicating that the 
higher order terms can be useful. 

One should, however, remark that 116 Sn lies in the middle of the sdg neutron shell. Fitting the 
isotopes in the lower half of the shell ( 104 Sn to 116 Sn) and the isotopes in the upper half of the shell 
( 118 Sn to 130 Sn) separately, one finds that both theories give indistinguishably good results in both 
regions. Therefore Q-deformation can be understood as expressing higher order correlations which 
manifest themselves in the form of particle-hole asymmetry. It is also known that a strong subshell 
closure exists at N=64 (which corresponds to 114 Sn). The presence of this subshell closure can also 
affect the neutron pair separation energies in a way similar to the one shown by the data. 

In I 3 ! a fit of the neutron pair separation energies of the Pb isotopes from 186 Pb to 202 Pb has also 
been attempted. In this case both theories give indistinguishably good fits. This result is in agreement 
with the Sn findings, since all of these Pb isotopes lie in the upper half of the pfh neutron shell. 
Unfortunately, no neutron pair separation energy data exist for Pb isotopes in the lower part of the pfh 
neutron shell. 

Concerning the values of T obtained in the case of the Sn isotopes (T = —0.0454, T = —0.0447), one 
observes that they are slightly smaller than the value (T = —0.0488) which would have been obtained 



by considering the neutrons up to the end of the sdg shell as lying in a single-j shell. This is, of course, 
a very gross approximation which should not be taken too seriously, since it ignores the fact that most 
properties of nuclei can be well accounted for by the valence nucleons alone, without being affected by 
the closed core. In the case of the Pb isotopes mentioned above, however, the best fit was obtained 
with T = —0.0276, which is again slightly smaller than the value of T = —0.0317 which corresponds to 
considering all the neutrons up to the end of the pfh shell as lying in a single-j shell. 

In summary, we have shown that pairing in a single-j shell can be described, up to first order 
corrections, by two Q-oscillators, one describing the J = pairs and the other corresponding to the 
J / pairs, the deformation parameter T = In Q being related to the inverse of the size of the shell. 
These two oscillators can be used for forming an suq(2) algebra. A Hamiltonian giving the correct 
pairing energies up to first order corrections in the small parameter can be written in terms of the 
Casimir operators of the algebras appearing in the uq(2) D uq(1) chain, thus exhibiting a quantum 
algebraic dynamical symmetry. The additional terms introduced by the Q-oscillators serve in improving 
the description of the neutron pair separation energies of the Sn isotopes, with no extra parameter 
introduced. 

In the previous section a generalized deformed oscillator describing the correlated fermion pairs of 
J = exactly has been introduced. This generalized deformed oscillator is the same as the one giving 
the same spectrum as the Morse potential (sec. 37), up to a shift in the energy spectrum. The use 
of two generalized deformed oscillators for the description of J = pairs and J ^ pairs in a way 
similar to the one of the present section is a straightforward task, while the construction out of them of 
a closed algebra analogous to the suq(2) obtained here is an open problem. The extension of the ideas 
presented here to the case of the BCS theory is an interesting open problem. 



35.3. Other approaches 



A g-deformed version of the pairing theory was assumed by Petrova t^j and Shelly Sharma 
with satisfactory results when compared to experimental data. The present construction offers some 
justification for this assumption, since in both cases the basic ingredient is the modification of eq. 
(35.4). It should be noticed, however, that the deformed version of eq. (35.4) considered in EEEHH is 
different from the one obtained here (eq. (33.24)). A basic difference is that the deformed theory of 
©HI reduces to the classical theory for q — > 1, so that g-deformation is introduced in order to describe 
additional correlations, while in the present formalism the Q" oscu lato r s involved for Q — > 1 reduce to 
usual harmonic oscillators, so that Q-deformation is introduced in order to attach to the oscillators the 
anharmonicity needed by the energy expression (eq. (35.6)). 

Continuing along the same line Shelly Sharma and Sharma derived Random Phase Approximation 
(RPA) equations for the pairing vibrations of nuclei differing by two nucleons in comparison to the 
initial one and applied their method to the study of the + states of the Pb isotopes, which offer a 
good example of pairing vibrations in nonsuperconducting nuclei. Furthermore, using deformed quasi- 
particle pairs coupled to zero angular momentum they developed a deformed version of the quasi-boson 
approximation for + states in superconducting nuclei and tested it against a schematic two-level shell 



ras been developed by Avancini and Menezes 
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model E^3. Another deformed two-level shell model 
g-deformed boson mappings have been developed in 

In addition a g- deformed version of the many-body BCS approximation for a pure pairing force has 
been developed 1^3, using g-deformed fermions satisfying su g (n) anticommutation relations. A set of 



quantum BCS equations and a g-analog of the gap equation have been derived 



36. Anisotropic quantum harmonic oscillators with rational ratios of frequencies 

3-dim anisotropic harmonic oscillators with rational ratios of frequencies (RHOs) are of current 
interest beca use of their relevance as possible underlying symmetries of superdeformed and hyperde- 
formed nuclei l 342 tr 43 l. In particular, it is thought that superdeformed nuclei correspond to a ratio 

of frequencies of 2:1. while hyperdeformed nuclei correspond to a 3:1 ratio. In addition they have been 
recently connected 0@ to the underlying geometrical structure in the Bloch-Brink a-cluster model 
and possibly to the interpretation of the observed shell structure in atomic clusters (see sec. 38), 



especially after the realization that large deformations can occur in such systems The 2-dim RHO 
is also of interest, since its single particle level spectrum characterizes the underlying symmetry of 
"pancake" nuclei 



RHOs are examples of maximally superintegrable systems 



us 



in N dimensions. Superintegrable 



The two-dim 



systems in N dimensions have more than N independent integrals (constants of motion). Maximally 
superintegrable s\ 



anisotropic harmonic oscilla- 
assical and the quantum mechanical 
have also been considered. While at 
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tors have been the subject of several investigatio ns^ Ja q tb at the_c 
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level. The special cases with frequency ratios 1:2 l 143 H 364 i and 1:3 
the classical level it is clear that the su(N) or sp(2N,R) algebras can be used for the description of the 
N-dimensional anisotropic oscillator, the situation at the quantum level, even in the two-dimensional 
case, is not as simple. 

In this section we are going to prove that a generalized deformed u(2) algebra is the symmetry algebra 
of the two-dimensional anisotropic quantum harmonic oscillator, which is the oscillator describing the 
single-particle level spectrum of "pancake" nuclei, i.e. of triaxially deformed nuclei with u x » u y , uj z 
[El. The method can be extended to the 3-dim RHO in a rather straightforward way. 

36.1. A deformed u(2) algebra 

Let us consider the system described by the Hamiltonian: 



H = - \Px +Py + — + — 

2 \ m l n l 



(36.1) 



where m and n are two natural numbers mutually prime ones, i.e. their great common divisor is 
gcd(m, n) = 1. 

We define the creation and annihilation operators e23 



+ x/m — ip T 
V2 

t _ y/n-ip y 
These operators satisfy the commutation relations: 



x/m + ip x 
V2 

y/n + ip y 
V2 ■ 



a, a' 



m 



— , other commutators = 0. 
n 



Further defining 



U=±{a,<J}, W=^{b,tf}, 
one can consider the enveloping algebra generated by the operators: 



S = ~(U -W), H = U + W. 
These genarators satisfy the following relations: 

[S Q ,S±] = ±S ± , [H,3i] = Q, for i = 0,±, 

**-n(-^)n(-^). 



(36.2) 
(36.3) 

(36.4) 
(36.5) 

(36.6) 
(36.7) 

(36.8) 

(36.9) 
(36.10) 



The fact that the operators Si, % = 0, ± are integrals of motion has been already realized in E^l. 

The above relations mean that the harmonic oscillator of eq. (36.1) is described by the enveloping 
algebra of the generalization of the u(2) algebra formed by the generators S , S + , £_ and H, satisfying 
the commutation relations of eq. (36.8) and 

[£_, 3 + ] = F m>n {H, S + l)- F m , n (H, S Q ), (36.11) 

where 

171 ( 2k — 1 \ n ( 2P — 1 \ 

F m , n (H, S ) = J] [H/2 + S - -^-J n [H/2 - S + J . (36.12) 

In the case of m = 1, n = 1 this algebra is the usual u(2) algebra, and the operators S , S± satisfy the 
commutation relations of the ordinary u(2) algebra, since in this case one easily finds that 

[S-,S+] = -2S . (36.13) 

In the rest of the cases, the algebra is a deformed version of u(2), in which the commutator [S-, S + ] is 
a polynomial of 5*0 of order m + n — 1. In the case with m — 1, n = 2 one has 

[S-, S + ] = 3S 2 -HS -^ + ^-, (36.14) 

4 16 

i.e. a polynomial quadratic in So occurs, while in the case of m = 1, n = 3 one finds 

S + ] = -4S 3 + 3HS 2 - 7 -S -^ + j, (36.15) 
i.e. a polynomial cubic in So is obtained. 

36.2. The representations 

The finite dimensional representation modules of this algebra can be found u sing th e concept of the 
generalized deformed oscillator (sec. 12), in a method similar to the one used inl2£3'L£3 for the study of 
quantum superintegrable systems. The operators: 

A* — S + , A = S-, J\f = S — u, u— constant, (36.16) 

where u is a constant to be determined, are the generators of a deformed oscillator algebra: 

~N,A ] ]=A\ [M,A\ = -A, A ] A = <$>(H,N), AA^ = ®{H, Af+ 1). (36.17) 

The structure function $ of this algebra is determined by the function F m n in eq. (36.12): 

®(H, AT) = F mtn (H, AT + u) = J[ [H/2 +Af + u- J J] [H/2 -M-u + -^—j . (36. 18 

The deformed oscillator corresponding to the structure function of eq. (36.18) has an energy dependent 
Fock space of dimension N + 1 if 

$(£,0) = 0, <$>{E,N +1) = 0, $(E,k)>0, for jfe = 1, 2, . . . , N. (36.19) 

The Fock space is defined by: 

H\E,k>= E\E,k>, N\E,k>=k\E,k>, a\E,0>= 0, (36.20) 



A j \E,k>= yJ$(E,k + l)\E,k + l >, A\E,k >= yJ$(E,k)\E,k- 1 > . (36.21) 
The basis of the Fock space is given by: 

\E,k >= (A ] ) k \E,Q >, k = 0,l,...N, (36.22) 



where the "factorial" [A;]! is defined by the recurrence relation: 

[0]! = 1, [k]\ = $(E,k)[k-l]\ . 



(36.23) 



Using the Fock basis we can find the matrix representation of the deformed oscillator and then the 
matrix representation of the algebra of eqs (36.8), (36.12). The solution of eqs (36.19) implies the 
following pairs of permitted values for the energy eigenvalue E and the constant u: 



where p = 1,2, ... ,m, 



„ AT 2p-l 2q-l 
E = N + — + — . 

2m 2n 

q = 1,2, ... ,n, and 

1 (2p- 1 2q-l 



u 



2m 



2n 



N 



(36.24) 



(36.25) 



the corresponding structure function being given by: 

2p-l 2k — 1 



x 



n *+ 



k=l 



2m 



2m 



^rf.r 2g-l 2£-l 
11 [N-x+-^ + 



i=i 



2n 



2n 



r(mx+p) r ((N — x)n + q + n) 



(36.26) 

m m n n T (mi + p — m) Y ((JV — x)n + q) 

where T(x) denotes the usual Gamma-function. In all these equations one has N — 0, 1, 2, . . ., while the 
dimensionality of the representation is given by N + 1. Eq. (36.24) means that there are m ■ n energy 
eigenvalues corresponding to each iV value, each eigenvalue having degeneracy N + 1. (Later we shall 
see that the degenerate states corresponding to the same eigenvalue can be labelled by an "angular 
momentum".) 

It is useful to show at this point that a few special cases are in agreement with results already 
existing in the literature. 

i) In the case m = 1, n = 1 eq. (36.26) gives 



while eq. (36.24) gives 



in agreement with Sec. IV. A of 



$(E,x) = x(N + 1-x) 
E = N + 1, 



(36.27) 
(36.28) 



ii) In the case m = 1, n = 2 one obtains for q = 2 



$(E, x) = x(N + l-x)(N+--x 



E = N + 



while for q — 1 one has 



$(£, x) = x(N + l-x)(N+--x 



E = N+~. 

4 



(36.29) 



(36.30) 



These are in agreement with the results obtained in Sec. IV. F of E^l for the Holt potential (for 5 = 0). 
iii) In the case m = 1, n = 3 one has for q = 1 



®{E, x) = x{N + 1- x)(n + -- x^ (n + - 
while for q = 2 one obtains 



E = N + 



${E, x) = x(N + l-x)(N+--x) [N + 



x 



E = N+1, 



(36.31) 



(36.32) 



and for q = 3 one gets 



${E, x) = x(N + 1 - x) [N + 



iV H x 

3 



E = N + 



(36.33) 



These are in agreement with the results obtained in Sec. IV.D of for the Fokas-Lagerstrom potential. 

In all of the above cases we remark that the structure function has forms corresponding to various 
versions of the generalized deformed parafermionic algebra of eq. (18.1), the relevant conditions of eq. 
(18.2) being satisfied in all cases. It is easy to see that the obtained algebra corresponds to this of the 
generalized parafermionic oscillator in all cases with frequency ratios 1 : n. 

The energy formula can be corroborated by using the corresponding Schrodinger equation. For the 
Hamiltonian of eq. (36.1) the eigenvalues of the Schrodinger equation are given by: 



E 



1 ( 1 

m V 2 



+ 



1 



n 



where n x = 0, 1, . . . and n y = 0,1, . . .. Comparing eqs (36.24) and (36.34) one concludes that: 

iV = [n x /m\ + [n y /n] , 
where [x] is the integer part of the number x, and 

p = mod(n x , m) + 1, q = mod(n y , n) + 1. 



(36.34) 



(36.35) 



(36.36) 



The eigenvectors of the Hamiltonian can be parametrized by the dimensionality of the representation 
N , the numbers p, q, and the number k = 0, 1, . . . , N. k can be identified as [n x /m]. One then has: 



H 



Sf) 



N 

(p,q) 

s 4 



N 

(p,q) 



k + 



Ar 2p-l 2q-l 

N + — h — 

2m 2n 

2p-l 2q-l 



N 

(p,q) 



2m 



2n 



N k x 

(p,q) ' / 

N 



-N 



N 



N 

(p,q) 



k + 1 



(p,q) ' 
k-i 



N 

(p,q) 



(36.37) 

(36.38) 
(36.39) 
(36.40) 



36.3. The "angular momentum" quantum number 



It is worth noticing that the operators So, S± do not correspond to a generalization of the angu- 



lar momentum, Sq being the operator corresponding to the Fradkin operator S, 
corresponding "angular momentum" is defined by: 

L = -i (S+ - S-) . 

The "angular momentum" operator commutes with the Hamiltonian: 

[H,L } = 0. 



q 128 
xx ^yy — 
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The 



(36.41) 



(36.42) 



Let \£ > be the eigenvector of the operator L corresponding to the eigenvalue i. The general form of 
this eigenvector can be given by: 



N i k ci. 



N 

(p,q) 



k 



(36.43) 



In order to find the eigenvalues of Lq and the coefficients c k we use the Lanczos algorithm E^3, as 
formulated in EH. From eqs (36.39) and (36.40) we find 



L \i >-- 



N % k c 

>=eJ2 



k=0 



N 

(p,q) 



v/[fc]! 

N 

E 



N 

(p,q) 



k=l 



N 



k 



(36.44) 



From this equation we find that: 

c k = (-l) k 2- k / 2 H k (t/V2)/M, U 2 = Y, 2~ n H 2 n (£/V2)/n\ 



N 

E 

n=0 



where the function H k (x) is a generalization of the "Hermite" polynomials (see also l 82 U 
the recurrence relations: 

if_i(a;) = 0, H (x) = 1, 
H k+1 (x) = 2xH k (x) - 2^ t>q) {k)H k - 1 {x), 
and the "angular momentum" eigenvalues i are the roots of the polynomial equation: 

H N+1 (£/V2) = 0. 



(36.45) 

), satisfying 

(36.46) 
(36.47) 



(36.48) 



Therefore for a given value of N there are N + 1 "angular momentum" eigenvalues £, symmetric around 
zero (i.e. if i is an "angular momentum" eigenvalue, then — i is also an "angular momentum" eigenvalue). 
In the case of the symmetric harmonic oscillator (m/n = 1/1) these eigenvalues are uniformly distributed 
and differ by 2. In the general case the "angular momentum" eigenvalues are non-uniformly distributed. 
For small values of N analytical formulae for the "angular momentum" eigenvalues can be found E3. 
Remember that to each value of N correspond m ■ n energy levels, each with degeneracy N + 1. 

In order to have a formalism corresponding to the one of the isotropic oscillator, let us introduce 
for every iV and (m, n, p, q) an ordering of the "angular momentum" eigenvalues 



pL, m ,n,p, q ^ where L = N and fi= -L,-L + 2,...,L-2,L, 



by assuming that: 



t L m,n,p, q < t L,m,n, P , q [f ^ < ^ 



the corresponding eigenstates being given by: 



\L,fi;m,n,p,q) = J2 ' 

fc=0 



) k H k (£^ m ' n 'P'^/V2) 



Mj2 k l 2 [k}\ 



N 



L 

E 

fc=0 



N 

(p,q) 



(36.49) 
(36.50) 

(36.51) 



The above vector elements constitute the analogue corresponding to the basis of "spherical harmonic" 
functions of the usual oscillator. The calculation of the "angular momentum" eigenvalues of eq. (36.49) 
and the coefficients di, d 2 , ■ ■ ■ , g?l+i in the expansion of eq. (36.51) is a quite difficult task. The existence 
of general analytic expressions for these quantities is not obvious. The first few "angular momentum" 
eigenvalues are given by: 



nl,m,n,p,q _i_ 

t ±1 — ± 



and 



1 r(m + p) r(n + q) 
\\ m m n n r(p) T(q) ' 

p2,m,n,p,q _ g 



T(m + p) T(2n + q) T(2m + p) T(n + q) 



\m m n n \ T(p) T(n + q) T(m + p) T(g) 



(36.52) 

(36.53) 
(36.54) 



For L > 2 the analytic expressions of the angular momentum eigenvalues and the coefficients dk are 
longer, but their calculation is a straightforward task. Numerical results for these quantities in the cases 
of frequency ratios 1:2 and 1:3 are given in_^'Eli. 

for details) one finds out the following points: 



. , i.e. "two copies" of the u(2) 
i.e. "three copies" of the u(2) 



After working out a few examples (see F°lf 71 lf 72 l 

i) In the basis described by eqs. (36.16)-(36.19) it is a trivial matter to distinguish the states 
belonging to the same irrep for any m : n ratio, while in the Cartesian basis this is true only in the 1:1 
case. 

ii) In the 1:2 case the irreps have degeneracies 1, 1, 2, 2, 3, 3, 4, 4, 
degeneracies 1, 2, 3, 4, ... are obtained. 

iii) In the 1:3 cases the degeneracies are 1, 1, 1, 2, 2, 2, 3, 3, 3, . . 
degeneracies are obtained. 

iv) It can be easily seen that the l:n case corresponds to "n copies" of the u(2) degeneracies. 

v) Cases with both m, n different from unity show more complicated degeneracy patterns, also 
correctly reproduced by the above formalism. In the 2:3 case, for example, the degeneracy pattern is 1, 
1,1,1,1,2,1,2,2,2,2,3,2,3,3,.... 

vi) The only requirement for each energy eigenvalue to correspond to one irrep of the algebra is that 
m and n have to be mutually prime numbers. If m and n possess a common divisor other than 1, then 
some energy eigenvalues will correspond to sums of ir rep s. i.e. to reducible representations. 

vii) The difference between the formalism used in BUSS and the one used here is that in the 
former case for given m and n appropriate operators have to be introduced separately for each set of 
(p, q) values, while in the present case only one set of operators is introduced. 



36.4- Multisections of the isotropic oscillator 

In III the concept of bisection of an isotropic harmonic oscillator has been introduced. One can 
easily see that multisections (trisections, tetrasections, . . . ) can be introduced in a similar way 0. 
The degeneracies of the various anisotropic oscillators can then be obtained from these of the isotropic 
oscillator by using appropriate multisections. 

Using the Cartesian notation (n x ,n y ) for the states of the isotropic harmonic oscillator we have the 
following list: 

N=0: (00) 

N=l: (10) (01) 

N=2: (20) (02) (11) 

N=3: (30) (03) (21) (12) 

N=4: (40) (04) (31) (13) (22) 

N=5: (50) (05) (41) (14) (32) (23), 
where N = n x + n y . The corresponding degeneracies are 1, 2, 3, 4, 5, 6, ... , i.e. these of u(2). 

A bisection can be made by choosing only the states with n y =even. Then the following list is 



obtained: 


N=0 


(00) 


N=l 


(10) 


N=2 


(20) (02) 


N=3 


(30) (12) 


N=4 


(40) (04) (22) 


N=5 


(50) (14) (32) 



The degeneracies are 1, 1, 2, 2, 3, 3, ... , i.e. these of the anisotropic oscillator with ratio of frequen- 
cies 1:2. The same degeneracies are obtained by choosing the states with n. y =odd. Therefore a bisection 
of the isotropic oscillator, distinguishing states with mod(n y , 2) = and states with mod(ra 2/ ,2) = 1, 
results in two interleaving sets of levels of the 1:2 oscillator. 

By analogy, a trisection can be made by distinguishing states with mod(n y , 3) = 0, or mod(n y , 3) = 
1, or mod(n y ,3) = 2. One can easily see that in this case three interleaving sets of states of the 1:3 
oscillator, having degeneracies 1, 1, 1, 2, 2, 2, 3, 3, 3, ... , occur. 

Similarly a tetrasection can be made by distinguishing states with mod(n y , 4) = 0, or mod(n y , 4) = 
1, or mod(n y , 4) = 2, or mod(n y , 4) = 3. The result is four interleaving sets of states of the 1:4 oscillator, 



having degeneracies 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, ... . 

By bisecting n x and trisecting n y one is left with six interleaving sets of states with degeneracies 1, 
1, 1, 1, 1, 2, 1, 2, 2, 2, 2, 3, 2, 3, 3, ... , i.e. degeneracies of the 2:3 oscillator. 

By bisecting (or trisecting, tetresecting, etc) both n x and n y one is obtaining the original u(2) 
degeneracies of the isotropic oscillator. 

It is therefore clear that the degeneracies of all m : n oscillators can be obtained from these of the 
isotropic oscillator by appropriate multisections. In particular: 

i) The degeneracies of the 1 : n oscillator can be obtained from these of the 1:1 (isotropic) oscillator 
by n-secting n y or n x . 

ii) The degeneracies of the m : n oscillator can be obtained from these of the 1:1 oscillator by 
m-secting n x and n-secting n y . 

36.5. Connection to 

For the special case m — 1, n — 2 it should be noticed that the deformed algebra received here 



coincides with the finite W algebra W3 P 75 ^ 76 ^ 77 ! ] 378 ! , The commutation relations of the W3 algebra 
are 

[H w , E w ] = 2E W , [H w , F w ] = -2F W , [E w , F w ] = + C w , (36.55) 

[C w , Ew] = [C w , F w ] = [C w , H w ] = 0, (36.56) 
while in the m — 1, n — 2 case one has the relations 

[Af, A ] ] = A\ [Af, A] = -A, [A, A ] ] = 35 2 - ^ - HS + (36.57) 
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[H, A ] ] = [H, A] = [H, S ] = 0, (36.58) 

with So = jV + u (where u a constant). It is easy to see that the two sets of commutation relations are 
equivalent by making the identifications 

F w = oA\ E w = pA, H w = -2S + kH, C w = f(H), (36.59) 

with 

P° = li k= 1 -, f( H ) = -~H* + ~. (36.60) 



36.6. Discussion 

In conclusion, the two-dimensional anisotropic quantum harmonic oscillator with rational ratio of 
frequencies equal to m/n, is described dynamically by a deformed version of the u(2) Lie algebra, the 
order of this algebra being m + n — 1. The representation modules of this algebra can be generated 
by using the deformed oscillator algebra. The energy eigenvalues are calculated by the requirement of 
the existence of finite dimensional representation modules. An "angular momentum" operator useful 
for labelling degenerate states has also been constructed. The algebras obtained in the special cases 
with 1 : n ratios are shown to correspond to generalized parafermionic oscillators. In the special case 
of m : n = 1 : 2 the resulting algebra has been identified as the finite W algebra W3 . Finally, it is 
demonstrated how the degeneracies of the various m : n oscillators can be obtained from these of the 
isotropic oscillator by appropriate multisections. 

The extension of the present met hod to the three-dimensional anisotropic quantum harmonic os- 
cillator is already receiving attention r 79 lr 80 l, since it is of clear interest in the study of the symmetries 
underlying the structure of superdeformed and hyperdeformed nuclei EE1IE1. 

37. The use of quantum algebras in molecular structure 



The techniques developed in this article can be applied in very similar ways in describing properties 
of diatomic and polytomic molecules. A brief list will be given here. 



1) Rotational spectra of diatomic molecules have been described in terms of the su 9 (2) model 



147,381,382 



3N l 



As in the case of nuclei, g is a phase factor (g = e tT ). In molecules r is of the 
order of 0.01. The use of the su g (2) symmetry leads to a partial summation of the Dunham expan- 
sion describing the rotational-vibrational spectra of diatomic molecules Molecular backbending 



(bandcrossing) has also been described in this framework t£a Rotational spectra of symmetric top 



molecules have also been considered E2J>t^a in the framework of the su g (2) sy mme try. Furthermore, two 



g-deformed rotators with slightly different parameter values have been used t ^l fo r the description of 
AI = 1 staggering effects in rotational bands of diato^ 
of AI = 2 staggering effects in diatomic molecules see 
nuclei see sec. 22.) 

2) Vibrational spectra of diatomic molecules have been d escribed in terms of g-deformed anharmonic 



molecules (see also 1223). (For a discussion 
0J ~'" For a discussion of staggering effects in 



oscillators having the su g (l,l) 
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mmetry, as well as in terms of generalized 
These results, combined with 1), lead to 



„£ 26 1392IJ393IJ394 

Hi^l. A two-parameter deformed anharmonic oscillator 



lor the u g (2) D o,(2) 
deformed oscillators similar to the ones used in 
the full summation of the Dunham expansion 
with u gp (2) D o qp (2) symmetry has also been considered l 3 -^. 

3) The physical content of the anharmonic oscillators mentioned in 2) has been clarified by con- 
structing WKB equivalent potentials (WKB-EPs) lIS and classical equivalent potentials 0, similar to 
the ones of sec. 13. The results have been corroborated by the study of the relation between su g (l,l) 
and the anharmonic oscillator with x 4 anharmini cities 1^1. The WKB-EP corresponding to the suJl.l) 
anharmonic oscillator has been connected to a class of Quasi-Exactly Soluble Potentials ( QES Ps) e£3. 

4) Generalized deformed oscillators giving the same spectrum as the Morse potential eHj and the 
modified Poschl-Teller potential ^2§, as well as a deformed oscillator containing them as special cases 

11 ^ have also been constructed. In addition, q-deformed versions of the Morse potential have been 



given, either by using the 
the usual Morse potential 



,(2,1) symmetry 



or by solving a g-deformed Schrodinger equation for 
£3. For the sake of completeness it should be mentioned that a deformed 
oscillator giving the same spectrum as the Coulomb potential has also bee n cons tructed 0. 

^hdeformed version of the vibron model for diatomic molecules E°J3e2a has been constructed 
in a way similar to that described in sec. 31. 
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6) For vibrational spectra of polyatomic molecules a model of n coup led generalized deformed oscil- 



lators has been built containg the approach of Iachello and Oss Elllell as a special case. Furthermore 
a system of two Q-defor med oscillators coupled so that the total Hamiltonian has the suq(2) symmetry 
(in the way described in c!3) has been proved to be appropriate for the description of vibrational spec- 
tra of triatomic molecules eU. In addition, a 3-dimensional anisotropic g-deformed harmonic oscillator 
has been used cl3 for the description of vibrational spectra of triatomic molecules. A description of a 
tetrahedral molecule (methane) in terms of g-deformed oscillators has been given in 0. 

7) Quasi-molecular resonances in the systems 12 C+ 12 C and 12 C+ 16 have been described in terms 



of a g-deformed oscillator plus a rigid rotator cil 



A review of several of the above topics, acco mp anied by a detailed and self-contained introduction to 



quantum algebras, has been given by Ravchev c±3. A recent review on algebraic methods in molecular 



spectroscopy has been given by Kellman I 



applied in molecular spectroscopy can be found in the books by Iachello and Levine 
and Van Isacker 0. 



while more extensive presentations of algebraic techniques 

123 and by Frank 



38. The use of quantum algebras in the structure of atomic clusters 



Metal clusters have been recently the subject of many investigations (see p ^ 21 }! 422 ! for relevant 
reviews). One of the first fascinating findings in their study was the appearance of magic numbers 



423 ,[424 , 42"4|426 



427 428 429 



, analogous to but different from the magic numbers appearing in the shell structure 
of atomic nuclei EH This analogy led to the early description of metal clusters in terms of the Nilsson- 
Clemenger model E 3 !, which is a simplified version of the Nilsson model [||]|Eij'IZi]'izl of atomic nuclei, 



in which no spin-orbit interaction is included. Further theoretical investigations in terms of the jellium 



model l 432 " 433 ! demonstrated that the mean field potential in the case of simple metal clusters bears great 
re Woods-Saxon potential of atomic nuclei, with a slight modification of the "wine 
The Woods-Saxon potential itself looks like a harmonic oscillator truncated at a 



similarities to__t 
bottle" type 



13 1 



certain energy value and flattened at the bottom. It should also be recalled that an early schematic 
explanation of the magic numbers of metallic clusters has been given in terms of a scheme intermediate 
between the level scheme of the 3-dimensional harmonic oscillator and the square well EH Again in this 
case the intermediate potential resembles a harmonic oscillator flattened at the bottom. 

On the other hand, as we have seen in sec. 28, a g-deformed version of the 3-dimensional harmonic 
oscillator has been constructed Ell, taking advantage of the u ? (3) ~) ap g (3) symmetry. The spectrum of 
this 3-dimensional g-deformed harmonic oscillator has been fo iind El 3 to reproduce very well the spec- 
trum of the modified harmonic oscillator introduced by Nilsson t^>EZ3, without the spin-orbit interaction 
term. Since the Nilsson model without the spin orbit term is essentially the Nilsson-Clemenger model 
used for the description of metallic clusters it is worth examining if the 3-dimensional g-deformed 
harmonic oscillator can reproduce the magic numbers of simple metallic clusters. 

It has indeed been found e23 that the 3-dimensional g-deformed harmonic oscillator with u g (3) D 
sOq(3) symmetry correctly predicts all experimentally observed magic numbers of alkali metal clusters 

§to 1500, which is the expected limit of validity for theories based on the filling of electronic shells 
. This indicates that u g (3), which is a nonlinear deformation of the u(3) symmetry of the spherical 
(3-dimensional isotropic) harmonic oscillator, is a good candidate for being the symmetry of systems of 
alkali metal clusters. 



39. Outlook 

Nobody likes binding himself by statements concerning the future. However, we attempt to give 
here a partial list of open problems, roughly following the order of the material in this review: 

1) The list of physical systems which can be classified under a generalized deformed su(2) symmetry 
(sec. 17) or under a generalized deformed parafermionic oscillator scheme (sec. 18) can be enlarged. 
Self-similar potentials and isospectral oscillator Hamiltonian systems could probably be related to these 
symmetries. 

2) The description of B(E2) values in terms of the su g (2) model attempted in sec. 21 takes into 
account only the kinematical deformation effects. In order to take into account dynamical deformation 
effects, one has to build a larger algebra, of which the quadrupole operators will be members. The case 
of the u g (3)Dso g (3) algebra has been considered in sees 28, 29. A detailed study is needed in order 
to see if the experimental data support the modifications predicted by these models. In particular the 
su g (2) prediction about B(E2) values increasing with increasing angular momentum J, supported by the 
predictions of other models as well (see sec. 21), requires further testing against detailed experimental 
data. 

3) The construction of Clebsch-Gordan coefficients for the subclass of generalized deformed su(2) 
algebras for which this could be possible (see sees 17, 25) is an open problem. 

4) The sug(3) D so q (3) decomposition for su g (3) irreps other than the completely symmetric ones 
(see sees 27-29 for the current state of the art) remains an open problem, the solution of which is 
necessary for developing a deformed version of the su(3) limit of the Interacting Boson Model. 

5) Realizations of multi-level shell models in terms of deformed bosons (see sees 33-35 for some one- 
and two-level cases) should be further pursued. 

6) The symmetry algebras of the various 3-dim anisotropic harmonic oscillators with rational ratios 
of frequencies should be worked out, since they are of interest in relation to superdeformed and hyper- 
deformed nuclei, and possibly to deformed atomic clusters (see sees 36, 38 for references). A deformed 
u(3) algebra shou ld occur i n this case, which could serve as the basis for building a deformed analog of 
the Elliott model Ell §12 ill suitable for superdeformed nuclei. 



7) In molecular physics (sec. 37) the study of vibrations of highly symmetric polyatomic molecules 
(including fullerenes) by these techniques is of interest. 

8) In the structure of atomic clusters (sec. 38) the 3-dimensional g-deformed harmonic oscillator 
(sec. 28) appears to be an appropriate basis for further investigations. 

In addition the following final comments are also in place: 

a) The usefulness of any "new theory" in nuclear structure is usually judged by its ability to predict 
some novel excitation mode. For example, the prediction of the l ow-l ying (around 3 MeV) 1 + mode 
in even nuclei has been a success of the Interacting Boson Model-2 t^l. In the present case, the su g (2) 



model prediction that B(E2) values in deformed nuclei increase as a function of angular momentum, and 
do not reach a saturation value as the su(2) symmetry suggests, is an interesting point which requires 
further experimental testing, as already mentioned in sec. 21. 

b) Any "new theory" in nuclear structure oughts to acquire a microscopic justification, by establish- 
ing links to the underlying fermionic degrees of freedom. The mappings of sees 33 and 34 are a first step 
in this direction, although they deal only with fermion pairs of zero angular momentum in a single- j 
shell. The generalization of these mappings to pairs of nonzero angular momentum and to multi-j level 
schemes is a formidable but challenging task. 

c) Quantum algebras, being nonlinear extensions of Lie algebras, are specially suited for describing 
small perturbations in systems characterized by Lie symmetries. Several examples in nuclei, molecules, 
and atomic clusters have already been mentioned. Future work should clarify if the success of quantum 
algebras is limited within these bounds or if there is some deeper physical reason behind their success. 
For the time being, quantum algebras have been proved to be a useful tool for the description of small 
deviations from the usual Lie symmetries. 

d) An interesting question is if quantum algebras are not only suitable for describing small deviations 
from Lie symmetries, but in addition can bridge different Lie symmetries. Preliminary work in this 
direction has been reported in sec. 30. 
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